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Recently a variety of technologies have been proposed for improving surface properties of poly- 
mers. Among them is surface grafting of polymers. Although this is a rather new technology, 
polymer surface grafting offers versatile means for providing existing polymers with new 
functionalities such as hydrophilicity, adhesion, biocompatibility, conductivity, anti-fogging, 
anti-fouling, and lubrication. This review article describes various methods of grafting and 
grafted surface characterizations. Medical and non-medical applications connected with this 
polymer surface grafting are also presented referring to recent publications. 

Keywords: Surface modification. Graft polymerization. Graft chain. Adhesion, Surface charac- 
terization 
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1 

Introduction 

Most of the untreated surfaces of polymers used in industry are not hydrophilic 
but hydrophobic. It is, therefore, difficult to bond these nonpolar polymer sur- 
faces directly to other substances like adhesives, printing inks, and paints 
because they generally consist of polar compounds. On the other hand, polymer 
surfaces generally adsorb proteins when brought into direct contact with a bio- 
logical system, resulting in cell attachment or platelet aggregation. The protein 
adsorption and attachment of biological components trigger a subsequent 
series of mostly adverse biological reactions toward the polymeric materials. 
Therefore, the technologies for surface modification of polymers or regulation 
of the polymer surface interaction with other substances have been of prime 
importance in polymer applications from the advent of polymer industries. 
Some of the technologies have been directed to introduction of new function- 
alities onto polymer surfaces. The new functionalities introduced include 
improved surface hydrophilicity, hydrophobicity, biocompatibility, conductivi- 
ty, anti-fogging, anti-fouling, grazing, surface hardness, surface roughness, 
adhesion, lubrication, and antistatic property. Theoretically, there is a large dif- 
ference in properties between the surface and the bulk of a material and only 
the outermost surface is enough to be taken into consideration when the sur- 
face properties are concerned. However, this is not the case for polymer surfaces, 
as the physical structure of the outermost polymer surface is generally not fixed 
but continuously changing with time due to the microscopic Brownian motion 
of polymer segments. The polymer surface generally has high segmental mobil- 
ity even at room temperature, in contrast to the rigid surface of metals and 
ceramics. This propensity for motion of polymer segments suggests that a poly- 
mer surface cannot be described as a depthless, two-dimensional plane, but is 
in reality a region with some thickness. Therefore, we should look at the poly- 
mer surface in a direction not only parallel to the plane but also perpendicular 
to the surface plane in order to understand well the nature of polymer surface. 
The surface of copolymers and blended polymers is likely to be heterogeneous 
in the horizontal chemical composition, but these polymers have a distribution 
of chemical composition also in the vertical direction. This vertical distribution 
is called a depth profile. 

Although a variety of technologies have been proposed for improving surface 
characteristics, surface modification of polymers by grafting is a rather new 
technology. It offers versatile means for incorporating new functionalities into 
existing polymers. However, in spite of the potentially wide applications of such 
surface grafting technology, this has been applied only to a few cases in indus- 
try, probably because the basic studies required for the applications are still in 
their infancies. Another reason may have been that such a grafted polymer sur- 
face is relatively expensive to produce and is used mostly in aqueous environ- 
ments which are rarely encountered in the conventional industrial applications 
of polymers. Examples of application in such aqueous environments include 
marine science, biotechnology, and biomedical engineering. As these fields have 
become active only in recent years, it is not surprising that polymer modifica- 
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tion by surface grafting has not yet become popular in industry. Recently, Hoff- 
man summarized surface modification technologies by means of physical, 
chemical, mechanical, and biological methods [ 1 ] . 

Characterization of the structure of grafted surface is crucial, because abetter 
understanding of the grafted surface will be a key factor for further developing 
this new technology. In this article, various methods of grafting and surface char- 
acterizations, and applications associated with grafting from recent literature 
publications will be reviewed. 

2 

Creation of Grafted Surfaces 

There are in principle two methods for producing grafted surfaces, as schemat- 
ically illustrated in Fig. 1; direct coupling reaction of existing polymer molecules 
to the surface and graft polymerization of monomers to the surface. Each has its 
own advantages and disadvantages, as demonstrated below. 

2.1 

Polymer Coupling Reactions 

If the polymer surface to be modified possesses reactive groups capable of com- 
bining other components, such as water-soluble polymer molecules, surface mod- 
ification can be readily conducted by chemical coupling reaction. Numerous syn- 
thetic reactions are available for this purpose. Bergbreiter [2] reviewed various 




Fig. la, b. Creation of grafted surfaces by: a direct polymer coupling reaction; b graft polymer- 
ization 
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technologies and analytical methods associated with chemical modification of 
polymer surfaces. Kramer [3] proposed expressions for kinetics of grafting of a 
pure end-functional polymer melt to a reactive interface assuming that the free 
energy of reaction was very large and negative. Kishida et al. directly immobi- 
lized poly(ethylene glycol) (PEG) chains onto a cellulose surface through esteri- 
fication [4] . Prior to coupling reaction, they derivatized the terminal hydroxyl 
group of PEG molecules to carboxylic acid using succinic anhydride. This PEG 
molecule with the terminal carboxylic acid was then chemically immobilized to 
the hydroxyl group on the cellulose surface using carbodiimide in non-aqueous 
media like toluene. Tezuka et al. [5] immobilized block copolymers onto poly- 
(vinyl alcohol) (PVA) and polyurethane (PU) surfaces. Polystyrene-polydimethyl- 
siloxane (PS-PDMS) block copolymers containing vinylsilane or diol at the end 
of the chain (“macromer”) were synthesized by living polymerization of styrene 
and dimethylsiloxane. Then the “macromer” was allowed to react chemically 
with the PVA or PU surface. Han et al. prepared PEO-grafted PU beads [6]. PU 
beads were first treated with hexamethylene diisocyanate in toluene with stan- 
nous octoate and subsequently chemically grafted with PEO in benzene in the 
presence of the same catalyst. Because a triethoxy silane group strongly bonds 
to a silanol group in the surface of silicon oxide, Clark et al. [7] converted a car- 
boxy-terminated polystyrene to a triethoxy silane end polymer and carried out 
its surface grafting onto a cleaned silicone wafer. Ester exchange reaction was 
performed for grafting of 2-ethylhexyl acrylate and methyl methacrylate poly- 
mer [8] and the resulting copolymer was exchanged with PEG monomethyl ester. 

2.2 

Graft Polymerization 

A variety of methods for graft polymerization onto different substrate surfaces 
will be described below from relatively recent reports. The methods include chem- 
ical graft polymerization and grafting with the use of high-energy radiation or 
oxidizing agents. For oxidizing a substrate surface, ozone, acidic treatments, and 
high-energy radiations in air can be used. Ionizing and nonionizing radiations 
as well as plasma treatment are usually selected for the high-energy sources. Stro- 
bel et al. [9] compared five gas-phase surface oxidation processes, i. e., corona dis- 
charge, flame, remote air plasma (Fig. 2), ozone, and combined UV/ozone treat- 
ment onto PP and poly(ethylene terephthalate) (PET) films. They revealed that 
flame, corona, and remote-plasma processes could readily oxidize the polymer 
surfaces, the 0/C ratio being greater than 0.1 when treated for 0.04, 0.05-0.5, and 
0.2 s, respectively. Of the five techniques, the frame treatment appeared to give 
the shallowest oxidation depth near the outermost surface region. In case that 
oxidized species is employed for initiation of graft polymerization, it should be 
noted that any method is not always applicable to all polymeric materials. For 
instance, fluorinated polymers such as polytetrafluoroethylene (PTFE) do not 
undergo significant oxidation by ozone and corona treatment, but active species 
for graft polymerization are effectively produced by glow discharge treatment of 
these polymers. Fluorinated polymers often suffer serious chain scission upon 
exposure to vacuum ultraviolet (VUV) radiation, but PET and polyethylene (PE) 
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Air flow 




Fig. 2. Schematic diagram of the remote-plasma treater (Reproduced with permission from 
Strobel et al,, J Adhesion Sci Technol 9; 365 Copyright (1995) VSP) 



do not [10]. Ichijima et al. [11] compared the concentrations of peroxides intro- 
duced onto a poly( methyl methacrylate) (PMMA) substrate by different initia- 
tion methods, including ozone exposure, UV irradiation, corona discharge, and 
glow discharge. They evaluated their effect on the yield of graft polymerization 
of acrylamide (AAm) onto the oxidized PMMA surface. The most suitable initi- 
ation method for graft polymerization of AAm onto the PMMA substrate was 
found to be UV irradiation. 

2.2.1 

Direct Chemical Modification 

Bergbreiter and Bandella [12] prepared a pH sensitive surface by immobilizing 
poly(acrylic acid) (PAAc) chains onto a PE surface. According to their technolo- 
gy, PE film surface was first functionalized by phenylpyrenylm ethyl group, fol- 
lowed by graft polymerization of t-butyl acrylate monomer and subsequent 
hydrolysis to form PAAc. Pyrene groups were found to localize near the interfa- 
cial region, and the average degree of polymerization (DP) of acrylate polymer 
was estimated to be 20-30. They also reported an oxidation method by chromic 
sulfonic acid [13]. The esterified PE surface was then UV-irradiated in order to 
introduce radicals on the substrate, followed by graft polymerization of acry- 
lonitrile. 

The surface of polymeric materials such as polypropylene (PP), PS, polyacry- 
lonitrile (PAN), and nylon was oxidized by immersing them in aqueous solution 
of oxidizing agents such as potassium peroxy disulfate under nitrogen purging 
at high temperatures [ 14] . Graft polymerization of water-soluble monomers such 
as AAm, methacrylic acid, and 3-aminopropyl methacrylate has been frequent- 
ly performed in aqueous solution with the use of ceric ion, for instance, at 
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2x10“^ mol/1 ceric ammonium nitrate in nitric acid at 50°C. A dewaxed raw jute 
fiber was surface-grafted with methyl methacrylate (MMA) polymer using an 
initiator from the combination of IO4 and Cu^^ ions [15]. The initiation of graft- 
ing was attributed to radicals generated from cellulose through complexing 
between aqueous periodate ions and the anhydroglucose unit of cellulose mole- 
cule, followed by oxidation of the latter at the 1,2-glycol position. Grafting effi- 
ciency of 20-25% was obtained when an optimal polymerization condition 
([Cu^+] =0.001 mol/1 and [IO4 ] =0.005 mol/1) was used. A redox system consist- 
ing of K 2 S 20 g-Na 2 S 203 -Cu^+ was used for surface graft polymerization of 
methacrylic acid (MAA) onto nylon 6 fibers and grafting took place mostly on 
the non-crystalline region of the substrate [16]. PET fibers were graft polymer- 
ized with a monomer mixture of A Am and MAA using AIB N [ 1 7] . The gr aft yield 
of the AAm polymer increased with the increasing fraction of MAA in the 
monomer mixture. 

2 . 2.2 

Ozone 

PET, PE, PS, PP, and PC films were pretreated with ozone to introduce peroxy 
groups onto the substrate surface and AAm,AAc, and sodium styrenesulphonate 
(NaSS) monomers were graft polymerized by a near-UV induced technique [18]. 
The authors also performed surface graft polymerization of these monomers 
onto a polyaniline film. Although the pristine polyaniline film was effectively sur- 
face grafted by the near-UV induced method, combination of this method with 
ozone pretreatment enhanced the graft yield [19]. Buchenska [20] modified 
polyamide (PA6) knitted fabrics by graft polymerization of AAm and the hydroxy 
peroxide groups were introduced onto the fabrics surface by direct oxidation with 
oxygen in the temperature range from 120 to 160 °C prior to graft polymeriza- 
tion. The graft density obtained was 0.14xl0'^mol/g and subsequent treatment 
of the grafted fabrics with hydrazine or 3-bromopropionic acid improved both 
hydrophilic and antistatic properties of the fabrics. 

2.2.3 
y -Rays 

y-Ray induced graft polymerization can be carried out by the following meth- 
ods: 1. simultaneous irradiation and grafting through in situ formed free radi- 
cals, 2. grafting through peroxide groups introduced by pre-irradiation, and 3. 
grafting initiated by trapped radicals formed by pre-irradiation. Stannett [21] 
had reviewed the past 35 years’ development of grafting via ionizing radiation. 
He recommended the use of electron beam radiation for polymer surface mod- 
ification. Done et al. [22] modified an industrial PP film by graft polymerization 
of AAc. The y-ray irradiated PP film was immersed in a degassed aqueous solu- 
tion of 50 wt% monomer and 0.1 wt% Mohr’s salt at 85 °C for 2 h. Mohr’s salt was 
added into the solution to prevent homopolymerization of the monomer and the 
peroxides yielded upon irradiation were quantitatively measured with the DPPH 
method. Poly(vinyl chloride) (PVC) powders were graft polymerized with AAc 
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after the introduction of hydroperoxides onto the substrate surface by y-ray irra- 
diation in air at a dose rate of 84 MR/h from a 2100 Ci ®°Co source [23]. Percent 
graft yields were 100-200%, indicating that grafting was not localized in the sur- 
face region. The presence of ethanol, a chain transfer reagent, on the other hand, 
remarkably reduced the graft yield. Polymers of A Am and 2 -hydroxy ethyl 
methacrylate (HEMA) were grafted onto an ethylene-propylene copolymer rub- 
ber (EPR) by the radiation technique to improve the water uptake, wettability, 
and biocompatibility [24] . Eang et al. [25] conducted graft polymerization of vinyl 
acetate onto EPR by y-irradiation to percentage graftings as high as 40%. Jan et 
al. [26] modified a porous polyfvinylidene fluoride) membrane to introduce 
positive charge by y-irradiation in the presence of vinyl-triphenyl-phosphoni- 
um bromide monomer. In addition, PTFE-FEP copolymer [27] , hollow fiber [28] , 
PDM rubber [29] , and PP fiber [30] were modified by graft polymerization with 
the radiation method. 

2.2.4 

Electron Beams 

Polymers of MMA, AAc, and MAA were grafted onto an ultrahigh molecular 
weight polyethylene (UHMWPE) fiber surface after pretreatment with electron 
beam irradiation [31]. Sundell et al. [32] pretreated a PE film with electron 
beams to facilitate the graft polymerization of vinyl benzylchloride onto the sub- 
strate. The inner surface of porous PE hollow fiber had also been modified by 
grafting of glycidyl methacrylate (GMA) polymer after electron beam irradia- 
tion [33]. 

2.2.5 

Glow Discharge 

Low-temperature plasmas have been extensively applied to modify surface prop- 
erties of polymer. Recently, surface modification by means of macromolecular 
plasma chemistry was reviewed by Denes [34], who also described the precise 
principle of plasma chemistry in detail. He pointed out that it is extremely diffi- 
cult to understand and control the reaction mechanism of plasma because the 
relatively wide energy range (0.5-5 eV) of the glow discharge generates a very 
large number of charged and neutral molecular fragments. Various energy 
sources, such as UV, VUV, electron beams, and ozone are also involved in the plas- 
ma treatment if oxygen or air is present. Major applications of low-temperature 
plasma include plasma polymerization, which deposits a crosslinked thin poly- 
meric layer on the substrate surface, and plasma treatment, which causes inten- 
sive oxidation in the surface region of the substrate. 

It should be stressed that polymer surface can also be modified by graft poly- 
merization utilizing free radicals or peroxides generated by the plasma treatment, 
similar to the effect of irradiation with high-energy radiations. The main differ- 
ence between plasmas and high-energy radiations is the density of polymer rad- 
icals generated by exposure of a polymer substrate to these high-power sources. 
Sheu et al. [35] immobilized poly( ethylene oxide) (PEO) surfactants on hydropho- 
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Fig. 3. Schematic diagram of argon glow discharge treatment for the immobilization of sur- 
factants (Reproduced with permission from Sheu et al.. Contact angle, wettability and adhe- 
sion, p 531 Copyright (1993) VSP) 

bic polymer surfaces using the crosslinking by activated species of inert gases 
(CASING) technique. Low density polyethylene (LDPE) was first coated with alkyl 
PEO surfactants and then the substrate surface was crosslinked by Ar plasma treat- 
ment (see Fig. 3). Non-fouling (protein resistant) polymer surfaces were report- 
ed to be obtained by this method. A PEO-PPO-PEO tribrock copolymer surfac- 
tant was also utilized for the CASING technique [35]. Hsie and Wu [36] exposed 
a PET film to the glow discharge of Ar plasma for 0.5-5 min at a gas pressure of 
0.5 torr and a radio -frequency of 13.56 Hz. They then performed surface graft 
polymerization of AAc onto the plasma-treated film both in the vapor and liq- 
uid phases without exposing to air. Freshly synthesized PE powders were also 
treated with glow discharge, followed by surface graft polymerization of styrene 
without exposing to air [37] . ESR study revealed that graft polymerization start- 
ed from the alkyl radicals introduced by the glow discharge and that more fresh- 
ly synthesized powders were more amenable to form active radicals. 

They also studied the melting characteristics of the powders with surface- 
grafted MMA polymer by DSC [38]. Grafting of Kevlar 49, or poly(p-phenylene 
terephthalamide), fiber with GMA and AAm polymers without exposing to air 
was studied by Yamada et al. [39] . According to their grafting system, Kevlar yarn 
was suspended in an upper plasma-irradiation vessel and irradiated for 1 min at 
an Ar gas pressure of 20 Pa, an rf frequency of 13.56 MHz, and an rf power of 
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100 W. The yarn was then immersed immediately in a GMA monomer solution 
placed at the lower side of the vessel (Fig. 4). They measured the adhesive strength 
between a resin matrix and the grafted Kevlar fiber. A similar method was 
employed by Yamaguchi et al. [40] for the graft polymerization of porous PE with 
methylacrylate monomer without contacting with air. It was found that radicals 
were formed in the porous substrate. Textile fabrics including cotton, cellulose 
acetate, and cupra ammonium cellulose were subjected to graft polymerization 
of HEMA, AAm,iY-isopropylacrylamide (NIPAM), AAc, and 2-methoxyethyl acry- 
late (MEA) monomers in the absence of air after exposing them to 13.56 MHz rf 
Ar plasma [4 1 ] . In most cases, the graft yield obtained was 40-80%. Lai et al. [42] 
improved the adhesive strength of a silicone rubber against a conventional adhe- 
sive tape. The silicone rubber was pretreated with O 2 plasma and exposed to air 
in order to introduce peroxide groups onto the substrate surface. The pretreated 
rubber was immersed in aqueous monomer solution of AAm or AAc. There exist- 
ed an optimum period of time for glow discharge treatment and graft polymer- 
ization with respect to both peroxide formation and adhesive force. Castner et 
al. [43] modified a silicone rubber by graft polymerization of AAc and charac- 
terized the surface changes using XPS, SIMS, and ATR-FTIR. Peroxides intro- 
duced onto the rubber surface were determined with the DPPH method. Lee and 
Shim [44] prepared a pH-sensitivepoly(vinylidene fluoride) membrane by graft- 
ing with AAc polymer. The film was treated with Ar plasma for up to 3 min and 
then exposed to air, followed by graft polymerization of AAc in deaerated 20 wt% 
aqueous solution at 60 °C. The degree of polymerization (DP) for the homopoly- 
mer ranged from 32 to 94, depending on the plasma exposure time. A Kapton 
(poly(Ar,N’-oxydiphenylene pyromellitimide)) film pretreated with Ar plasma 




Fig. 4. Schematic representation of apparatus for plasma-graft polymerization (Reproduced with 
permission from Yamada et al., J Appl Polym Sci 60; 1847 Copyright ( 1 996) John WUey 8t Sons, Inc.) 
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was grafted with vinylimidazole polymer after exposing the plasma-treated film 
to air [45 ] . Grafting was also performed in benzene solution at 40-80 °C for 24 h 
to have a graft yield around 0.8 wt% in order to improve adhesion between the 
Kapton film and copper [46]. Qiu et al. [47] treated PU film with plasma at an rf 
of 2.45 GHz and a power of 300 W, exposed to air for 10 min, and then subjected 
the film to graft polymerization of PEG methacrylate (PEGMA) monomer. Hsi- 
ue et al. [48] pretreated a silicone rubber and a TPX film with Ar plasma and then 
subjected them to surface-grafting of polyHEMA to enhance cell adhesion onto 
the surfaces. These authors used also AAc [49], HEM A, and PEO for surface graft 
polymerization [50, 51]. Argon plasma treated PET, PE, and PP films were sur- 
face-grafted with N,N-dimethylacrylamide (DMAA) polymer to render the sur- 
face lubricious [52] . When the films were pretreated with benzoyl peroxide pri- 
or to graft polymerization, the graft amount was further increased. Plasma- 
induced graft polymerization was also performed on membranes of PAN and 
polysulphone using AAc as the monomer [53]. 

2.2.6 

Corona Discharge 

Both glow and corona discharges are currently used for surface modification of 
polymeric materials. Corona discharge treatment is by far simpler than glow dis- 
charge treatment, because the former technique does not require evacuation of 
the discharge system. However, corona treatment can result in more damage to 
the substrate polymers and the discharge conditions are more difficult to con- 
trol than those of plasma treatment. The latter treatment is generally carried out 
under reduced gas pressure. 

The surface of a low density PE sheet was modified by graft polymerization 
of ethylene glycol methacrylate with the use of corona discharge pretreatment 
[54] .A knife-type electrode, 1.5 mm away from the PE sheet surface, was connect- 
ed to an rf generator of 50 W and 100 kHz. After corona treatment, grafting was 
performed at 80 °C for 4 h under continuous nitrogen bubbling to yield a comb- 
like PE surface. Various types of PP materials were subjected to surface treatment 
with corona discharge and the hydrophilicity and hydroperoxide yield were eval- 
uated [55]. Isotactic PP was found to be more easily oxidized due to high crys- 
tallinity. The same authors also performed surface graft polymerization of AAm 
onto the pretreated PP surface. Lee et al. [56] grafted a corona-treated PE sheet 
with AAc polymer and subsequently subjected the surface to chemical reaction 
with various ionic reagents. Corona-treated PE films were grafted with various 
water-soluble polymers and the amounts of cell adhesion onto the grafted sur- 
faces were compared [57]. 

2.2.7 

UV Irradiation 

The use of UV irradiation appears to be an excellent method for surface graft- 
ing of polymers because of its simplicity and cleanliness. Additional reasons for 
the suitability of the photochemical method for surface grafting of polymers are 
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as follows: 1) photochemically produced triplet states of carbonyl compounds 
can abstract hydrogen atoms from almost all polymers so that graft polymeriza- 
tion may be initiated; 2) high concentrations of active species can be produced 
locally at the interface between the substrate polymer and the monomer solu- 
tion containing a sensitizer when UV-irradiation is applied through the substrate 
polymer film; 3) in addition to the simplicity of the procedure, the cost of ener- 
gy source is lower for UV radiation than for ionizing radiation. 

Ranby et al. [58] suggested possible applications of surface modification to syn- 
thetic polymer fabrics by photo-induced graft polymerization. They performed 
graft polymerization of AAc and AAm onto the surface of fibers and films of PP 
and HDPE utilizing photo-initiators, such as benzophenone, 4-chlorobenzophe- 
none, and hydroxycyclohexylacetophenone. Ranby [59] also discussed the adhe- 
sive properties of surface-modified materials. In order to functionalize the cellu- 
lose surface, graft polymerization of GMA was carried out in a Pyrex tube con- 
taining a cellulose sample, GMA monomer, and H 2 O 2 [60, 61]. The H 2 O 2 was 
decomposed by photo irradiation to produce hydroxyl radicals which might 
extract the hydrogen atom from polymer substrates to yield polymer radicals 
capable of initiating grafting. Graft yield depended on the species of photosensi- 
tizers used. They also performed graft polymerization using an initiator con- 
taining a cerium salt. A similar method was employed for graft polymerization of 
AAm on an ethylene-vinyl alcohol copolymer film [62] and for graft polymeriza- 
tion of AAc and 4-vinyl pyridine on a PE film [63] . In the latter case, the film was 
previously coated with benzophenone solution. Edge et al. [64] graft polymerized 
different monomers in the vapor phase onto a polyetherimide surface in the pres- 
ence of a benzophenone solution in acetone. The identical method was also applied 
to the graft polymerization of some monomers onto low density PE and PET sur- 
faces both in air and in N 2 , and the surface properties were compared after char- 
acterizing with contact angle and XPS measurements [65]. PE, PP, PET, and PVA 
fabrics were surface-modified by photo-induced graft polymerization of AAc, 
AAm, glycidyl acrylate, and vinyl pyridine to improve the dyeing property of fab- 
rics [66]. A block copolymer membrane consisting of PE, PS, and polybutadiene 
was subjected to graft polymerization of AAm in the presence of a small amount 
of isopropanol [67]. Water-soluble benzophenone derivatives were utilized to ini- 
tiate UV-induced graft polymerization onto the membrane and the reaction 
mechanisms were discussed. Ulbrecht et al. [68] modified an ultrafiltration mem- 
brane with grafted PEG. Lee et al. [69] performed UV-induced graft polymeriza- 
tion of AAc, M AA, and NIPAM onto a polyamide membrane with combination of 
plasma treatment. Permeability of riboflavin depended on temperature when the 
NIPAM graft polymerized membrane was subjected to flow testing. A PE plate 
surface was surface-grafted with MAA polymer in the liquid phase after coating 
with a benzophenone (sensitizer) and poly( vinyl acetate) (PVAc) mixture [70,71]. 
After drying, the coated plate was immersed in an aqueous solution of MAA at 
60 °C under a nitrogen atmosphere. Nakayama and Matsuda [72] performed sur- 
face photo -grafting using photo -reactive groups. Two different surfaces were graft- 
ed: a PS film whose surface was partially derivatized with N,N-diethyl-dithiocar- 
bamyl groups prior to UV irradiation, and a PET film which was precoated with 
a photo -reactive copolymer cast from the toluene solution prior to UV exposure. 
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Photo -induced graft polymerization was also performed without using any 
photosensitizers and degassing procedure. Kang et al. [73] graft polymerized 
water-soluble monomers such as AAc, NaSS, and DMAA onto a PIPE film sur- 
face by near-UV irradiation with the combination of Ar plasma pretreatment. 
The plasma-pretreated film was exposed to air and then placed in an ampoule 
containing deaerated aqueous monomer solution, followed by exposure to a high- 
pressure mercury lamp. Grafting was also performed without the degassing pro- 
cedure using riboflavin which consumed the dissolved oxygen in the course of 
photo irradiation. The surface of electroactive polymers such as polyaniline [74] , 
polyimide [75], and polypyrrole film [76] was modified by similar graft poly- 
merization of various water-soluble and ionic monomers. To improve the adhe- 
sion between epoxy molding compound (EMC) and ball grid array (EGA) sub- 
strate in the packaging of microelectronics systems, the surface of the epoxy- 
based EGA was graft polymerized with GMA [77] ,as depicted in Fig. 5. The authors 
also modified the surface of PE, PET, and PS films, and characterized the surface 
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Fig. 5. Schematic diagram for the process of surface graft polymerization of GMA onto GA and 
subsequent reaction with amine curing agent in EMC 
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chemical structure [78]. Quaternary amines were introduced onto a PET film 
through photo-induced graft polymerization [79] . The UV-induced surface graft 
polymerization method was applied to UHMWPE and PET fibers for improved 
adhesion with an epoxy resin [80]. 

A laser beam was used for graft polymerization of AAc onto a tetrafluo- 
roethylene-perfluoroalkyl vinyl ether copolymer film [81]. The film placed in 
contact with AAc solution was irradiated with KrF laser through the film to excite 
the film/solution interface. Surface composition of the grafted film determined 
by XPS revealed an extensive loss of fluorine atom and an increase of oxygen 
atom in addition to the presence of a Cls line shape, similar to that of AAc 
monomer. Mirzadeh et al. [82] used pulsed laser beam for the graft polymer- 
ization of AAm on a rubber surface in the presence of a photosensitizer, ben- 
zophenone, or AIBN. 

Vacuum ultraviolet (VUV) irradiation was utilized for the graft polymeriza- 
tion of a mesogen-containing monomer onto fluorine-containing polymeric 
materials during the preparation of an anisotropic liquid crystal layer [83]. The 
grafted substrate comprised of two layers - a fluorine-containing polymer which 
provided good physical-mechanical properties, and a grafted chain layer which 
exhibited the usual features of nematic liquid crystal. 

3 

Characterization of Grafted Surfaces 

Besides hydrophilicity and wettability, the atomic composition, overall molecu- 
lar size, shape, mobility and conformational relaxation time of graft chains are 
of importance to the exploration physical and chemical structures of grafted sur- 
faces. For these characterizations, a variety of analytical methods including spec- 
troscopy, scattering, ellipsometry, and hydrodynamics have been proposed, in 
addition to computer simulation. Small angle neutron scattering is one of the 
most powerful techniques in determining the overall size of polymer chains and 
the spatial distribution of the chain segments [84]. The intensity and phase of 
scattered neutrons depend on the position and structure of the nucleus. On the 
other hand, neutron reflection, which is a relatively new technique, provides a 
unique characterization method for the atomic composition as a function of the 
depth of thin films. Fourier transform infrared spectroscopy coupled with atten- 
uated total reflection (ATR-FTIR), Auger spectroscopy, and X-ray photoelectron 
spectroscopy (XPS) have become classical methods, but still provide very valu- 
able information regarding the constituent elements and chemical structure 
near the surface region. The chemical derivatization techniques are also useful 
for surface analysis. 

3.1 

Structure of Graft Polymer Chains 

Graft chains or long-chain polymer molecules attached by one end to a surface 
or the interface in contact with solvents are often called “polymeric brushes”. 
A number of studies on the structure of such polymer chains have been 
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undertaken using various theoretical means and computer simulations. The 
techniques employed include self-consistent field (SCF) theory [85-93], mole- 
cular dynamics (MD) [94, 95], and Monte Carlo (MC) simulation [96-100]. 
Detailed theoretical calculations revealed that long flexible polymer chains 
exhibit stretching under high graft densities, resulting in much greater height 
of graft chains than the radius of gyration. Irvine et al. [92] developed an SCF 
treatment of star-shaped polymer graft chains tethered to an impenetrable 
surface and used the model to calculate the near-surface equilibrium segment 
density profiles. At a moderate graft density, graft star molecules exhibited 
maximal stretching. Although these analyses promote an understanding of 
molecular structures under a variety of conditions or various environments, 
e. g., under 0 point [100] and with or without solvent [95], it is not easy to prove 
experimentally the predicted molecular structure because of the difficulty in 
preparing the model graft chains. In addition, grafted surfaces are very diffi- 
cult to characterize in detail because of the extremely low densities of graft 
chains on the surface. For instance, graft density amounts to only 0.1 pg/cm^, 
if graft chains with a molecular weight of 1x10^ are fixed as a monolayer at a 
frequency of one chain per 100 nm^ on a surface. Conventional analytical 
methods are not sensitive enough to allow us to determine separately the 
length and the number density of graft chains fixed on a surface at such a low 
concentration, although several methods are available to evaluate the overall 
graft density which is a product of the chain molecular weight and the num- 
ber density of the graft chains. 

Concentration profiles of PS graft chains were studied by a neutron reflection 
method [101]. Graft chains consisted of end-functionalized deuterated poly- 
styrene. Although the measurement was not performed under the true equilib- 
rium conformation, the observed metastable state was in good agreement with 
that predicted from the SCF theory. The kinetics of the penetration of graft 
chains into the polymer matrix was also investigated. 

Using angle-resolved XPS, the chemical composition and structure of film sur- 
faces grafted with water-soluble polymers were investigated [78] . The XPS results 
showed that, in the case of substantially high grafting, the graft polymer chains 
penetrated or became partially submerged beneath a thin surface layer much 
richer in the substrate polymer. The surface structure was further confirmed 
using dynamic water contact angle measurements. 

3.2 

Thickness of Graft Chains 

Some structural studies have been performed using simplified graft model chains. 
A long polymer chain occasionally possesses a number of functional sites which 
act as attachment points for the polymer chain to bind to a substrate surface as 
illustrated in Fig. 6. When a silica surface is placed in a melt or a solution of PDMS 
for a sufficiently long period of time,some segments become permanently bound 
to the solid surface. This type of graft layer is called pseudobrush or Guiselin brush. 
Cohen Addad et al. [102] determined the dry thickness of the pseudobrush pre- 
pared from various molecular weights of PDMS and concluded that the dry thick- 
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Fig. 6. Idealized representation of a “pseudobrush” or Guiselin brush 



ness h had a linear correlation between and the monomer size a. h is propor- 
tional to the size of a mobile chain and N is a number of segments in one poly- 
mer chain. They determined the dry thickness by washing out all the unbound 
polymers by a good solvent followed by drying and weighing the residual poly- 
mer. Raudino and Zuccarello [103] derived a simple equation expressing the mean 
height of graft chains, having a statistical distribution of the internal conforma- 
tions and hydrogen bonding, by combining statistical and quantum mechanical 
methods. 

Wyart et al. [104] prepared a series of pseudobrushes by simultaneous adsorp- 
tion of hydroxy-terminated PDMS onto a silicon wafer surface and reported that 
the dried thickness of these layers obeyed the scaling law of where 

(|)o denotes the initial volume fraction of occupied polymer chains resulting from 
the interdigitation through a rubber matrix [105]. This relation was confirmed 
by other independent measurements such as neutron scattering [106], ellipsom- 
etry, and X-ray reflectance [107]. They also showed the existence of a maximum 
value for the 90° peel off adhesive energy at <))o^^^=0.4 ~ 0.5. Cosgrove et al. [108] 
determined the hydrodynamic thickness by photon correlation spectroscopy. 
Silica particles having radii ranging from 0.08 to 0.2 pm were grafted with PS of 
Mw of 12 000-24000. The graft density of PS was from 30 to 65 pg/cm^. The 
hydrodynamic radius of particles was calculated from the Stokes -Einstein equa- 
tion using the measured diffusion coefficient and the hydrodynamic thickness 
was obtained from the difference of the radii between the surface-grafted and 
the bare silica particles. The root mean square of the layer thickness could be fit- 
ted to the Alexander-de Gennes brush model, indicating that the chains were 
highly extended. 

If a neutral layer of graft chains is present near a charged surface, the elec- 
troosmotic motion of the electrical double layer through the interface will be 
retarded. The reduction in electrophoretic mobility due to the neutral polymer 
chains assumes the hydrodynamic layer thickness. Janzen et al. [109] measured 
the electrophoretic mobility as a function of ionic strength for a PEG -grafted lipid 
bilayer. The mobility conformed to a model in which the PEG chains formed a 
surface layer of polymer with viscous drag arising from electro osmotic flow and 
provided the average layer thickness. 

The hydrodynamic thickness of graft chains both in nonsolvated and solvat- 
ed states was experimentally quantified by Webber et al. [110]. When a poly(2- 
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vinylpyridine)/polystyrene (PVP/PS) diblock copolymer was adsorbed in the 
pores of well-characterized mica membranes from toluene, the PVP block 
anchored on pore wall while the PS block imitated a terminally attached graft 
chain. They calculated the thickness of polymer chains from pressure drop as a 
function of flow rate through the membrane utilizing the modified Hagen- 
Poiseuille equation. Reversible extension of the polymer chains was observed 
when the solvent was altered from toluene to heptane or from heptane to 
toluene during the more than one month run. They further demonstrated that 
there was no shear thinning of the polymer layer over the shear rate range of 
lO^-lO^s-i. 

The flow rate retardation due to the attached graft chains was also observed 
using a porous membrane [111]. In this work, a porous poly(vinylidene fluo- 
ride) membrane was surface-grafted with PNIPAM which is soluble in water but 
has a lower critical solution temperature (LOST) around 31-33 °C. The flux of 
pure water through the grafted membrane varied by more than ten times 
between the temperatures above and below the LCST. The temperature sensi- 
tivity was reversible and reproducible. In this experiment, the thickness of graft 
chains could not be determined because no precise data for the pore size and 
density were available. The approximate layer thickness of water-swollen graft 
chains was directly measured using PNIPAM-grafted PU surface [112]. A virgin 
PU film became thicker under a temperature rise from 26 to 50 °C due to the 
thermal expansion, whereas the PNIPAM-grafted PU film became thinner 
under the same temperature rise because of the phase transition of the graft 
polymer chains. The thickness was of the order of a micrometer when estimat- 
ed from the difference of the thickness change. The thickness of grafted DMAA 
and AAm polymer chains was estimated from the hydrodynamic flow using a 
PU tube (ID=0.65 mm, length=200 cm) [113]. The inner wall of PU tube was 
graft-polymerized with water-soluble monomers by an ozone pretreatment 
method. The flow rate of water through the PU tube in the range of laminar flow 
(Reynolds’ number less than 10) was measured at a constant temperature. Under 
the simple assumption that the flow retardation was caused by the narrowing of 
the inner diameter of the tube in the presence of the surface graft layer, the graft 
thickness was estimated to be in the order of several micrometers. The flux 
change through a porous filter was also intensively discussed by de Gennes 
[114]. The increase in flux through a thin film composite membrane was report- 
ed by Mukherjee et al. [115] using a surface-modified membrane by chemical 
fluorination. 

With respect to the dynamic structure of graft chains, direct and relevant exper- 
iment is difficult to perform, although some approaches have been proposed from 
theoretical analysis. Fytas et al. [116] reported on the dynamic behavior of brush- 
like chains using evanescent-wave dynamic light-scattering (EWDLS) which has 
an advantage in obtaining a time dependence of the deformation within a time 
scale of 10“^ to 10^ s. They used an asymmetric PEO-PS block copolymer com- 
posed of two distinct blocks; a long PS block and a shorter PEO fragment. Dis- 
solving the block copolymer in toluene, the PEO end-group was anchored onto 
a quartz prism surface, while PS chain imitated a stretched graft chain having 
the average layer thickness ranging from 45 to 130 nm. The evanescent wave prop- 
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Fig. 7. Schematic diagram of the evanescent wave dynamic light-scattering instrumentation 
(Reproduced with permission from Fytas et al., Science 274; 2041 Copyright (1996) American 
Association for the Advancement of Science) 



agated through the graft chains under the condition of total internal reflection 
(see Fig. 7). 

Few studies were performed on the thickness of ionic graft chains. vonGoel- 
er and Muthukumar [117] evaluated the height of a polyelectrolyte in solvents of 
varying salt concentrations in the high grafting regime. Integrating the electro- 
static potentials over the directions parallel to the surface but neglecting the den- 
sity fluctuation, they found that the dimensionless parameters and 
were valid for describing the long-ranged electrostatic interaction, where v, k, 
and I were the Coulombic interaction of strength, the inverse Debye screening 
length, and the Kuhn length, respectively. For instance, the chain height h was I 
]yi/2for large values of both and and h was with Kx<3 for large 
yN^l 2 was small. They also discussed the phase diagram of the 

scaling behavior and the transition of the polyelectrolyte in a poor solvent [118]. 
Experimental studies on the interaction between polyelectrolyte graft chains are 
also rare. 
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4 

Applications 

4.1 

Non-Medical Application of Grafted Surfaces 

4.1.1 
Adhesion 

When two chemically identical rubbers are brought into close contact with each 
other, a population of mobile chains at the rubber-rubber interface promote the 
adhesion between them, de Gennes [119] discussed the penetration of mobile graft 
chains in a rubber having an identical chemical structure. Such graft polymer 
chains are called “adhesion promoters” or “connectors”. Raphael and de Gennes 
[120] calculated the adhesion energy GpuHout of the connectors grafted on one block 
and free at the other end, and thus being pulled out without any chemical rupture. 
They also compared the adhesion energy with that associated with chemical bond 
Gscission ^od arrived at the following conclusion: Gscission/ Gpull out is nearly equal to 
Ux/Uv, where Ux is a typical van der Waals energy between two adjacent monomers 
and Uy represents the chemical bond energy. Some experimental results did not 
satisfy this relation, probably because the chains employed in the experiment were 
too short or because of the formation of bundles due to poor solvent. 

Ultra-high modulus fibers such as aramid and carbon fibers have been cur- 
rently utilized for composite material fabrication. Ultra-high modulus polyeth- 
ylene (UHMPE) fiber is also applicable for composite fabrication because of the 
light weight in addition to its high modulus, vibration damping, and resistance 
to chemicals. However, this fiber has drawbacks such as poor interfacial adhe- 
sion with the polymer matrix of the composite because of highly hydrophobic 
nature of the fiber surface. 

Li and Netravali [121] modified an UHSPE fiber surface using allylamine plas- 
ma deposition to improve its adhesion to epoxy resin. The shear strength of the 
modified surface increased by a factor of 2-3 when measured by pull-out test 
using a U-shaped epoxy resin. The epoxy resin and a curing agent were mixed 
and molded into U-shaped form without disturbing the fiber (see Fig. 8). 




A 



B 



Fig.S.a Silicone rubber mold, b Epoxy resin specimen with embedded fiber ready for single- 
fiber pull-out test (Reproduced with permission from Li and Netravali, J Appl Polym Sci 44: 333 
Copyright (1992) John WUey & Sons, Inc.) 
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The surface of Kevlar 49 fiber was chemically surface modified with epoxy 
groups containing reagents such as epoxy tetraglycidyl ether of diphenyldiamino 
methane [ 122] . It was found that the modified Kevlar fiber had enhanced adhesion 
to the epoxy matrix resin. Penn and Jutis [123] attached amine-terminated pen- 
dant groups to the surface of aramid fiber by chemical reaction. A single filament 
pull-out test used to assess the effect of the pendant groups on the fiber-epoxy 
matrix bond strength revealed that the presence of the pendant groups increased 
the adhesive performance. Ranby et al. have performed graft polymerization of A Ac 
[ 124] , glycidyl acrylate, and GMA [125] onto the surface of low density polyethyl- 
ene with the use of UV irradiation in the vapor phase utilizing benzophenone as 
a photosensitizer. The surface of Kevlar fiber was also modified by Mori et al. [ 126] 
with the use of UV irradiation technique. To improve the Kevlar surface adhesion 
to epoxy resin, graft polymerization of AAm and GMA was performed onto the 
fiber surface. Following plasma treatment and the subsequent exposure to air, the 
fiber was placed in the monomer solution and subjected to UV irradiation. ATR- 
FTIR and XPS measurements revealed that the graft chains were present in the sur- 
face region of Kevlar fiber with the presence of epoxy groups on the surface of 
PGMA-grafted fiber. They also performed similar surface graft polymerization 
onto a UHMPE fiber pre-treated with Ar plasma to improve the wettabihty and 
adhesion [127]. A corona-treated PE sheet was surface-grafted with GMA andsub- 
jected to the measurement of adhesive shear strength against a cured epoxy resin 
[128]. The adhesive strength of the grafted PE was nearly twice as high as that of 
the corona-treated PE. They ascribed the enhanced adhesion between the surface- 
grafted PE and the epoxy resin to covalent bonding of the epoxy groups on the 
grafted surface to the amine moieties in the epoxy resin (see Fig. 9). 




Fig. 9. Schematic representation of the structure of the interface between the GMA-graft poly- 
merized PE and the cured epoxy resin 
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4.1.2 

Adhesive Interaction in Aqueous Media 

Polymer surfaces modified by graft polymerization of water-soluble monomers 
underwent substantial adhesion to another surface when they were brought into 
contact in the presence of water, under pressure and subsequently dried [129]. 
The surface having a larger graft density generally exhibited stronger adhesion 
when the water present at the interface dried up. Substantial interaction occurred 
almost instantaneously upon contact when one surface was grafted with an anion- 
ic polymer and the other grafted with a cationic polymer. The interaction between 
similarly charged surfaces was weak in the presence of water, probably because 
of the electrostatic repulsion operating between the similarly charged groups. 
However, even a similarly charged film pair exhibited adhesive force when water 
was completely dried up, probably because of entanglement of graft chains (see 
Fig. 1 0). They concluded that the adhesion strength depended on the surface graft 
density, the microstructure of the grafted surface, the nature of the interfacial 
interaction, and the adhesion (drying) time. Similar adhesion characteristics 
were observed by Kang et al. [ 130] . 

Wong et al. [131] measured directly the interaction potential between a teth- 
ered ligand and its receptor in aqueous media. Using a surface force apparatus, 
the interaction force-distance profile was determined between streptavidin immo- 
bilized on a lipid bilayer and biotin tethered to the distal end of lipid-PEG. Both 
lipid bilayers containing streptavidin and biotin were absorbed onto the surface 
of mica having a specific curvature. Both cationic and anionic polymer grafted 
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Fig. 10. Schematic representation of interaction between surfaces with graft chains 
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films exhibited a significant and instantaneous adhesive interaction in aqueous 
media towards solid surfaces, such as surface-modified cellulose film, soda glass, 
and untreated PET film, so long as the countersurface had oppositely charged 
zeta potentials [132]. It should be noted that none of surface pairs exhibited a 
significant attractive force, regardless of the sign of their zeta potentials, unless 
their surface was not graft polymerized with anionic or cationic polymers, as can 
be seen in Table 1 . It maybe postulated that the diffused polyelectrolyte graft chain 
layer can easily and almost directly bind to the oppositely charged surface region 
of the opposed substrate. In other words, the surface which showed a strong attrac- 
tive force against a surface-grafted film may have a true surface potential the sign 
of which is opposite to that of graft chain. Detailed measurement of adhesion to 
a well-defined surface would make it possible to estimate the surface potential 
of solid materials not corresponding to the zeta potential. 

More detailed electrostatic interaction in aqueous media between two surfaces 
grafted with ionic polymer was measured with the use of atomic force microscopy 
(AFM) [133]. The AFM tip surface was modified by graft polymerization of a 
cationic monomer after being coated with a cyanoacrylate polymer. A PET film 
which was the counter-surface of the tip for the AFM measurement was also mod- 
ified by surface graft polymerization of ionic monomers. Appreciable adhesion 
was clearly observed in water between the grafted tip and the PET surface, while 
a repulsive interaction was noticed between the tip and the PET film when both 
were grafted with identically charged monomers. Addition of KCl to the medium 
exhibited a remarkable reduction in the adhesive interaction. These interactions 
were attributed to the Coulombic force between the grafted ionic polymer chains. 

4.2 

Medical Applications of Grafted Surfaces 

Most biomedical materials are used in constant contact with living systems, such 
as blood, cells, and tissues. Since the material surface can undergo unfavorable 
biological responses when in contact with a recipient, most of the conventional 
surfaces need to be modified so that the materials can function as designed. 

In the case of medical applications, some specific requisites usually dictate the 
modification technique. Among several methods investigated, grafting technique 
has advantages over others in several ways. They include easy and controllable 
introduction of newpolymer chains with a high surface density and precise local- 
ization of the chains at the surface, while keeping the bulk properties unchanged. 
Furthermore, covalent attachment of polymer chains onto a polymer surface can 
avoid their delamination in aqueous media, and thus ensures long term stabili- 
ty of the introduced chains, in contrast to physically coated polymer chains. All 
these aspects give the rationale for applying the grafting process for surface 
modification. Accordingly, medical applications of graft polymers have attract- 
ed considerable attention in recent years. 

In this section, we will highlight the use of the grafting technique for design- 
ing polymeric biomaterial surfaces that exhibit non-fouling property, selective 
protein adsorption, enhanced tissue adhesion, and minimum frictional damage 
to mucosa membranes. 
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4.2.1 

Non-Fouling Surfaces 

Prevention of protein adsorption is crucial to blood contacting devices includ- 
ing catheters, dialyzers, vascular grafts,blood containers, and oxygenators. When 
a material surface is brought into contact with blood, adsorption of serum pro- 
teins takes place as an initial event, which triggers successively the thromboge- 
nesis and complement activation cascade via the classical pathway. Since these 
unfavorable foreign-body responses would lead to serious clinical problems, 
much effort has been devoted to the creation of non-fouling surfaces. 

There is experimental evidence to support the versatility of surface grafting 
for minimizing protein adsorption as well as thrombogenesis. The grafting effect 
lasts for about a month or more, long enough for accomplishing the specific pur- 
pose of the devices, such as dialysis and oxygenation. However, the long-term 
durability of the grafting effect is still questionable [134, 135]. For example, vas- 
cular grafts and catheters for blood access require a non-fouling effect for a long 
duration of time if implanted permanently. 

For creating non-fouling surfaces, various water-soluble polymers have been 
used for surface grafting. They include nonionic, hydrophilic polymers such as 
PAAm, PDMAA, PEG, EVA, and PHEMA. A polymer containing phosphoryl- 
choline, a cell membrane constituent having zwitterions, was also used for the 
same purpose. Grafting of these polymers can be achieved by surface graft poly- 
merization [136, 137],coupling reaction [4, 138-142], surface segregation [143, 
144], and surface physical interpenetration [145]. Current advances in surface 
graft polymerization for non-fouling surfaces will be reviewed below. The details 
of other techniques are not within the scope of this section. 

Fujimoto et al. demonstrated, through an ex vivo adsorption experiment 
using radio -labeled immunoglobulin G (IgG),that IgG adsorption onto a PU film 
was remarkably reduced by surface grafting of PAAm through glow discharge 
treatment [146] and ozone oxidation [147]. Interaction of the PU surface with 
platelets was greatly reduced by this modification, when assessed by an ex vivo 
arterio-venous (A-V) shunt experiment in rabbits [146, 147]. Ruckert and 
Geuskens [148] reported that surface graft polymerization of PVP effectively 
reduced adsorption of fibronectin onto a styrene-(ethylene-co-butene)-styrene 
triblock copolymer film. It was further reported that the surface of PU catheters 
with tethered PDMAA chains remained unfouled even after 3 weeks of implan- 
tation in the rabbit vena cava inferior [149]. Recently, Kishida et al. [150] used 
reverse transcription-polymerase chain reaction (RT-PCR) to examine the 
expression of interleukin-l(3(IL-ip) mRNA secreted by macrophage-like cells 
(HL-60) cultured on a grafted PE surface as an index of inflammatory stimula- 
tion. They observed that cells cultured on the PAAm-grafted surface expressed 
a low level of IL- 1 (3 mRNA, indicating the non-fouling ability of the grafted sur- 
face. 

Rejection of protein adsorption to the outermost grafted surface is attributed 
to a steric hinderance effect due to the tethered chains. A grafted surface in con- 
tact with an aqueous medium, a good solvent of the chains, has been identified 
to have a diffuse structure [57, 151, 152]. Reversible deformation of the tethered 
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chains due to invasion of mobile protein molecules into the layer would lead to 
a repulsive force which is governed by the balance of entropic elasticity of the 
chains and osmotic pressure owing to the rise in the segment concentration. The 
overlapped repulsive force would prevent the direct contact of protein molecules 
with the substrate surface. 

It is interesting to note that the extent of rejection of protein adsorption is relat- 
ed to the polymer graft density [ 146]. As shown in Fig. 1 1, the amount of adsorbed 
IgG increased gradually beyond the threshold (approximately 1 0 pg/cm^), although 
the number of platelets adhered was kept at a low level. This result suggests that, 
when the graft layer is much larger in thickness than the protein radius and has 
a high water content, the protein molecules will readily migrate into the graft lay- 
er to be adsorbed, similar to protein molecules retained in the PAAm hydrogel in 
electrophoresis. Figure 12 illustrates schematically this protein sorption into a 
thick graft layer. Such sorption of plasma proteins into polymer matrices should 
be avoided, since the sorbed proteins would trigger the cascade of blood coagu- 
lation, similar to the plasma proteins in the stagnant, nonflowing blood. 

Extensive grafting gives rise not only to an increase in the thickness of graft 
layer, but also to an increase in the volume fraction of graft segments in the lay- 
er. The latter also has a significant influence on protein adsorption. The water con- 
tent of graft layer, which is a measure of volume fraction of graft segments in aque- 
ous media, can be varied by changing the condition of surface graft polymeriza- 
tion. A well-defined model experiment using several hydrogels made of PVA, 
PAAm, PVP, PEG, and PHEMA [153] provided much information for under- 
standing the interactions of grafted surfaces with blood components. The repre- 
sentative result obtained in this study is shown in Fig. 13. Apparently, all the hydro- 
gels have an optimum water content at 80-95 wt% to minimize platelet adhesion, 
regardless of the chemistry of network chains. A graft layer with a water content 
higher than this threshold would not reject platelet adhesion. Kulik and Ikada 




Fig.11. Adsorption of '^^I-labeled IgG and adhesion of platelet to the surface of PAAm-graft- 
ed PU films as a function of graft density: (A) IgG (virgin, ▲), (O) platelet (virgin, •) 
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Fig. 12. Interaction of blood with surface having grafted water-soluble chains 




Fig. 13. Effect of water content of hydrogels on platelet adhesion; PVA (□), PAAm (O), methoxy- 
PEG methacrylate (A), PVP (•), and PEG (A) 
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[153] attributed the phenomenon to the entrapment of platelets in the partially 
broken network structure of large pore size present at the surface of hydrogels. In 
contrast, when the water content of the hydrogels is below 80-95 wt%, the pene- 
tration of platelets through the hydrogel network seems to be greatly limited due 
to the high volume fraction of network chains. When the water content of the 
hydrogel becomes lower than 80 wt%, attraction of platelets by the outermost sur- 
face of the hydrogel will predominate probably due to the large free energy at the 
interface between the aqueous phase and the hydrogel surface. It was further 
demonstrated that carboxylic acids fixed to the hydrophilic polymer network of 
hydrogels served to reduce the number of platelets adhered. In contrast, cationic 
groups were found to have opposite effects on platelet adhesion [ 153] . This attrac- 
tive interaction of cationic surface is operative not only for cell adhesion [56, 57] 
but also for protein adsorption [154]. These ionic effects may have relevance to 
ionic charges distributed over the surface of cells and protein molecules. As 
demonstrated by Kato et al. [154], protein adsorption onto polymer surfaces 
grafted with highly ionic polymer chains is primarily governed by the electro- 
static attraction and repulsion between the charged surfaces and charged proteins. 
The ionically grafted polymer surface prevents adsorption of proteins with sim- 
ilar charge, but accelerates adsorption of proteins with opposite charge. 

PEG molecules which are relatively nontoxic and capable of reducing the 
interactions between blood components and man-made materials, can also be 
tethered to a polymer surface through surface grafting. This has been achieved 
by free-radical polymerization of methacrylate monomers carrying a pendant 
PEG chain [47, 155-158]. The surface of a PU film was subjected to UV-induced 
graft polymerization of methoxy-PEG methacrylate monomers with 4, 9, and 
23 units of ethylene glycol (EG) [ 156]. As shown in Fig. 14, the monomer with the 
shortest PEG length of only 4 EG units was most inert toward the blood compo- 
nents when the graft yield was largely reduced by the use of high concentrations 




Fig. 14. Platelet adhesion to the PU films grafted with methoxy-PEG methacrylate at different 
concentrations of chain transfer agent. The number of EG units of the monomer; 4 (•), 9 (O), 
and 23 (□) 
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of a chain transfer agent. In addition, extraordinarily high graft yields were not 
effective in preventing protein adsorption, similar to the case of very low graft 
yields which were insufficient to reduce the interactions with proteins and 
platelets. Plasma-induced surface graft polymerization of PEG methacrylate was 
carried out onto bioresorbable poly(D,L-lactide) films [158]. Fibrinogen adsorp- 
tion and platelet adhesion decreased upon grafting. Nevertheless, there was no 
significant effect on protein adsorption and platelet adhesion, so far as the num- 
ber of EG units was maintained from 1 to 10. 

In an attempt to synthesize hemocompatible dialysis membranes, a phos- 
phorylcholine-containing methacrylate polymer was grafted to a cellulose mem- 
brane in a heterogeneous system using cerium ammonium nitrate [159, 160]. 
Kang et al. [161] and Qiu et al. [47] reported further immobilization of heparin, 
a naturally-occurring anticoagulant, onto surfaces grafted with acryloylbenzo- 
triazole and ethylene glycol methacrylate polymers, respectively. It seems that 
the heparin-immobilized surfaces inhibit thrombus and fibrin formation by 
binding to the proteins involved in the clotting cascade. 

As stated above, the non-fouling property is primarily based on the 
hydrophilicity of the graft chains. In this respect, multifunctionally hydroxylat- 
ed PAAm, in which up to six hydroxyl groups were attached to each of the side 
chain terminal [162], and a glucose-containing methacrylate polymer, in which 
glucose residues were carried by each side chain through an ester bond [163], 
will be interesting candidates as these graft chains and capable of generating non- 
fouling surfaces. 

4.2.2 

Physiologically-Active Surfaces 

Various reactive groups introduced onto polymer surfaces by direct grafting of 
functional polymers [73, 164, 165] or post-derivatization [166] of graft chains 
can be used to bind biologically active macromolecules, including proteins, poly- 
saccharides, and deoxyribonucleic acids, onto the substrate surface. For exam- 
ple, primary amines, derivatized from amide groups of PAAm graft chains via 
Hofmann degradation as shown in Fig. 15, can be utilized for chemical immobi- 

PAAm-CONHs 

i +CIO- 

e 

PAAm-CON-CI 

I 

PAAm-N=C=0 
I +OH- 

PAAm-NH 2 

Fig. 15. Derivatization of the amide group of a PAAm graft chain into primary amine via Hof- 
mann degradation 




Surface Modification of Polymers by Grafting 



29 



lization of proteins onto a substrate [166], The amount of covalently immobilized 
proteins increased almost proportionally with the increase in graft yield, sug- 
gesting that protein molecules penetrated into and immobilized in the graft lay- 
er [ 164] . Such a surface possesses specific bioactivity associated with the immo- 
bilized biomacromolecules, allowing highly specific molecular recognition to be 
achieved. An example is selective adsorption of pathogenic proteins to the mod- 
ified surface from blood. This selective adsorption has an important application 
as immunoadsorbents. 

Immunoadsorption, an advanced therapeutic modality, focuses on detoxifi- 
cation of patient blood rich in high-molecular -weight pathogenic substances, 
mostly abnormal autoantibodies such as rheumatoid factors in rheumatoid 
arthritis (RA) and anti-DNA autoantibodies in systemic lupus erythematosus 
(SLE). Detoxification of these pathogens will be accomplished through extra- 
corporeal perfusion of the patient plasma or whole blood over an affinity col- 
umn made of immunoadsorbents. These adsorbents perform their function 
through the same mechanism as conventional affinity adsorption, where proteins 
in the liquid phase are adsorbed on the specific ligands immobilized onto an 
insoluble support. 

A major problem associated with the current immunoadsorption is the low 
capacity of adsorbents, which can probably be attributed to the materials used 
as a solid support matrix. To solve this problem, an attempt was made to syn- 
thesize immunoadsorbents utilizing a solid support made of super fine PET 
microfibers [167, 168]. The use of such a fibrous support has great advantages 
over the conventional matrices, because this fiber is very large in specific surface 
area, excellent in mechanical strength, and biosafe. 

The feasibility study of the immunoadsorbent synthesized through grafting 
of PAAc chains onto the PET fiber and the subsequent protein immobilization 
was conducted [167]. Three kinds of serologically active ligand proteins, includ- 
ing protein A, a specific antigen, and a specific antibody, were bound to the sur- 
face of PAAc-grafted PET fibers using l-ethyl-3-{3-dimethylaminopropyl) car- 
bodiimide (EDC). These adsorbents were tested for selective removal of IgG, spe- 
cific polyclonal antibodies, and specific antigens. The binding capacity and 
specificity of the adsorbents were also evaluated in phosphate buffer solution and 
serum. The results are shown in Table 2, which provides useful information for 
applying immunoadsorbents to clinical systems. 

An immunoadsorbent applicable to blood detoxification in patients with SLE 
has been synthesized by the grafting technology utilizing PET microfibers [168]. 
Double-stranded DNA was chosen as the ligand to be immobilized because the 
circulating pathogen in this case is polyclonal antibodies specific to DNA. The 
work was devoted to developing a method for DNA immobilization onto the 
grafted surface, because few publications have reported covalent immobilization 
of native DNA directly onto polymeric supports. Ample and stable immobiliza- 
tion of calf thymus double-stranded DNA, a highly anionic macromolecule, was 
carried out onto the PET fiber surface following the tethering of cationic graft 
chains with tertiary amines. A large capacity and a high specificity of anti-DNA 
antibody binding were found for the sera obtained from lupus mice and a human 
patient. 
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4.2.3 

Slippery Surfaces 

Polymer surfaces modified by immobilization of highly-hydrophilic graft chains 
exhibit lubrication in the hydrated state [169-173] . This property has spurred the 
effort to make the surface of tissue-contacting tubular devices such as catheters, 
cannulae, endoscopes, and cystoscopes lubricious. Most of the outer surfaces of 
body orifices, such as mouth, esophagus, nose, urethra, and vagina, into which 
tubular devices are often inserted are covered with a mucosa layer. This layer 
consists mainly of hydrophilic glycoproteins producing highly viscous sub- 
stances overlaying the endothelial cell basement. Lubricating the surface of the 
tubular devices will enable painless insertion, precise operation, and protection 
of the tissues from injuries. 

Following ozone oxidation, the surface of a PU film was graft -polymerized 
with DMAA and the coefficient of kinetic friction (pk) for the fully hydrated, 
grafted films of two different graft densities was determined against a cleaned 
steel plate in distilled water as a function of the sliding velocity [ 174] . It was found 
that grafting of PDMAA effectively reduced the frictional force. 

The same strategy was adopted to make the surface of a cystoscope lubricious 
[ 175] . It was covered with a PDMAA-grafted PU sheath of the same diameter as 
the cystoscope and the efficacy of the lubricated cystoscope was evaluated using 
rabbits by an in vivo test simulating cystoscope operation. As shown in Fig. 16, 
the maximal resistance force on the cystoscope was decreased by grafting. A his- 
tological study proved that urethral damage caused by rubbing against cysto- 
scope was effectively reduced by this lubrication technique. 




Cycle of reciprocating movements 



Fig. 16. Change in maximal resistance force for the cystoscope model during in vivo test. Per- 
pendicular bar shows standard deviation 
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4.2.4 

Tissue-Adhesive Surfaces 

4.2.4.1 
Soft Tissues 

It is often demanded that the surface of polymeric biomaterials should exhibit 
permanent tenacious adhesion to soft connective and dermal tissues. However, 
conventional non-porous, polymeric materials will be encapsulated by a fibrous 
membrane generated de novo by surrounding fibroblasts, when subcutaneous- 
ly implanted into the living body in contact with soft connective tissues. This is 
a typical foreign body reaction of the living system to isolate foreign materials 
from the host inside the body. On the other hand, it should be noted that the small 
gap present between a percutaneously-implanted device and the surrounding tis- 
sue provides a possible route for bacterial infection because of the lack of micro- 
scopic adhesion at the interface. 

In contrast, when collagen, which is the most abundant extracellular compo- 
nent of connective tissues, is covalently immobilized onto an implant surface uti- 
lizing reactive groups introduced by surface graft polymerization and then 
implanted into the body, a fibrous membrane formed strongly adheres to the col- 
lagen-immobilized implant surface [176, 177]. Consequently, there remains no 
gap between the implant and the soft tissue. It was also reported that collagen 
immobilized onto a PE surface through PAAc graft chains had an inhibitory 
effect toward the tumor promoting activity [178]. 

Following the corona-discharge treatment of a silicone surface and the sub- 
sequent graft polymerization of AAc, type I atelocollagen was immobilized onto 
the grafted surface with the use of water-soluble carbodiimide [176, 177]. As 
depicted in Fig. 17, the immobilization reaction involves two steps, i.e., activa- 
tion of carboxylic acids and the following nucleophilic substitution with prima- 




EDC 



I +(CH3)2N-(CH2)3-N=C“N-CH2-CH3 





collagen 






Fig. 17. Collagen Immobilization onto the PAAc-grafted polymer surface with the use of car- 
bodiimide 
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ry amine moieties carried by collagen molecules. This method facilitated ample 
immobilization of collagen up to the order of 10 pg/cm^ under mild conditions 
in aqueous media. This amount seems to be sufficient for the complete coverage 
of an implant surface with collagen. The feasibility study of this two -step metho d 
for covalent immobilization of collagen onto a PAAc-grafted surface was recent- 
ly conducted by Lee et al. [179]. A simpler synthetic route for fabricating colla- 
gen-immobilized devices has been proposed recently [180]. The idea is based on 
the electrostatic attraction between different polyelectrolytes, through which 
collagen bearing a net positive charge forms a complex with an anionically graft- 
ed, polymeric surface. This method also provides stable immobilization of col- 
lagen. 

The force required to pull out a collagen-immobilized silicone device was 
measured after three weeks of percutaneous implantation of this device into the 
rabbit back [177]. The results revealed that collagen immobilization remarkably 
enhanced the adhesion strength at the interface. Histological observation 
revealed that the immobilized collagen did not activate epidermal down-growth, 
but promoted the intimate adhesion at the material-tissue interface, although a 
relatively severe inflammatory reaction was observed in the initial stage of 
wound healing. As shown in Fig. 1 8, firm adhesion mediated by the immobilized 
collagen plays an important role as a barrier against infection. 

The rationale for immobilizing collagen onto a biomaterial surface for its soft 
tissue adhesion is given by the fact that fibroblasts proliferate under either direct 
binding to an RGB (Arg-Gly- Asp) motif carried by the collagen molecule or indi- 
rect attachment mediated by tissue fibronectin bound to collagen in the living 
connective tissue. Plausibly, a collagen-immobilized polymer surface serves as a 
bed, similar to the natural extracellular matrix. Empirical confirmation of this 
mechanism was given by Tamada and Ikada [ 1 81 ] . The immobilized collagen was 
shown to promote the in vitro proliferation of fibroblasts and its metabolic activ- 
ity in the initial stage of subculture. 




Fig. 18. Bacterial infection observed during percutaneous implantation of silicone devices into 
rabbits; (O) non-treated and (•) collagen-immobilized 
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Collagen immobilization onto a silicone rubber surface which underwent 
plasma-induced graft polymerization of PHEMA was also reported to have 
improved attachment and growth of corneal epithelial cells onto the rubber sur- 
face [48,50,182]. 

4.2.4.2 
Hard Tissues 

A flexible polymeric material having an ability to bond to hard tissues was syn- 
thesized through plasma-induced graft polymerization of a phosphate-contain- 
ing monomer, methacryloyloxyethylene phosphate [183]. Potential applications 
of such a material are in artificial ligaments, tendons, intervertebral discs, and 
periodontia where strong adhesion to bones is needed. The dihydrogen phos- 
phate of the monomer has an affinity for hydroxyapatite (HAP) which is the main 
inorganic component of bone tissues. When immersed in a simulated body solu- 
tion of pH 7.4 containing calcium and phosphate ions at their supersaturated con- 
centrations with respect to HAP, a polymer surface having the phosphate-con- 
taining chains promoted deposition of a thin HAP layer onto its surface. Figure 1 9 
illustrates the schematic model proposed for deposition of an HAP layer. It was 
further demonstrated that bonding of the deposited HAP layer to the substrate 
surface was significantly improved by grafting of the phosphate-containing 
polymer chains. The study in simulated body environment led to the conclusion 
that the grafted surface has potentiality to bond firmly to bone tissues. 

As in the case of strong fixation, rapid establishment of bone-implant adhe- 
sion is an important factor influencing the clinical consequences. It seems that 
this effect can be achieved by accelerating the HAP crystallization at the implant 
surface upon implantation. In this connection, an HAP layer was deposited in 
vitro onto a polymer surface grafted with the phosphate-containing polymer pri- 
or to implantation [ 184] . The pre-deposited HAP layer was fabricated by a solu- 
tion-mediated wet method. The pre-deposited HAP served as a seed for HAP 
crystallites to grow rapidly, consuming calcium ions from the surrounding aque- 
ous phase at a rate of 7 pg/cm^/day, when brought into the physiological envi- 
ronment. Furthermore, enhanced calcium deposition at the HAP pre-deposited 




Fig. 19.a Proposed structure of the HAP firmly attached to the polymer surface having phos- 
phate-containing graft chains, b Binding between the graft chains and the HAP crystal 
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Implantation time (Week) 

Fig. 20. Percentage of the PE rod surfaces in direct contact with the bone substrate established 
during the implantation into rat femora: (O) non-treated, (□) with phosphate-containing graft 
chains, and (•) with the HAP pre-deposited layer. Errorbars represent standard deviation 



surface was also demonstrated in osteoblast cell culture [185]. The bone-bond- 
ing ability of polymer surfaces having phosphate-containing graft chains and a 
pre-deposited HAP layer was evaluated in vivo [186]. The surface-modified spec- 
imens were prepared using a PE rod and PET fabrics as substrates and implant- 
ed into the rat femur. As shown in Eig. 20, the results confirmed the bone-bond- 
ing ability of these surfaces and the accelerated HAP apposition by the pre- 
deposited HAP. 
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1 

Introduction 

Crystalline polymers generally comprise a variety of phase structures, including 
the crystalline and noncrystalline components. The noncrystalline component 
is thought to be in a supercooled state of the molten state of polymers, either in 
the rubbery or glassy state. The question is whether a somewhat ordered non- 
crystalline component exists or not, due to the coexistence with the crystalline 
component. E lory reported in 1 949 that the boundary between the crystalline and 
amorphous regions of most long-chain molecules will not be well-defined as is 
typical of monomeric systems [1]. Subsequent theoretical analyses, involving 
several methods, have quantitatively established the existence of an interfacial 
region that comprises a transition phase from the crystalline to the amorphous 
phase [2-11]. Over the past few years, a variety of experimental methods have 
confirmed the expectation of such an interphase [12]. These methods involve 
broad-line NMR [13-16] high-resolution solid-state ^^CNMR [17-31] Raman 
spectroscopy [32-35] , and small-angle X-ray and neutron scattering, among oth- 
ers [36-42]. 

To elucidate the phase structure in detail it is necessary to characterize the 
molecular chain conformation and dynamics in each phase. However, it is rather 
difficult to obtain such molecular information, particularly of the noncrystalline 
component, because it is substantially amorphous. In early research in this field, 
broad-line NMR analysis showed that linear polyethylene crystallized from 
the melt comprises three components with different molecular mobilities; solid, 
liquid-like and intermediate molecular mobility [13-16]. The solid component 
was attributed to molecules in the crystalline region, the liquid component to 
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that in the noncrystalline amorphous and the intermediate component to that in 
an intermediate region between the crystalUne and amorphous regions. A simi- 
lar approach was also used on samples crystallized from dilute solution and 
some characteristics of the phase structure of the samples was elucidated. How- 
ever, since broad-line NMR can characterize substances only by overall mol- 
ecular mobility, further detailed information with regard to molecular chain 
conformation and alignment in each component could not be obtained. Hence, 
some objection to the broad-line *H NMR analysis could remain [49]. 

Rather recently, we have studied the solid-state structure of various polymers, 
such as polyethylene crystallized under different conditions [17-21], poly 
(tetramethylene oxide) [22], polyvinyl alcohol [23], isotactic and syndiotactic 
polypropylene [24, 25], cellulose [26-30], and amylose [31] with solid-state high- 
resolution NMR with supplementary use of other methods, such as X-ray dif- 

fraction and IR spectroscopy. Through these studies, the high resolution solid- 
state NMR has proved very powerful for elucidating the solid-state structure 
of polymers in order of molecules, that is, in terms of molecular chain confor- 
mation and dynamics, not only on the crystalline component but also on the non- 
crystalline components via the chemical shift and magnetic relaxation. In this 
chapter we will review briefly these studies, focusing particular attention on the 
molecular chain conformation and dynamics in the crystalline-amorphous 
interfacial region. 



2 

Solid-State NMR [43] 

2.1 

High-Resolution NMR Spectrum 

In this chapter we review the study of the solid structure of many crystalline poly- 
mers mainly with high-resolution solid-state NMR, so we will briefly sum- 
marize here its principles for the convenience of the reader. For this purpose we 
first consider the Hamiltonian of an ensemble of nuclei possessing spin in a sta- 
tic field Bq. The Hamiltonian to be considered of this spin system can be written 
as: 

H = Hz + Hs + Hd + Hj 

Here Hz represents the direct interaction energy of spins with the external mag- 
netic field Bq. This term, the so-called Zeeman energy, gives the main resonance 
line as a 6-function for each nucleus with spin. Hs, chemical shift Hamiltonian, 
describes the indirect interaction of spins with Bq via electrons, and gives differ- 
ent shifts to the main resonance peak that can distinguish each constituent nucle- 
us in a molecule by difference of the electromagnetic environment. Ho and Hj 
describe the direct and indirect (via electrons) dipolar interactions between 
spins, respectively. The former gives the line width to the main resonance line and 
the latter, the so-called J-coupling or spin-coupling, brings splitting or satellites. 
In addition to these terms there are some other terms, e.g. quadrupole Hamil- 
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tonian. In this chapter we only consider NMR for spin systems whose nuclei 
have spins of either zero or 1/2 and, therefore, such further terms are omitted. 

In or NMR measurements in solution, the direct dipolar interaction Hp 
actually disappears because, due to rapid molecular motion, the interspin (inter- 
nuclear) vectors are rapidly space-averaged within the time scale of the mea- 
surements. Hence, Hs and Hj are detectable as sharp lines or splittings in a high- 
resolution or spectrum and they can be related to the detailed molecular 
structure or conformation of the substance being investigated. In a solid, how- 
ever, the directions of the internuclear vectors are stationarily fixed in time even 
if they are distributed randomly in space. Then, Ho gives a very wide line width 
to the main resonance line and completely masks all lines due to H$ and Hj. There- 
fore, in usual NMR measurements, one observes only a very broad resonance 
spectrum which is determined by Hz and Ho- Some information is, of course, 
obtainable, even if the spectrum is very broad,because the line shape that is deter- 
mined by Ho reflects the geometrical alignment of nuclei and the line width 
reflects the degree of molecular motion. However, a spectrum in which the terms 
due to Hs are detectable cannot be obtained and, therefore, is termed the broad- 
line NMR spectrum. 

To obtain high-resolution spectrum for solid matter, we have to eliminate the 
effect of Ho and furthermore average that of Hs, because Hs is a tensor of rank 
2. In the solid-state NMR, the high-resolution spectrum is obtained by com- 

bining three techniques, ^H dipolar decoupling (DD), magic angle spinning 
(MAS), and cross-polarization (CP). DD is a technique to eliminate the effect of 
Ho, which is achieved by applying a resonant oscillating field Bi of sufficient 
amplitude on ^H spins. If the resonant oscillating field is applied, the dipolar 
field between and ^H and between *H spins themselves now changes the 
direction rapidly, so that the dipolar interactions Ho between and ^H, and 
between ^H spins themselves, are removed. Furthermore, the interaction 
between nuclei is negligible because of their low natural abundance (ca. 
1.1%). Hence the magnetization, without any contribution from Ho > can be 
detected. With use of ^H dipolar decoupling (DD), the spectrum which contains 
only the contribution from Hz and Hs can be obtained. The spectrum thus 
obtained can be called the Hs spectrum, i.e. chemical shift spectrum. However, 
high-resolution spectrum still cannot be obtained in general,because the chem- 
ical shift Hamiltonian Hs has anisotropy as a tensor quantity of rank 2 as men- 
tioned above. The elimination of this anisotropy is usually achieved by rotating 
samples rapidly around an axis inclined at an angle of 54° 44' to the static field 
Bo during measurement. Here 54° 44' is the so-called magic angle that fulfills the 
relationship: 



1 - 3cos^0 = 0. 

This process is called magic angle spinning (MAS). If the rate of spinning is larg- 
er than the line width caused by the anisotropy of H$ in the frequency scale, the 
anisotropy is removed and a sharp resonance line appears at the average of the 
principal values of Hs as a tensor quantity. Thus, a high-resolution NMR spec- 

trum can be obtained by using two procedures, DD and MAS. 
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In addition to MAS and DD explained above, cross polarization (CP) is usu- 
ally used in present solid-state high-resolution NMR measurements. This is 

a process to enhance magnetization by cross polarization (CP) between ^^C 
and spins. The ^^C magnetization is enhanced by the equilibrium magneti- 
zation of abundant spins. This procedure is achieved by contacting the ^^C 
spin system to the system while Zeeman energies of both ^^C and spins 
against the respective sub-field Bj are equalized. A '^C magnetization as large 
as 4 times its equilibrium magnetization is obtained by this process. Further- 
more, in this case, the time demanded for one measuring pulse sequence is 
determined by the spin-lattice relaxation time Tih of 'H that is usually much 
shorter than the Tic of ^^C. Hence, the time required for obtaining a spectrum 
with a sufficient signal-noise ratio can be largely shortened by a factor of 1/100 
or sometimes 1/1000 by using a CP in comparison to the measurements with- 
out CP. However, it should be noted that the CP efficiency may differ according 
to the molecular mobility of samples. Sometimes, we could not use CP for the 
purpose of characterizing different phases with different molecular mobilities 
in samples. 

In Fig. 1, the pulse sequences frequently used in this work are schematically 
shown. The pulse sequence I is used for obtaining DD/MAS ^^C NMR spectrum 
without the CP procedure. The *^C magnetization which appears in the Bq direc- 
tion for a time of tc is turned through 90° by a resonant field and the free indue- 



I. Single Pulse Sequence for DD/MAS 





II. Pulse Sequence for CP/M AS 




lU. Pulse Sequence for Longitudinal and Transverse Relaxation 

90’ pulse train gp* 

n n 



Fig. 1 . Pulse sequences for high-resolution NMR. I. Single pulse sequence for DD/MAS with- 
out CP. II. Pulse sequence for CP/MAS. III. Pulse sequence for longitudinal and transverse relax- 
ation. See the text for detail 
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tion decay (FID) is observed in the xy plane (perpendicular to Bq) under 
decoupling (DD ) with MAS throughout the measurement. Here, if the waiting time 

that is actually the repetition time of the pulse sequence, is set to be as long as 
5 times the Tic, thermal equilibrium spectrum can be obtained [44]. Here we 
sometimes express this pulse sequence as (90 °-FIDdd- i£)n , where n is the rep- 
etition number of the measurement. Sequence II is that required to obtain the 
CP/MAS NMR spectrum. The *H magnetization which appears in the Bq 
direction for a time of is turned through 90° and locked in the ^H resonant field 
and the magnetization obtained by the cross polarization with the ^H system 
is observed under DD. The cross polarization is performed by contacting both 
spin systems by applying the resonant field while Zeeman energies of both ^H 

and spin systems in relation to respective resonant fields are equalized. 

In addition to pulse sequences I and II, various pulse sequences were used in 
this work. Figure 1, shows a pulse sequence for obtaining the partially relaxed 
spectrum in longitudinal and transverse directions (Pulse sequence III). After 
saturating by a 90° pulse train, the magnetization appearing longitudinally 
is turned 90° and allowed transverse relaxation for r, and the FID is observed 
under DD. Using this pulse sequence, the longitudinal and transverse relaxations 
of the different components, particularly shorter Tic components in polymers, 
can be examined. 

2.2 

Spin Relaxation 

High-resolution NMR spectra give detailed information of molecular align- 
ment of polymers via chemical shift that distinguishes each carbon in the mol- 
ecule. Furthermore, studies of the magnetic relaxation provide information 
about molecular motion with relation to each carbon. In the case of organic sub- 
stances and polymers that contain only *H and as spin having nuclei, the 
relaxation is conducted by the time fluctuation of the magnetic dipolar interac- 
tions between and ^H spins. Hence, it provides information about molecular 
motion via the time fluctuation of interspin vector (i.e. internuclear vector), con- 
necting and ^H nuclei. 

The relaxations in the directions parallel and perpendicular to Bq, i.e. the lon- 
gitudinal and transverse relaxations, are often measured in present NMR. The 

relaxation parallel to Bq is called the longitudinal relaxation or spin-lattice relax- 
ation, or simply Ti relaxation, since it is carried out by energy exchange between 
the spin system and the lattice. On the other hand, the relaxation perpendicular 
to Bo is called transverse relaxation or spin-spin relaxation, or simply T 2 relax- 
ation, since it is conducted by the spin exchange inside the spin system. Hereafter, 
we designate the spin-lattice and spin-spin relaxations of the magnetization 
as Tic and T 2 C, and those of ^H as Tin and T 2 H. respectively. These are given as 
functions of the spectral density which is defined as the Fourier transform of the 
time-dependent part of the dipolar interaction. Since this time-dependent part 
is a function of the orientation of the internuclear vector, it is generally called the 
orientation function. That is Tic and T 2 C are given by the spectral densities, 
J (fi>)m>that are defined as Fourier transforms of the correlation functions of the 
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different terms of the orientation function, and can be defined as in Eqs. ( 1 ) and 
(2), respectively. 

1/Tic = (1/20)K [Jo(tt>H-<»c) + 3Ji(n}c) + 6 J 2 (tUH+ft>c)] (1) 

I/T 2 C = (1/40)K [4Jo(0) + Jo(tt)H-<»c) + 6 J 1 (6}H)+3Ji(mc) + 6 J 2 (n>H+ft>c)] (2) 

Here the spectral densities are defined as: 

J(®)m = J exp(-i®T)dT/ 1 (0)^ | 

Fm(t) are different terms, designated by suffix m, of the orientation function in 

the dipolar interaction and (t) (0) ®re their correlation functions, coh and 

coq are the resonant frequencies (Larmor frequencies) of 'H and They are 
defined by magnetogyric ratios, yu and yc, of *H and and the field intensity 
Bo as oiH = -ThBo, «c = -TcBo. 



K is a constant defined as 

K = (n,/47tfrly}:nr-^ 



where fiQ is the magnetic transmittance of vacuum and h Planck's constant/27T. 
r is the internuclear distance between and *H. If the correlation functions 
evolve exponentially with one time constant as Tc as 

fm(T)F„,(0) = |F(0)'| exp(-|T|/Tc) (3) 

the longitudinal and transverse relaxation times are given as 



1/Tic =(1/20)K 
(4) 
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where is the relaxation time of the motion. The theory with relation to the 
motion dictated by only one relaxation time (correlation time) is called single 
correlation-time theory. 

As revealed from Eqs. (1) and (2), or their candid forms (4) and (5), the lon- 
gitudinal relaxation is determined by the spectral densities in the order of can 
coc, whereas the transverse relaxation involves the contribution from the zero 
frequency component Jo(0). In the case of solid matter, rc is generally very long. 
Hence, the transverse relaxation is predominantly determined by the zero fre- 
quency component Jo(0). In Eq. (5), for example, the zero frequency term (the 
first term) dominates the other terms that are reciprocally proportional to tq for 
oo^T^ >> 1. Tic increases as tc increases (i.e. as the material under consideration 
becomes solider), whereas T 2 C decreases infinitely as rc increases. Eor example. 
Tic is generally in an order of several tens ~ several hundreds of seconds for the 
crystalline component and in an order of a few tenths of a second for rubbery 
components of polymers. On the other hand, T 2 C is of an order of a few tens of 
microseconds for the crystalline or glassy component and a few milliseconds for 
the rubbery component of polymers. In this work. Tic and T 2 C are used for char- 
acterizing different components in crystalline polymers. 



3 

Linear Polyethylene Crystallized from the Melt 
3.1 

Introduction and Approach by Broad-Line NMR 

Linear polyethylene in the solid-state has been most extensively studied as a typ- 
ical representative of a crystalline polymer by various methods because of the 
simplicity of its molecular and crystal structure. Figure 2-A shows a broad-line 
NMR spectrum in the first derivative form of linear polyethylene at 20 °C. 
The sample is a molecular weight fraction with a viscosity-averaged molecular 
weight of 90,000, which was crystallized isothermally at 130 °C from the melt. 
As pointed out in the last section, one can only see in this broad-line NMR 
spectrum the contribution of the dipolar interaction Ho between the spins 
themselves, since all contributions from H$ or Hj that provide detailed infor- 
mation about the molecular alignment in the sample are masked by Ho- How- 
ever, the broad line shape narrows and approaches a Lorentzian function as the 
molecular mobility increases. Hence, one can characterize the structure of the 
sample by examining the line shape and line width in terms of molecular 
mobility. At first glance the sample may comprise two different components 
with different line widths, one being a rather broad Gaussian and the other a 
narrow Lorentzian function, that may correspond to the crystalline and amor- 
phous phases, respectively. However, if this kind of spectrum is analyzed into 
two components, assuming a two-phase structure of crystalline and amorphous 
phases as shown in Fig. 2-B, (either straight line decomposition [45] or sym- 
metrical decomposition [46,47] the broad component exceeds the crystalline 
fraction determined by X-ray diffraction analysis or density measurements. 
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Fig. 2. a 60 MHz broad-line *H NMR spectrum of bulk polyethylene at 20 °C. b Two-component 
analyses, dotted line: straight line decomposition [45] dashed line: symmetrical decomposi- 
tion [46, 47] 



This suggests that the broad component includes some contribution from a 
rigid noncrystalline component. 

Bergmann et al. broke down such spectra for polyethylene into three compo- 
nents, i.e. broad, intermediate, and narrow components [13, 14]. The methylene 
groups of the polymer were divided into three classes: rigid, hindered-rotation- 
al, and micro -Brownian mobile methylene groups, and these were considered to 
contribute to the broad, medium, and narrow components, respectively. The 
mass fraction of the broad component thus estimated was in good accord with 
the crystalline fraction obtained from X-ray diffraction analysis and density 
measurements. The hindered-rotational and micro -Brownian mobile methylene 
groups were assigned to two noncrystalline phases, the former with less mobil- 
ity and the latter with a liquid-like mobility. The relaxation processes, called a, (3 
and Y> observed in the dielectric and mechanical measurements were well 
explained as corresponding to the processes associated with a local molecular 
movement in the crystalline region, and micro-Brownian and limited molecular 
motions in the noncrystalline region, respectively. 

We have examined linear polyethylene samples, covering a very wide range of 
molecular weight, that were crystallized from the melt by a similar three-com- 
ponent analysis of the broad-line ‘H NMR spectra and obtained very important 
information about the phase structure of the samples [15,16]. Figure 3 shows the 
molecular-weight dependence of the mass fractions of the three components at 
room temperature. As can be seen, although the narrow component is negligible 
for samples of molecular weights smaller than 31,800 it appears clearly at a mol- 
ecular weight of 44,900 and increases as the molecular weight increases. The 
mass fraction of the medium component is as small as 0.05-0.08 for samples 
whose narrow component is actually devoid but it gradually increases with the 
appearance of the narrow component and reaches a rather high level for samples 
with higher molecular weights. In other words, the samples with very low mole- 
cular weight are only made up of crystalline lamellae and noncrystalline inter- 
lamellar material without a liquid-like amorphous component. The amorphous 
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Fig. 3. Mass fractions of crystalline, rubbery, and crystalline-amorphous interfacial phases of 
bulk polyethylene as a function of molecular weight Mv. O: data by broad-line *H NMR analy- 
sis. •: data by high-resolution ’^C NMR analysis 

component appears after the mass fraction of the interlamellar material reach- 
es a chosen level. 

These studies provide important information about the interlamellar mater- 
ial of crystalline polyethylene as a function of molecular weight. The samples with 
molecular weight less than about 30,000 are composed of the crystalline and 
interfacial regions without an interzonal region with liquid-like character. With- 
in this molecular weight range, the average ratio of the lamellar thickness to the 
extended molecular chain length, ^/x, changes from unity to about 1/4 [48]. 
Hence, a molecular chain will participate in a crystal lamella almost always only 
once, or at least less than three times, so that molecular chain folding will be very 
scarce in forming the lamellar structure. The noncrystalline molecular chains 
that are excluded from the crystalline region will be somewhat severely restrict- 
ed by the crystalline molecular chains in conformation and mobility. The non- 
crystalline component with random conformation with liquid-like mobility will 
not appear. 

As the molecular weight increases from 45,000 to 100,000, the ratio ^/x varies 
from 1/4 to 1/7. The number of molecular chains that penetrate several crystal 
lamellae, or participate repeatedly in a crystal lamella by molecular chain fold- 
ing, will gradually increase. The amorphous region, in which the molecular con- 
formations are distributed over all allowable conformation and rapidly trans- 
formed each other, can appear in this molecular weight range. The transition 
region from the crystalline to this amorphous region (the crystalline-amorphous 
interfacial region) will comprise molecular chains whose conformation and 
motion are severely restricted, including folded molecular chains, entanglement 
and molecular end groups. However, since broad-line ^H NMR can characterize 
the phase structure only in terms of the molecular mobility through the line shape 
and line width or the second moment of the spectrum, it is difficult to obtain fur- 
ther detailed information about molecular alignment or conformation and 
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hence some objection to our analysis could remain [49] . However, we were able 
to obtain more detailed information about the phase structure, not only for these 
samples but also for many crystalline polymers, by use of solid-state high-reso- 
lution NMR [17-31]. We will discuss the above-cited results obtained by 
NMR later. 

3.2 

Approach by High-Resolution NMR 

High-Resolution NMR Equilibrium Spectrum. Figure 4 shows the DD/MAS 
NMR spectrum at room temperature for a linear polyethylene of molecular 
weight of 3.0 x 10® that was crystallized isothermally at 130 °C from the melt. 
This spectrum was obtained at a field strength of 4.7 T (Larmor frequency of 
ooc = 50 MHz) by a single pulse sequence (ti/2 - FIDdd - 1700s)n that is schemat- 
ically shown in Fig. l-I. In this pulse sequence, since r^ was set at 1700 s longer 
than 5 times the longest Tic of this sample, the spectrum expresses the thermal 
equilibrium state. In the spectrum two distinct resonance lines can be recognized 
at 33 and 31 ppm relative to the chemical shift of tetramethylsilane (TMS); des- 
ignated hereafter as resonance lines I and II, respectively. Since the 33 ppm val- 
ue corresponds to the average of the principal values of the chemical shift ten- 
sor for trans-tmns methylene sequences of polyethylene or n-paraffins in 
orthorhombic crystalline form [50-52] and the 31 ppm value is close to that 
observed for polyethylene in solution, we simply assign these lines to the 
orthorhombic crystalline and the noncrystalline component, respectively. This 
spectrum expresses a thermal equilibrium state as mentioned above. Hence, it 
reflects strictly the contribution from all components in the sample. However, 
since simple examination of this spectrum does not provide further informa- 
tion, we have examined the relaxation phenomenon of each line. 

I 

1 



II 




, I . , , , I , , 

35 30 

ppm from TMS 



Fig. 4. 50 MHz DD/MAS *^C NMR spectrum of bulk polyethylene with a viscosity-average mol- 
ecular weight of 3.0 X 10® at room temperature. The spectrum was obtained by pulse sequence 
I with a repetition time,Xf = 17,000 s. The chemicalshift is based on that of tetramethylsilane 
(TMS) 
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Spin-Lattice Relaxation. In order to determine whether each resonance line 
comprises a single component, we first measured the spin-lattice relaxation time 
Tic by the pulse sequence developed by Torchia [53] or by the standard satura- 
tion-recovery pulse sequence. The Tic values thus obtained were 2560, 263 and 
1.7 s for resonance line I and 0.37 s for line II. As reported by several investiga- 
tors, the line at 33 ppm is associated with three different Tic values [17, 54, 55]. 
This means that this line is contributed to by three components with different 
molecular mobilities. However, since each component was represented by a sin- 
gle Lorentzian line shape at 33 ppm, they are all assignable to methylene groups 
in the orthorhombic crystalline form or in the trans-trans conformation. The 
component with a Tic of 1.7 s can be assigned to methylene groups with a some- 
what pronounced molecular motion in the vicinity of the crystalline stem end. 
Two longer Tic's, 2560 and 263 s, can be assigned to thicker and thinner lamel- 
lar crystalline components. On the other hand, the noncrystalline component 
which appears at 31 ppm is associated with a single Tic of 0.37 s. This simply 
implies that the noncrystalline component comprises a single phase in as much 
as judged only by Tic. 

In order to determine the content of this noncrystalline line further, we exam- 
ined in more detail the behavior of the spin-lattice relaxation. Figure 5 shows the 
partially relaxed spectra in the course of the inversion recovery pulse sequence 
( 180 °-t- 90 °-FIDdd- 10 s)i 2 o with varying x values. The magnetization that was 
recovered for 10 s in the z direction was turned to negative z direction by 180° 
pulse and the magnetization recovered in z direction for varying x was measured 
in the xy plane under DD. The spectra at different steps of the longitudinal 
relaxation were obtained by Fourier transform and are shown in Fig. 5. In these 
spectra the contribution from the crystalline components with Tic's of 2,560 and 
263 s are eliminated due to the lack of time for recovery at each pulse sequence. 
Therefore, we observed preferentially the relaxation process of the noncrys- 




Fig. 5. Partially relaxed spectra of the bulk polyethylene with My of 3.0 x 10®, taken by inver- 
sion recovery pulse sequence ( 180 °-t- 90 ‘’- 10 s)i 2 o with different relaxation times x’s 




Phase Structure of Polyethylene and Other Crystalline Polymers by Solid-State NMR 



53 



talline component at 31 ppm. See the spectrum at r = 0.4 s, an enlargement of 
which is shown at the right Fig. 5. The line at 33 ppm due to the crystalline com- 
ponent with Tic of 1.7 s is still in the negative z direction, while the line at 31 ppm 
due to the noncrystalline component has already appeared in the positive z direc- 
tion. It is to be noted here that a downfield shoulder is visible in the line at 31 
ppm. This suggests the presence of another line due to the noncrystalline com- 
ponent with a different chemical shift. Although the noncrystalline component 
comprises a single phase judged only by T ic, two noncrystalline components with 
the same Tic and with somewhat different chemical shifts must be assumed. Since 
the two noncrystalline phases assumed here are associated with the same Tic 
value, they should have the same molecular mobility in the Tic relaxation time 
frame but not necessarily the same in the other relaxation time frames because 
of different molecular chain conformation (different chemical shift). 

Spin-Spin Relaxation. In order to determine the content of these two noncrys- 
talline components, we next examined the transverse relaxation by a pulse 
sequence shown schematically in Fig. l-III. The partially recovered magnetiza- 
tion in the z direction for a of 3.5 s was followed by transverse relaxation for 
varying time Tj , and the FID was observed under 'H DD. The spectra at differ- 
ent steps of transverse relaxation were thus obtained, and the result is shown in 
Fig. 6. The spectrum at X( = 0.5 ps (a dead time in this pulse sequence) is assumed 







H 




. . ■ 1 , . . , I . , , , 

35 30 

ppm from IMS 

Fig. 6. Partially relaxed spectra of the bulk polyethylene with My of 3 x 10®, taken by pulse 
sequence II with = 3.5 s.A,B, C,D,E,F, G, and H show the transversally relaxed spectra for 
Tt = 0.5, 20, 60, 100, 140, 500, 2000, and 4000 ps, respectively 
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to contain the contributions from all noncrystalline components and a part of 
the crystalline component with Tic = 0.7 s.As is seen, with increasing T( the con- 
tributions from the crystalline and noncrystalline components with shorter 
transverse relaxation times quickly disappear, leaving only the contribution 
from the noncrystalline component with a longer transverse relaxation time. 
Thus, at = 100 ps only a single peak remains at 31.0 ppm. With further increas- 
ing T(, this single line gradually decreases in intensity but the line shape remains 
essentially unchanged. 

The peak height of the noncrystalline resonance line at 31 ppm during the 
transverse relaxation is plotted against the relaxation time r^ in Fig. 7. It can be 
seen that the overall decay curve (the bottom line in the figure) can be clearly 
resolved into two parts, a rapid decay within 50 ps and a subsequent slow decay. 
The initial slope of the slow decay yielded T 2 C = 2.4 ms, and the T 2 C of the rapid 
decay was estimated to be 44 ps by the usual decay analysis as shown at the top 
of the figure. This indicates that the two noncrystalline components are associ- 
ated with T 2 c’s of 44 ps and 2.4 ms. Together with the longitudinal relaxation 
phenomena shown in Fig. 5, it can be concluded that the 31 ppm resonance com- 
prises two lines, a downfield line with T 2 C of 44 ps and an upheld line with T 2 C 
of 2.4 ms, both associated with the same Tic of 0.37 s. The latter line can be 
assigned to the amorphous component in the rubbery state, since a T 2 C of 2.4 
ms is typical of rubbery polymers. On the other hand, the downfield line can be 
assigned to a noncrystalline component whose molecular motion is severely 
restricted, although its motion in the Tic time frame is the same as the rubbery 
amorphous component. 



t? / ms 

0 0.05 0.1 015 




Fig. 7. Transverse relaxation of 31 ppm resonance line: the peak height at 31 ppm in Fig. 6 is 
plotted against transverse relaxation time Ti 
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Line Shape Analysis of Thermal Equilibrium Spectrum. It has been shown that the 
crystalline region comprises three components with different T i c values and the non- 
crystalline region consists of two regions with different molecular mobilities, i.e. rub- 
bery and rigid noncrystalline components. Based on this fact we analyzed the ther- 
mal equilibrium spectrum shown in Fig. 4. For this purpose, we first determined the 
elementary line shape of each component. Firstly it was found that the elementary 
line shape of the crystalline component can be taken as a Lorentzian function at 33 
ppm. The elementary line shapes of the two noncrystalline components at ca. 3 1 ppm 
were obtained by analyzing the actual spectra during the transverse relaxation shown 
in Fig. 6; the process is shown in Fig. 8. Here A is a longitudinally relaxed spectrum 
for of 3.5 s (spectrum A at T( = 0.5 ps in Fig. 6), which comprises the contribution 
of a crystalline component with Tic of 1.7 s and full contribution fi-om the two non- 
crystalline components. The elementary line shape of the rubbery amorphous com- 
ponent was obtained from the spectra after the contributions from the crystalline com- 
ponent and the rigid noncrystalline component with T 2 C of 44 ps had disappeared. 
As mentioned above, the partially relaxed spectra at if > 100 ps have an essentially 
unchanged line shape at 31 ppm, typical of rubbery polymers. Hence, the partially 
relaxed spectrum at Tj = 140 ps was taken as the elementary line shape of the amor- 
phous noncrystalline component (B in Fig. 8). Spectrum C is a difference spectrum 
that was obtained by subtracting B from A. It comprises the contribution from the 
crystalline component with Tic of 1 .7 s and the rigid noncrystalline component. Since 
spectrum C can be broken down well into two Lorentzian line shapes at 33.0 and 31.3 
ppm, the latter line shape at 31.3 ppm was taken as the elementary line shape of the 
rigid noncrystalline components. Thus all elementary line shapes of the crystalline 
and two noncrystalline components were determined, and were found to be well 
approximated by Lorentzian line shapes with different line widths. 
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Fig. 8. Line shape analysis of the noncrystcJline resonance of the bulk polyethylene with My of 
3 X 10®. A and B correspond to A and E in Fig. 6, respectively, C is the difference spectrum 
(A-B), and D shows the line shape analysis. In D, the composite curve of the component lines 
is mostly superimposed on the experimental spectrum 
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On the basis of the elementary line shapes thus obtained, spectrum A has been 
analyzed in terms of three Lorentzian lines by least-squares fitting; the line width 
and the peak height of each component have been varied while keeping the peak 
positions of 33.0, 31.3, and 31.0 ppm. As shown in Fig. 8D, the composite curve 
of the three Lorentzian functions agrees well with the observed spectrum A. 
Based on this analysis, the total thermal equilibrium spectrum has been analyzed 
into three Lorentzian functions centered at 33.0, 31.3, and 31.0 ppm with use of 
the elementary line shapes thus obtained. Although the crystalline line at 33.0 
ppm comprises three components with different Tic’s, the line is assumed to be 
a single Lorentzian. The result is shown in Fig. 9, where the integrated intensity 
fractions of the respective components are also shown. As can be seen, the com- 
posite curve indicated by a dotted line reproduces well the observed spectrum. 
The integrated intensity fraction of the crystalline line is in good accord with the 
degree of crystallinity estimated by the broad-line 'H NMR analysis; the latter 
also coincides with the value obtained from density measurements. Therefore, 
the result obtained here confirms the conclusion previously obtained by broad- 
line NMR analysis, in terms of chemical shift and relaxation phenomena, that 
the linear polyethylene consists of a crystalline component and two noncrys- 
talline components with different molecular mobilities. 

The half-width of the crystalline component line was estimated to be 18 Hz. 
This value reflects the very stable orthorhombic crystalline phase of this sample. 
The component line shape centered at 3 1 .0 ppm represents the contribution from 
the amorphous phase in which the molecular conformations are changed rapid- 
ly over all permitted conformations. The relatively narrow line width estimated 
as 38 Hz is caused by the rapid molecular motion. The line centered at 31.3 ppm 
represents the noncrystalline phase in which the local molecular motion can 
occur in the same manner as in the amorphous phase (in Tic time frame), but a 
long-range molecular motion accompanying a conformational change related to 
a 10-20 methylene sequence is severely restricted. The wide line width as 85 Hz 



integrated intensity fraction 
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Fig. 9. Line shape analysis of the equilibrum NMR spectrum of bulk polyethylene shown 
in Fig. 4. A\ crystalline component centered at 33.0 ppm, B: crystalline-amorphous interfacial 
component at 31.3 ppm, C: amorphous component at 31.0 ppm. The composite curve of the 
component lines is mostly superimposed on the experimental spectrum 
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indicates that the molecular conformation is mostly fixed over all permitted con- 
formations. Considering the wide symmetric lineshape and the chemical shift, 
it can be concluded that this phase comprises the transition region from the crys- 
talline to the amorphous phase. 

In addition to Tic and T 2 c> the proton spin-lattice relaxation time, Tm, was 
also measured by observing the magnetization under DD that was 
obtained by a crosspolarization from *H magnetization that was recovered in an 
inversion recovery process. The Tih provides not only the proper conditions for 
the CP/MAS measurements of NMRbut also the information about the loca- 
tion of different phases via spin diffusion. As can be seen in Table 1, the two 
noncrystalline lines are associated with well-distinguishable Tm values. That is, 
distinctly different Tm's, 1.61 and 0.39 s, can be recognized. The former value 
indicates a restricted molecular motion in the Tih time frame. Table 1 also con- 
tains the data of another sample with molecular weight of 248,000. For this sam- 
ple, the Tih of the crystalline phase was measured to be 2.20 s. We note that this 
value is very close to 2.04 s for the line at 31.3 ppm, whereas it is significantly 
longer than 0.50 s for the line at 31.0 ppm. This shows that the noncrystalline 



Table 1. Linear polyethylene samples and their characterization 



Sample 








Crystalline 

components 


Noncrystalline 

components 


Molecular 


Crystallinity 




Monoclinic 


Orthor- 


Interfacial 


Rubbery 


weight 










hombic 






Mv 


(l-^)d' 




Parameter 


(34.4 ppm) (33.0 ppm) 


(31.3 ppm) (31.0 ppm) 


Bulk 
















crystals 
















3 X 1Q6 




0.620 


mass fraction 


= 0 


0.66 


0.18 


0.16 








half width (Hz) 




18 


85 


38 








Tih 




1.87 


1.61 


0.39 








Tic(s) 




2560,263,1.7 


0.37 


0.37 








T 2 c (ms) 






0.044 


2.4 


248000 


0.761 


0.755 


mass fraction 


0.06 


0.70 


0.16 


0.08 








half width (Hz) 


50 


18 


86 


37 








TmW 




2.20 


2.04 


0.50 








TicW 




2750,111,1.3 


0.41 


0.41 


Solution 
















crystal 
















91000 


0.834 


0.780 


mass fraction 


0.09 


0.68 


0.23 


~0 








half width (Hz) 


60 


20 


133 










Tm(s) 




1.90 


1.90 










TicW 




220,21,2.0 


0.46 










T 2 c (ms) 






0.04 





^ Crystallinity from density measurements. Crystallinity from broad-line *H NMR. 
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component, which comprises the crystalline-amorphous interphase, and whose 
contribution to the spectrum appears at 31.3 ppm, is located adjacent to the crys- 
talline component so that spin diffusion occurs between these two phases. In 
Table 1 the data of the sample crystallized from dilute solution are also listed. We 
will discuss these data in the following section. 

Molecular Weight Dependence of Phase Structure. Similar line shape analysis 
was performed for samples with molecular weight over a very wide range that 
had been crystallized from the melt. In some samples, an additional crystalline 
line appears at 34.4 ppm which can be assigned to trans-trans methylene 
sequences in a monoclinic crystal form. Therefore the spectrum was analyzed in 
terms of four Lorentzian functions with different peak positions and line widths; 
i.e. for two crystalline and two noncrystalline lines. Reasonable curve fitting was 
also obtained in these cases. The results are plotted by solid circles on the data 
of the broad-line NMR in Fig. 3. The mass fractions of the crystalline, amor- 
phous phases and the crystalline-amorphous interphase are in good accord with 
those of the broad, narrow, and intermediate components from the broad-line 
NMR analysis. 

Provided the crystalline stem length is known and a stacked lamellar struc- 
ture is assumed, the thicknesses of the interphase and amorphous phase can be 
evaluated from these data by the Eqs. (6) and (7): 



^i=^cXi/2Xc (6) 

Ca ~ Cc Xa / X<; (7) 

where and are the thicknesses of the interphase and amorphous phase, and 
Xc , Xa, and x; designate the mass fraction of the crystalline, amorphous, and inter- 
phases, respectively. For the evaluation of and ^a> for samples with different 
molecular weights, we used the values for the crystalline stem length that were 
reported by Voigt-Martin and Mandelkern [56], and by Bassett et al. [57] for the 
linear polyethylene samples crystallized under the same condition. The result is 
shown in Fig. 10, where the thickness of each phase is plotted against molecular 
weight over a very wide range. Samples with very low molecular weights consist 
only of a lamellar crystalline phase and a crystalline-amorphous interphase (in 
this molecular weight range the interphase may be considered as a noncrys- 
talline overlayer, since the amorphous rubbery phase is not present). The thick- 
ness of the interphase gradually increases and levels off as the molecular weight 
exceeds 30,000 with the appearance of the rubbery phase. The thickness of the 
interphase stays unchanged at 34A in the range 30,000-100,000. It increases 
again with enhanced increasing of the thickness of the rubbery phase as the mol- 
ecular weight further increases above 100,000. 

The fact that the thickness of the interphase estimated here stays unchanged 
at 34 A in the molecular weight range of 30,000-100,000, while the mass frac- 
tion and thickness of amorphous phase change remarkably, is particularly 
meaningful. Flory et al. [6, 7] anticipated in 1984 based on their lattice theory 
that the methylene chains that emerge from the basal plane of lamellar crys- 
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Fig. 10. Thickness of the crystalline, amorphous, and crystalline-amorphous interfacial phas- 
es of bulk polyethylene as a function of molecule weight My 



tallites must traverse the crystalline-amorphous interphase as a transition 
region of a thickness of ca. 25A to reach the amorphous phase where methyl- 
ene chains have no correlation with the crystalline phase. This thickness of 34A 
is not so far from the theoretical expectation of 2sA, and since the latter value 
can be changed by parameters assumed in their calculation, the result obtained 
here confirms experimentally their theoretical expectation. As the molecular 
weight increases above 100,000 to 3 x 106, the thickness of the interphase 
increases and reaches 80A. This is because the stacked lamellar structure 
assumed is no longer adequate for samples with such large molecular weights. 
The samples studied here were obtained by the isothermal crystallization at 
130 °C, followed by slow cooling. In samples with very large molecular weights, 
additional crystallization may occur during the cooling and the structure will 
become more complicated so that the simple stacked lamellar crystalline struc- 
ture cannot be assumed. 

3.3 

Discussion 

It is evident that the noncrystalline component is distributed in two phases that 
are associated with the same Tic but different T 2 C values. What does this mean? 
In order to understand this phenomenon we have to refer to the theory of the 
relaxation reviewed in Section 2.2 [43]. Provided the internuclear vector between 
carbon and hydrogen nuclei involves only a single motion, that is, if each term 
of the correlation function of the dipole-dipole interaction between and 
spins evolves exponentially with one correlation time tc (relaxation time of 
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motion) according to Eq. (3), the T 2 C must be same when Tic is the same. This 
fact cannot be understood if these two phases are not assumed to involve more 
than two independent motions dictated by different tc values. 

Consider the relation between the correlation function and its spectral den- 
sity. Slower and faster decays of the correlation function (i.e. slower and faster 
motions) give narrower and wider distributions of the spectral density, respec- 
tively. Figure 1 1 (a) shows some decay curves of the correlation function for 
motions with different tc values and (b) indicates the distributions of their 
respective spectral densities that are obtained by Fourier transform of the decay 
curves. Here A, B, C in (a) are the decay curves with tca. tcb. tcc and A, B, C in 
(b) are the distribution of their respective spectral densities. As can be seen, the 
decay becomes slower and the spectral density distribution becomes narrower 
as the Tc increases. Assume here that tca « ^cb « tcc and the amorphous phase 
involves two independent motions dictated by tca and tcb whereas the crys- 
talline-amorphous interphase involves two motions dictated by tca and tcc- 
Here, tca characterizes a local molecular motion with relation to few carbon 
atoms in the main molecular chain, and tcb>tcc a somewhat long-ranged motion 
with relation to a conformational change of ca. 1 0-20 carbon atoms. In other 
words, it is assumed that somewhat long-ranged motion is different between the 
two phases but local motion is the same, the former is dictated by tcb or tcc and 
the latter by a common relaxation time tca* 



(a) correlation function 




V 



Fourier transform 




Fig. 11. a The free induction decay (correlation function) with different relaxation times tc’s, 
and b Fourier transforms (spectral density) 
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As can be seen from Eq. (1), Tic i® determined by the spectral densities at the 
frequencies of £l>c> and can ± me- Assume here that the distributions of the spec- 
tral densities are as shown in Fig. 1 l-(b), where the frequency range of con, coq, 
and £i>h ± coc is shown by the vertical arrow. Then the Tic's of both of the non- 
crystalline phases are determined by the local motion dictated by tca> indepen- 
dent of the other motions that have no meaningful spectral density at the fre- 
quency range indicated by the arrow. Then the two noncrystalline phases must 
have the same Tic. On the other hand, the T 2 C is determined by the spectral den- 
sity at zero for the solid matter as pointed out at the end of Sect. 2.2. Hence, the 
T 2 c's of the amorphous phase and the crystalline-amorphous interphase are 
respectively determined by the motion of B and C, because the spectral density 
at zero of motion A is negligible. Since the spectral densities at zero frequency of 
both the B and C motions are quite different, the two phases must have quite dif- 
ferent T 2 C values, [58]. 

To this end, we emphasize that the two phases assumed for the noncrystalline 
component are well-defined by the differences in T 2 c- In the crystalline-amor- 
phous interphase, a long-ranged molecular motion accompanying the confor- 
mational change of ca. 10-20 carbon atoms is very slow or almost inhibited. In 
this phase the trans-trans conformation may be somewhat abundant but gener- 
ally all permitted conformations are widely distributed. 



4 

Linear Polyethylene Crystallized from Dilute Solution 
4.1 

Introduction and Approach by Broad-Line NMR 

When linear polyethylene is crystallized from dilute solution, lozenge-shaped 
crystallites are observed under an electron microscope that are sometimes called 
single crystals [59-61]. The thickness of the lozenge-shaped crystallites is gen- 
erally smaller than 10 nm, whereas the width exceeds several pm, and the mole- 
cular chains are oriented approximately perpendicular to the wide lozenge faces. 
The problem here is the molecular chain conformation in the boundary between 
the crystalline lamellae or on the lamellar crystallite surface. Since the crystal 
stem length is much shorter than the extended molecular chain length for sam- 
ples with an ordinary molecular weight, the molecular chains have to participate 
in such a lozenge-shaped crystallite repeatedly by regular or irregular folding. 

Our broad-line NMR analysis showed that this type of sample generally 
consists of the phase structure of lamellar crystallites and noncrystalline over- 
layer with a negligible amount of the noncrystalline amorphous phase [16, 62]. 
In broad-line NMR spectra of solution-grown linear polyethylene samples, a 
narrow component that suggests the existence of a liquid-like amorphous phase 
is hardly recognized. In Table 2, the three-component analysis of the broad-line 
NMR spectra of linear polyethylene samples with different molecular weights 
that were crystallized isothermally from 0.08% toluene solution at 85 °C for 24 
hours under a nitrogen atmosphere is siunmarized. 
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The mass fraction of the narrow component that corresponds to the rubbery 
noncrystalline amorphous phase is as small as 0.003-0.006. The mass fraction 
does not increase appreciably with increasing temperature, but stays almost 
unchanged up to 70 °C. Hence, it is concluded that solution-grown samples do 
not actually comprise a rubbery amorphous phase. This conclusion is confirmed 
by high-resolution solid-state *^C NMR with more detailed information. 

4.2 

Approach by High-Resolution NMR 

Figure 12 shows the DD/MAS *^C equilibrium spectrum of a solution-grown 
polyethylene with a viscosity-average molecular weight of 91,000 at room tem- 
perature. In order to analyze this spectrum, and determine the content of the 
structural components of this sample, the spin-lattice and spin-spin relaxation 
phenomena were examined by similar techniques to those employed for the 
study of the bulk-crystallized samples. Firstly, the spin-lattice relaxation that 
suggests the existence of two different noncrystalline components was not seen. 
The Tic was estimated to be 0.46 s. Hence, the noncrystalline component was 
judged to comprise a uniform phase as far as the Tic value was concerned. To 
examine further the content of this noncrystalline component, the spin-spin 
relaxation was examined. That is, the partially recovered magnetization in the z 
direction for of 3.5s was relaxed in the transverse direction (i.e. in the xy plane) 
for a varying time Xt and the FID was observed under ^H DD. In the spin-spin 
relaxation pulse sequence, all noncrystalline magnetizations quickly disap- 
peared simultaneously within if = 100 ps. The T 2 C was estimated to be ca. 40ps. 
Therefore, the noncrystalline component of the solution-grown sample evident- 
ly forms a Uniphase in as much as judged by both Tic and T 2 c- The Tic of 0.46 s 
and T 2 C of ca. 40ps are almost equivalent to those of the crystalline-amorphous 
interphase of the bulk-crystallized polyethylene (cf. Table 1 ). This suggests that 
the molecular mobility and conformation of the noncrystalline component of 
the solution-grown sample are similar to those in the crystalline-amorphous 
interphase of the bulk-crystallized samples. 

We have examined solution-growth polyethylene samples which differ in mol- 
ecular weight over a wide range. However, it was found that this characteristic 
phase structure does not change appreciably in relation to the molecular weight. 
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Fig. 12. 50 MHz DD/MAS '^C NMR spectrum of solution-grown polyethylene with a viscosity- 
average molecular weight of 91,000 at room temperature. The spectrum was obtained by pulse 
sequence I (shown in Fig. 1) with the repetition time,T^ = 1,500 s 



4.3 

Discussion 

The lozenge-shaped crystallites of solution-grown samples are sometimes 
assumed to be formed by regularly folding molecular chains. The neighboring 
crystalline stems are assumed to be connected tightly by a regularly folded meth- 
ylene sequence such as a ggtgg conformation. However, the NMR results cited 
above do not allow such regular folding molecular chains to form the lozenge- 
shaped crystallite. If the crystal methylene sequence of the trans zig-zag confor- 
mation is connected by a ggtgg sequence, additional resonant lines assignable to 
the gt-tt, gg-tt, tg-gt, gt-gg methylene carbons are expected in the NMR spec- 

trum upheld to the resonance line assignable to the crystalline trans-trans meth- 
ylene sequence, i.e. in the range 30-22 ppm. Here the gt-tt, for example, desig- 
nates the methylene carbon both of whose sides are adjacent to gt and tt meth- 
ylene sequences, respectively. In fact, such resonance lines were distinctly recog- 
nized for crystalline cycloalkanes [c-(CH 2 )„] [63, 64]. I. Ando et al. have also 
reported the appearance of such resonance lines in their CP/MAS spectrum of 
solution-grown polyethylene [65]. We have examined critically the DD/MAS 
NMR spectrum of a polyethylene sample that was crystallized isothermally from 
dilute solution under strictly controlled conditions. However, we could not find 
any trace of such resonance lines. The resonance lines of the regularly folded 
methylene sequence become vague with increasing temperature above room 
temperature because of the broadening of each line even for well-crystallized 
C-(CH 2 ) 60 - 

For this cycloalkane, the Tic’s of the methylene carbons in the folded sequence 
are estimated to be 35-39 s at room temperature and become shorter with 
increasing temperature [66] . Hence, the Tic’s for the regularly folded methylene 
sequence in polyethylene, if it exists, will be appreciably shorter than 35 s at room 
temperature, and much shorter than ca. 220 s of the inner methylene carbons of 
polyethylene (cf. Table 1). Hence, if we measure the DD/CP NMR spectrum 
of the solution-grown polyethylene sample by a single pulse sequence with a rep- 
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etition time of ca. 50 s, the contribution from the folded methylene sequence 
should be recognizable because of the reducing intensity of the inner methylene 
carbons. However, we could not find any trace of such resonance lines. Hence, we 
have to suppose that I. Ando et al. took into account some noise in their CP/MAS 
NMR spectrum, because of insufficient signal/noise ratio, as the resonances 
from the folded methylene sequences, and that their experimental results must 
be revised. Thus, we conclude that the lozenge-shaped single crystal of polyeth- 
ylene is not formed by connecting the crystal sequence with the regularly fold- 
ed methylene sequence. 



5 

High-Pressure Crystallized Polyethylene 
5.1 

Introduction 

It is well-known that when linear polyethylene is crystallized from the melt at 
high pressure and high temperature, a structure of extremely thick crystal lamel- 
lae is produced, sometimes called an extended crystal. In the preceding sections, 
we examined the phase structure of linear polyethylene, crystallized either from 
the melt or dilute solution at atmospheric pressure. However, the phase structure 
of linear polyethylene is thought to differ widely according to its crystallization 
conditions, thermal and mechanical history, if the material is the same. In this 
section, we examine samples crystallized under high pressure with the combined 
use of NMR and electron microscopy [18]. The materials chosen were crystal- 
lized at high pressure, to produce samples of the highest crystallinity which are 
arguably the least complicated of melt-crystallized polyethylene, so enhancing 
the prospects of obtaining definitive conclusions. 



Ti b Ie3 .Chuiict4ZLM.'tioai ofpra.Min-cijo.tJK.tedliitirpoJ.'f ihyltriti-iivij'le- 



Simplt 






Cijo.faUiiidtf’ 


Crjo-til 

thlclliitutiun.) 




Douilf^ 


1HHMR= 1 


Fi4ctiDii2 


2S600 


1.22 


0J)71 


03J2S 


457= 


1400 


Fl4CtiDIi 1 


33000 


1.25 


0S>7Z 


0^44 


580= 


1400 




BfflOO 


Ki.4 


0J)28 


08Se 


S2^ 


1408 


HO20-54P 


231000 


8.2 


0^51 


0878 


12tf 


1428 


HImx IWO 


1.1 X 10^ 


- 


OSflP 


0834 




147.2 



*■ Eitniftittd ly g.pjc . 
b Ej-timAttd fcooKi. dtiiiitf 
■= HME. 

Ptdi thitnvio^aivi, m,*4ivndit4rittof lJ0Kfr(iin.. 

■= H-uivi'btr-i.Ttrigtcria-tilthiliinjj-ti dettuftiii.fd.b'f tiUctEoii UTLicio^^phf . 

^ Hum'b er-iTengtcr|a.tilthil!iieiti detenvuntd'bf iritthod. 








Phase Structure of Polyethylene and Other Crystalline Polymers by Solid-State NMR 



65 



5.2 

Experimental 

Six linear polyethylene samples were examined and the results are collected in 
Table 3. They were melted at 260 °C for 1 h at 4.95 kbar and then cooled slowly. 
Cooling rates over the range of the crystallization temperatures were 0.25 K/min 
for two molecular weight fractions with lower molecular weights (Fraction 1 and 
2) and 2 K/min for the other samples. As can be seen, these high-pressure crys- 
tallized samples are generally characterized by very thick lamellae with high- 
melting temperatures. However, their crystal thickness is still much shorter than 
their extended molecular lengths except for the two samples with lower molec- 
ular weight. Hence, the lamellar structures of these samples are thought to be 
formed by folding molecular chains in essentially the similar manner as samples 
crystallized at atmospheric pressure. 

Broad-line ^H NMR. The line shape analysis of broad-line 'H NMR spectra into 
the broad, medium, and narrow components was carried out for these samples, 
and they were assigned to the crystalline, interfacial, and rubbery amorphous 
phases, respectively. Their mass fractions are shown in Fig. 13 as a function of 
weight-average molecular weight. In Fig. 13, the data of samples crystallized at 
atmospheric pressure shown in Fig. 2 are also shown for comparison by a bro- 
ken line. When crystallized at atmospheric pressure, samples consist only of the 
crystalline phase and a noncrystalline overlayer in the molecular weight range 
lower than 20,000 as pointed out previously. Above this the amorphous phase 
with rubbery molecular motion appears and its mass fraction continues to 
increase with molecular weight; the noncrystalline overlayer also starts to 
increase above a molecular weight of 20,000. On the other hand, for the high- 
pressure crystallized samples, the rubbery amorphous phase does not appear 
below a molecular weight of 1 00,000, and is only just discernible in Rigidex 9 with 
a molecular weight of 1 32,000 (the mass fraction of the amorphous phase can be 




Fig. 13. Mass fractions of three phases of bulk polyethylene versus weight-average molecular 
weight, obtained by broad-line *H NMR analyses. Solid lines indicate high-pressure crystal- 
lized samples, while broken lines show samples which were isothermally crystallized under 
atmospheric pressure 
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read by the upper solid line with data shown by open circles as the distance from 
the upper horizontal coordinate axis). Beyond this molecular weight, the rubbery 
amorphous phase increases concomitantly with decreasing crystallinity but to a 
lesser extent in comparison with that found for the samples crystallized at atmos- 
pheric pressure; even above a molecular weight of 10® its mass fraction hardly 
exceeds 0. 1 0, and the crystallinity remains as high as 0.80. This seems to be char- 
acteristic of the high-pressure crystallized samples. 

High-Resolution NMR Spectrum. Figure 14 shows the DD/MAS NMR 
spectra of pressure- crystallized samples of Fraction 1 , Rigidex 9, HO20-54P, and 
Hifax at room temperature. For all samples, there are distinctly recognizable 
sharp peaks assignable to the orthorhombic crystalline component at ca. 33 ppm 
(Peak I) and the noncrystalline component at ca. 31 ppm (Peak II). For samples 
HO20-54P and Hifax, the resonances assignable to a monoclinic crystalline and 
noncrystalline components are further recognized at ca. 34.4 ppm (Peak HI) and 
31 ppm, respectively. For other samples, the existence of the noncrystalline con- 
tribution is a little vague due to the high crystallinity, but the existence was clear- 
ly confirmed by enhancing the noncrystalline contribution using a shorter r^ in 
the pulse sequence (90 °-FIDdd- t^),, for obtaining NMR spectrum. 

Spin-Lattice and Spin-Spin Relaxation. In order to examine the content of these 
crystalline and noncrystalline components, we examined the spin relaxation of 
each resonance line. Firstly it was found that the line due to the orthorhombic 
crystalline component at 33 ppm involves plural Tic’s for all samples, as summa- 
rized in Table 4. In relation to the Tic for each sample, very long Tic values are 
recognized for three higher molecular weight samples. These values are expected 




. I .... I .... I .... I 
40 35 30 25 

ppm from TMS 

Fig. 14. DD/MAS NMR spectra of different polyethylene samples crystallized under high 
pressure, obtained by a tt/ 4 single pulse sequence [7t/4-FlD-T^]; (a) Fraction 1, (fo) Rlgldex 9, 
(c) HO20-54P, (d) Hifax. The repetition time is indicated in each spectrum 
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from their very thick crystal lamellae as extrapolated values from the values 
reported for samples with thinner lamellar thicknesses [17, 67]. The Tic’s of two 
fractionated samples with lower molecular weights are also longer than 4,000 s. 
However, these Tic’s are still short, if their very thick crystal lamellae are consid- 
ered. The crystal thickness of these samples is equivalent to or longer than the 
extended molecular chain length. Hence the lamellar structure of these samples 
is thought to be formed mostly by fully extended molecular chains without chain 
folding so that some methyl molecular chain ends maybe incorporated inside the 
crystal lamellae. This will cause the relatively shorter Tic against the very thick 
crystal lamellae. For all samples there are recognized plural Tic values in a very 
wide range of several thousands, several hundreds,several tens,and a fewseconds 
as can be seen in Table 4. These plural Tic values seem to correspond to the wide 
distribution of the lamellar thicknesses observed by electron microscope. 

On the other hand, the T i c of the noncrystalline resonance at ca. 3 1 ppm is found 
to be in the range of 0.35-0.50 s for all samples. When judged only by the Tic, the 
noncrystalhne component forms a monophase for all samples. However, the non- 
crystalline line involves distinctly different T 2 C values at least for the three larger 
molecular weight samples. As an example, spin-spin relaxation behavior of the 

noncrystalhne line for Hifax is shown in Fig. 15. The spin-spin relaxation behav- 
ior was examined by the pulse sequence (90°-X( -FIDdd- ^e)n> with r^ = 10s and 
varying T(. The height of the peak of the noncrystalhne component at 31 ppm in 
the course of this spin-spin relaxation is plotted against r^. The overall decay curve 
(the upper line in Fig. 15) can be clearly resolved into two parts, a rapid decay and 
a subsequent slow decay. The slope of the slow decay yielded T 2 C = 2.02 ms, and 
the T 2 C of the rapid decay was estimated to be 0.020 ms by the usual decay analy- 
sis as shown at the bottom of Fig. 15. This indicates that the noncrystalline com- 
ponent comprises two different phases, both associated with the same Tic of 0.35 s 
and distinctly different T 2 c’s of 2.02 ms and 20 ps. The former is assigned to the 
amorphous phase with the rubbery molecular motion and the latter to the crys- 
talline-amorphous interphase. The Tic and T 2 C values for all samples are summa- 
rized in Table 4. 




Fig. 15. >3C spin-spin relaxation behavior of the noncrystalline component of Hifax 
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As can be seen from Table 4, two lower molecular weight samples actually com- 
prise only the crystalline and noncrystalline interlamellar material, devoid of 
amorphous phase. On the other hand, the larger molecular weight samples com- 
prise three phases; the crystalline, amorphous, and crystalline-amorphous inter- 
phase in a similar fashion to the atmospheric-pressure crystallized samples. 
However, we note that the Tze’s of the crystalline-amorphous interphase for two 
higher molecular weight samples is appreciably shorter than those of the atmos- 
pheric-pressure crystallized samples. This demonstrates that the molecular 
chain motion in the crystalline-amorphous interphase of these pressure-crys- 
tallized samples on a Tze time scale is more severely restricted. 

Lineshape Analysis and Phase Structure. Since the pressure-crystallized samples 
generally comprise a crystalline and two noncrystalline phases, we analyzed the 
equilibrium spectra of two higher molecular weight samples shown in Fig. 14, 
and the results are summarized in Table 5. If we assume the typical stacked 
lamellar structure for these samples, the thickness of each phase can be esti- 
mated by using the mass fraction obtained from the line shape analysis and the 
crystalline stem length obtained by other techniques. The thicknesses of the 
crystalline-amorphous interphase and the amorphous phase, and ^a, respec- 
tively, may be estimated by Eqs. (6) and (7) with Xc = Xd - 1 - xa where Xd and Xc 2 
denote the volume fractions of the orthorhombic and monoclinic crystalline 
phases. Since the values of the crystalline stem length or the lamellar thickness 
involve a wide distribution, and also depend on the method of estimation, the sit- 
uation is not so simple. Here, however, the weight-average crystal thickness, esti- 
mated from the nitric acid degradation/g.p.c. technique [68] ,was used. The results 
are shown in Table 5, together with the mass fractions of each phase. 

As can be seen, for both samples the crystalline mass fraction (orthorhombic 
plus monoclinic) agrees well with that estimated from the density measurements 
and broad-line ^H NMR analysis. The thicknesses of the interphases are 3.8 and 
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8.0 nm for HO20-54P and Hifax, respectively. The value obtained for HO20-54P 
is almost the same as the 3.4 nm reported for the atmospheric-pressure crystal- 
lized samples with normal molecular weights in the range 30,000-150,000 (cf. 
Fig. 10). On the other hand, the value for Hifax is significantly larger than the val- 
ues cited in the preceding sections for other types of samples as well as the val- 
ues of 1. 5-2.0 nm which were evaluated theoretically by Flory et al. [7]. An 
increase in the thickness of the interfacial phase was also noted for the atmos- 
pheric-pressure crystallized samples if the molecular weight exceeded 1 50,000, 
as mentioned previously. 

5.3 

Discussion 

As pointed out previously with relation to Fig. 13, the dependence of the phase 
structure on molecular weight shifts to the higher molecular weight side in com- 
parison with samples crystallized at atmospheric pressure. In particular, the rub- 
bery amorphous component appears above a molecular weight of 30,000 for 
atmospheric-pressure crystallized samples,but above a value of 1 00,000 for high- 
pressure crystallized samples. The significance of this appears to lie in the poly- 
dispersity of the samples. The atmospheric-pressure crystallized samples exam- 
ined in the previous section were fractions with a molecular weight polydiper- 
sity {MJMn) of ca. 2, so at M„ = 30,000 the number -averaged extended molec- 
ular chain length is approximately twice the number-averaged crystal thick- 
ness Below the limit of 2 for the ratio X„/^n> a molecular chain participates 
only once in a crystallite and so-called extended crystals are produced. If the ratio 
Xn/^n exceeds 2, a molecular chain participates repeatedly in a crystallite by 
chain folding and the amorphous phase is thought to appear. On the other hand, 
with regard to the pressure- crystallized samples, the rubbery phase could only 
just be detected by the broad-line NMR analysis even for the Rigidex 9 sam- 
ple with Mw = 132,000 (Mn = 13,200). The fit here does not at first sight appear 
so good, since the number-averaged chain length (115 nm) is not much bigger 
than the number-averaged crystal thickness (82 nm). However, during high- 
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pressure crystallization considerable molecular fractionation into discrete crys- 
tal populations occurs, and taking this into account there is good reason to sup- 
pose that the same presence or absence of repeat folding is the determining fac- 
tor. For the samples with higher molecular weights that involve the rubbery 
amorphous component, Xn/^n andXw/^w are ~ 2.2 and 7.5-26, respectively. It is 
then reasonable to assume, even in the pressure-crystallized samples, that a mol- 
ecular chain participates repeatedly in a crystallite or several crystallites by chain 
folding and the amorphous phase appears. 



6 

Ultra High Modulus Polyethylene 
6.1 

Introduction 

Recently, various techniques that produce highly oriented linear polyethylene 
with a ultra high modulus (hereafter, referred to as UHMPE) have been devel- 
oped. In this section, we will examine the structure of the UHMPE that was pre- 
pared by highly drawing a dried gel [69]. Even if bulk polyethylene is uniaxially 
highly drawn by a normal method at a temperature between the Tg and Tm, the 
phase structure is essentially similar to the undrawn sample. That is, it involves 
three phases of the crystalline and two noncrystalline phases, although the mass 
fraction and detailed content of each phase are somewhat different. However, 
UHMPE samples may have a particular phase structure. 

6.2 

Experimental 

The first samples examined were prepared by the method developed by Smith et 
al. [70, 71] and by Matsuo [72]. Sample films of a thickness of ca. 100 pm were 
obtained by drying a gel which was obtained by quenching a 0.4 g/ dl decalin solu- 
tion of linear polyethylene with a molecular weight of 3 x 1 0® from 140 °C in ice- 
water. The samples thus obtained could be drawn to a very high extent because 
of very few intermolecular chain entanglements. However, since they could not 
be drawn highly in one step, they were drawn 10 times at the first step in a sili- 
con oil bath at 145 °C at a rate of 1.6 times/min and then at the second step they 
were drawn so that the final draw ratio was 50, 100, and 150 times. 

DD/MAS NMR Spectra. Figure 16 shows the equilibrium NMR spectra of 
those samples that were obtained by the pulse sequence (7i/4-FIDDD-3Tic)n at 
room temperature. Here the sample A is the undrawn dried gel, B is the sample 
obtained from A by annealing at 145 °C for 4 minutes, and C, D, and E are sam- 
ples drawn 50, 100, and 50 times, respectively. For most spectra there is a recog- 
nized downfield resonance at ca. 33 ppm assignable to the orthorhombic crys- 
talline component and an upfield resonance at ca. 31 ppm assignable to the non- 
crystalline component. 
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Fig. 16. 50 MHz DD/MAS ’^C NMR spectra of highly drawn polyethylene samples with differ- 
ent draw ratios, obtained by a ti/4 single pulse sequence (ri/4-FlD-3T ic) 
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Spin-Lattice and Spin-Spin Relaxations. In order to determine the content of 
these crystalline and noncrystalline resonances, the longitudinal and transverse 
relaxations were examined in detail. It was first confirmed that the noncrystalline 
resonance of all samples is associated with Tic in an order of 0.45-0.57 s. Hence, 
the noncrystalline component of all samples comprises a monophase, in as much 
as judged only by Tic- However, it was found that the noncrystalline component 
of drawn samples generally comprises two phases with different T 2 C values; amor- 
phous and crystalline-amorphous interphases. The dried gel sample does not 
include rubbery amorphous material; it comprises the crystalline and rigid non- 
crystalline components. However, the rubbery amorphous phase with T 2 C of 5.5 
ms appears by annealing at 145 °C for 4 minutes. For the orthorhombic crystalline 
component, three different T ic values, that suggest the distribution of crystallite 
size, were recognized for each sample, as normal for crystalline polymers [ 1 7, 54, 
55] . The Tic and T 2 C of all samples examined are summerized in Table 6. 
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The longest Tic of sample increase with the increases in draw ration. How- 

ever, even for 150 fold drawn sample (E) the Tic does not exceed 1650 s. This 
seems to reflect the phase structure of this type of sample. In addition to the 
orthorhombic crystalline resonance at 32.9 ppm, the existence of downfield res- 
onance at 34.7-34.8 ppm assignable to the monoclinic crystalline component for 
all samples was confirmed by examining the relaxation phenomena. 

6.3 

Phase Structure and Discussion 

The mass fractions of all phases of each sample were obtained by line shape analy- 
sis of the equilibrium spectra, and the results are summarized in Table 7. 

As can be seen, the crystalline fraction (orthorhombic plus monoclinic) 
decreases from 0.70 to 0.54 at a draw ratio of 50 times and increases to 0.82 at the 
largest draw ratio. The amorphous phase that appears at a lower draw ratio 
decreases with increasing draw ratio, accompanying the decrease of the inter- 
phase. Such phase structure as elucidated here will reflect on the various macro- 
scopic properties of samples. 

Figure 17 shows the thermograms of these samples measured with a differ- 
ential scanning calorimeter at a rate of 1 0 °C/min. The temperature at which melt- 
ing is complete is extraordinarily high for most samples with the exception of 
the annealed undrawn sample (sample B) in comparison to the reported values 
for linear polyethylene. The melting temperature of sample A (undrawn dried 
gel sample) is 147.0 °C, which exceeds the equilibrium melting temperature that 
was obtained for this polymer by theoretical extrapolation from the melting tem- 
peratures of n-paraffins by Flory et al. [73] or by the dependence of the observed 
melting temperatures on the crystallization temperature by Mandelkern et al. 
[74]. Hence, the reason for such a high melting temperature could be attributed 
to the very high viscosity at the molten state of this sample due to the very high 
molecular weight. That is, complete randomization of molecular chains does not 
stop when melting is complete while heating at a rate of 1 0 °C/min due to the very 
high viscosity at the melt. Hence, the observed high melting temperature can be 
attributed to the decrease of entropy change of the fusion. 
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Fig. 17. DSC thermograms of highly drawn polyethylene samples with different draw ratios. A: 
undrawn, B\ undrawn, annealed, C; 50 times drawn, D; 100 times drawn, £; 150 times drawn 



As mentioned above, 0.39 fraction of the amorphous phase appears upon 
annealing at 145 °C for 4 minutes for this sample (see Table 7). In compliance 
with this change, the melting temperature of sample B decreases to 138 °C, which 
is slightly higher than the melting temperatures observed for well-crystallized 
samples of this polymer. However, upon drawing, the melting temperature 
increases extraordinarily and three well-distinguishable endothermal peaks are 
observed. For example, for the 50 fold drawn sample (sample C) the melting tem- 
perature increases to over 150 °C and two small additional endothermal peaks 
at 145.0 and 147.9 °C are observed in addition to the main peak at 137.9 °C before 
complete fusion. Three similar endothermal peaks in the fusion are reported for 
a fibrous crystal of ultra-high molecular weight polyethylene [75], a uniaxially 
drawn sample [76] and a drawn sample of single crystal mat [77]. These 
endothermal peaks were attributed from the low temperature side to be the melt- 
ing of randomly oriented crystallites, the solid-solid phase transformation from 
orthorhombic to hexagonal crystal form, and the melting of the hexagonal crys- 
tal [69, 70]. However, since the melting of the hexagonal crystal could not be 
defined for a highly drawn product of a single crystal mat, Furuhata et al. point- 
ed out that the intermediate peak could not be attributed to the solid-solid crys- 
tal transformation [71]. We note here that the intensity of the highest tempera- 
ture peak increases remarkably upon drawing (cf. the isotherm for 100 and 150 
fold drawn samples). Such an intensity increase cannot be understood by the rea- 
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soning pointed out above. Considering the increase in intensity on drawing, the 
origin of the highest temperature peak will be simply attributable to anisotrop- 
ic fusion process in the time scales of the heating rate of 10 °C/min. 

As can be seen from Table 7, although the dried gel does not contain amor- 
phous material, the amorphous phase appears by drawing at 145 °C and could 
not be removed even if drawn as many as 150 times. The amorphous phase is 
thought to be removed to produce excellent ultra-high modulus material. If the 
dried gel is drawn at a slightly lower temperature than 145 °C, the amorphous 
material may be completely removed. We next examined highly oriented poly- 
ethylene samples that were produced by another mode of drawing [78]. A dried 
gel was obtained by casting the solution of linear polyethylene with a molecular 
weight of 6 X 10^ in an aluminum frame cooled by ice water and drying. It was 
then solid-state coextruded at 1 10 °C by the method developed by Kanemoto et 
al. [79] . This solid-state extrusion corresponds to a 5 fold drawing. It was further 
drawn until 215 times at the maximum at 135 °C. The phase structure obtained 
by a similar NMR technique is summarized in Table 8. 

As can be seen, the amorphous phase is not present in all samples and the mass 
fraction of the interphase decreases as the draw ratio increases. By drawing as many 
as 150 times, the mass fraction of the interphase becomes as low as 0.06. A sufficient 
effect of drawing is obtained and the elastic modulus becomes as high as 190 GPa. 



7 

Linear Polyethylene with Randomly Distributed Ethyl Branches 
(Hydrogenated Polybutadiene) 

7.1 

Introduction 

In the previous sections, we have discussed linear polyethylene samples that were 
prepared under different conditions. The polymers, called polyethylene, differ 
widely in molecular structure, i.e. in molecular weight and its distribution, exis- 
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tence of foreign components other than methylene groups and long or short 
branches. According to such differences in molecular structure, the crystallinity 
as well as the phase structure in the solid state differ widely. In this section we 
will examine linear polyethylene that contains a small number of ethyl branch- 
es [21]. 

Examined are two samples that contain 2.2 mol % of ethyl side groups with 
molecular weights of 16,000 and 400,000 (samples P16 and P420, respectively). 
The concentration of ethyl branches was determined by standard high-resolu- 
tion NMR methods [80-82]. These samples were crystallized from the melt 
by quenching at -78 °C. 

7.2 

Phase Structure as Revealed from the Analysis of Resonance Lines 
of Main Chain Methylene Carbons 

Figure 18 shows the thermal equilibrium DD/MAS NMR spectra of samples 
pi 6 and p420. Two sharp resonance lines are observed at 32.4 and 30.5 ppm, that 
can be assigned to methylene carbons in the orthorhombic crystalline and non- 
crystalline phases, respectively. The weak resonance at 1 0.8 ppm can be assigned 
to the methylene carbon of the ethyl branches. Examining the longitudinal and 
transverse relaxation of the methylene resonances at 32.4 and 30.5 ppm, it was 
found that the crystalline resonance can be represented by a Lorentzian function 
though it is associated with three different Tic values and the noncrystalline res- 
onance can be broken down into two components with the same Tic and differ- 
ent T 2 c’s. Based on the relaxation phenomena, a line shape analysis of the equi- 
librium spectra of both samples was performed. Figure 1 9 shows the line shape 
analysis of the equilibrium spectrum of sample P420. 

As can be seen, the resonance of the main chain methylene carbons consists 
of large Lorentzian functions centered at 32.4 and 30.5 ppm and a rather wide 
Lorentzian function at 32 ppm. These functions represent well the orthorhom- 




Fig.18. Thermal equilibrium ’^C NMR spectra for two hydrogenated polybutadiene samples 
P420 and Pi 6, obtained by a r/2 single pulse sequence (ti/4-F1D-t^) with r^= 2000 s 
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Fig. 19. Line shape analysis of the equilibrium spectrum of P420. The large dotted Lorentzians 
centered at 32.4 and 30.5 ppm and rather wide dotted Lorentzian centered at 32 ppm represent 
the orthorhombic crystalline and noncrystcdline Eimorphous phases and crystalline-amor- 
phous interphase, respectively. The dotted curve that is almost completely superimposed on 
the experimental spectrum indicates the composite curve of the component Lorentzians. 
dashed Weakly Lorentzians at 39, 34, 28, and 26 ppm represent the contributions from the 
methine and methylene carbons (a and p to the methine and methylene in the ethyl side 
group), respectively 



bic crystalline and noncrystalline-amorphous phases and crystalline- amor- 
phous interphase, respectively. Weak resonances are also observed at 39, 34, 28, 
and 26 ppm, which can be assigned to the methine and methylene groups (a and 
(3 to the methine and methylene in the ethyl side group). These resonances, how- 
ever, are not pertinent to the present discussion. Thus it is evidenced by the relax- 
ation phenomena and the line shape analysis that the low density linear poly- 
ethylene samples examined comprise three phases similar to linear polyethyl- 
ene; the crystalline, amorphous, and crystalline-amorphous interphase. The 
mass fraction of each phase was obtained by this line-decomposition analysis. 
The thicknesses of the three phases were determined from the mass fractions 
thus obtained with the long spacings determined by the small-angle X-ray dif- 
fraction analyses, assuming a stacked lamellar structure. Raman spectroscopic 
analysis was previously carried out for these samples by Alamo et al. [83]. The 
phase structure of these samples was also obtained by line shape analysis of 
Raman spectra based on the elementary spectra of 1 00% crystalline and amor- 
phous states, and the results are summarized in Table 9. 

The half-widths of 37-39 and 78-88 Hz, respectively, for the crystalline and 
amorphous phases are significantly larger than 1 8 and 38 Hz for those of the bulk- 
crystallized linear polyethylene (cf. Table 1). This is caused by incorporation of 
minor ethyl branches. The molecular alignment in the crystalline phase is slight- 
ly disordered, and the molecular mobility in the amorphous phase will therefore 
be promoted. With broadening of the crystalline and amorphous resonances, the 
resonance of the interphase also widens in comparison to that of bulk-crystal- 
lized linear polyethylene samples. This shows that the molecular conformation 
is more widely distributed from partially ordered trans-rich conformation to 
complete random conformation, characteristic as the transition phase from the 
crystalline to amorphous regions. 
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It is noted that both the mass fractions and the thicknesses obtained by the 
NMR method are in very good agreement with the results reported by Alamo 
et al. using Raman spectroscopy [83]. The thicknesses of the interphases are esti- 
mated to be IsA for two samples, although the molecular weight and the degree 
of crystallinity differ widely. The results obtained here for the model ethylene 
copolymers give further support to the concept that there is a quantitatively 
defined, nonregularly structured interfacial region between the crystalline and 
amorphous phases. 

7.3 

Partitioning of Methyl Carbons Between Different Phases 

In the equilibrium spectra shown in Fig. 18, the resonance of the methyl carbon 
in the ethyl branches (referred to as IB 2 carbon, hereafter) can be seen in ca. 10 
ppm range. It is possible to obtain the information of the partitioning of the eth- 
yl side groups between different phases by analyzing this resonance. In our dis- 
cussion up to now we have tacitly assumed that the ethyl side groups do not enter 
the lattice, i.e. the crystalline phase remains pure. A central question associated 
with this class of copolymers is whether the side group enters the lattice and if so 
under what conditions, i.e. equilibrium or nonequilibrium. A variety of experi- 
mental results, involving many different techniques, have made clear that side 
groups greater than methyl are effectively excluded from the crystal lattice [84]. 
Among these many techniques, solid-state NMR methods have been used to 
study this problem. Although most of the NMR studies agree with the conclusion 
cited above, there are minor disagreements [85-87] . Studies with a hydrogenated 
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polybutadiene that contained 1.7 mol% ethyl branches led to the conclusion that 
a very small proportion, 5-10% of the total branch concentrations, was located 
within the crystalline region. In a study involving polydisperse, broad-composi- 
tion distribution ethylene-butene copolymers, it was concluded that about 9% of 
the total branches were located within the lattice [87]. The concentration of eth- 
yl groups entering the lattice that was deduced from these works is admittedly 
very small. However, since most of these conclusions were deduced without con- 
sideration of the presence of the crystalline-amorphous interphase, it is necessary 
to re-examine the problem by analyzing the resonance of the IB 2 carbon. 

Figure 20 shows the DD/MAS spectrum in the IB 2 resonance range for the sam- 
ple PI 6 that was obtained by a single-pulse sequence with a repetition time of 50 
s. This spectrum is considered to be the equilibrium IB 2 resonance, since the 50 
s is longer enough than 5 times the longitudinal relaxation time of the IB 2 car- 
bon. There is a recognized downfield broad shoulder in addition to the upheld 
main resonance centered at 10.8 ppm [88]. It was found by examining the longi- 
tudinal and transverse relaxations that both the downfield and upheld reso- 
nances are associated with the same Tic of ca. 0.1 s but with different T 2 c’s; ca. 
20 ps for the former, and greater than a few milliseconds for the latter. Based on 
this analysis, the equilibrium IB 2 resonance was broken down into two Lorentz- 
ian functions. The result is shown by the dotted curves in the Fig. 20. As can be 
seen, the composite curve of the two component Lorentzians coincides well with 
the observed spectrum. The intensity of the downfield resonance is very small 
relative to the upheld one. The integrated mass ratio of the downfield resonance 
to the upheld resonance was estimated to be 1:5. Here, the main upheld reso- 
nance can be definitely assigned to the IB 2 methyl carbons in the noncrystalline 
amorphous phase. The question remains as to the proper phase assignment of 
the downfield shoulder. Considering the fact that the T 2 C is as short as 20 ps, this 
resonance couldbe assigned to a branch methyl carbon located either in the crys- 
talline region or in the crystalline-amorphous interphase. The concentration of 
ethyl branches involved is, however, very small, being in the order of only 0.44 
mol%. We note here that the Tic of this component is 0.1-0.2 s, and it is the same 
as in the amorphous phase. 

In principle, a methyl carbon with such a short Tic couldbe assigned to the crys- 
talline phase if the methyl groups were an integral part of the crystal structure. In 




Fig. 20. Equilibrium spectrum and line shape analysis in the 1 B 2 carbon range for the sample 
Pi 6. The spectrum was obtained by a single pulse sequence with a repetition time of 50 s 
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polypropylene, for example, where this condition is satisfied, the Ti c for the methyl 
carbon in the crystal is much shorter than that for the crystalline methine and 
methylene ones, as will be described later [24]. The reason for such a short Tic is 
that there is a threefold rotational motion of the methyl group in a frequency range 
that has an influence on the spin-lattice relaxation time. However, such an influ- 
ence would not be expected if the methyl group were not an inherent part of the 
crystal structure, as is the case of the copolymers being studied here. Hence, based 
on the very short Tic values, the downfield resonance could be assigned to methyl 
carbon in the crystalHne-amorphous interphase. This assignment is consistent 
with other studies involving the location of the ethyl side groups [82]. 

Since a very large proportion of the methyl groups are located in the amor- 
phous region, it is not possible, or reasonable, to assign the downfield resonance 
completely to the crystalline region. A chain passing between the amorphous and 
crystalline regions has to traverse the interfacial region. Hence, some of the side 
groups would have to be located in this region. VanderHart et al. [85, 86] also 
studied similar copolymers by solid-state NMR. The basic spectra that they 
reported are very similar to those presented here. In their analysis of a hydro- 
genated polybutadiene sample containing 1.7 mol% ethyl branches (very simi- 
lar to the one studied here), they found the ratio of the IB 2 carbons of the down- 
field to upheld resonance to be 1 : 1 0. Considering the weak intensity of the down- 
field resonance, this result, and that of the present work, can be considered to be 
very similar. However, they assigned the downfield IB 2 methyl carbon solely to 
the crystalline phase because the T^iOf 8.2 ms for the IB 2 downfield line almost 
agreed with the T/^ of 8.7 ms for the crystalline methylene line. Since only two 
regions, the crystalline and amorphous phases, were considered in their work, it 
was logical to assign the downfield resonance to the crystalline phase. As was dis- 
cussed above, with the establishment of a crystalline-amorphous interphase and 
consideration of the T ic and T 2 c> the downfield resonance can be assigned to the 
interphase. The very close Tip values of the crystalline methylene resonance and 
the downfield line of the 2Bi resonance that were reported can also be under- 
stood logically as the downfield line being assigned to the crystalline- amor- 
phous interphase. Thus it is concluded that the ethyl branches are located in the 
crystalline-amorphous and amorphous phases in a ratio of 1 :5. Very minor eth- 
yl branches may be located in the crystalline region; however, it is not possible 
to estimate such minor concentrations using this NMR technique. 

8 

Poly(tetramethylene oxide) 

8.1 

Introduction 

Until now we have been concerned mostly with crystalline polyethylene. In this 
section we consider the solid-state structure of poly(tetramethylene oxide) [22] . 
Since the melting temperature of this polymer is 42 °C, we examined the struc- 
ture at temperatures below room temperature. The sample was prepared by ring- 
opening polymerization of tetrahydrofuran by using triethyloxonium hexa- 
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fluoroantimonate in methylene chloride at 0 °C. The number- average molecular 
weight determined by osmometry at 28 °C in toluene, was 67,000. The sample was 
melted at 65 °C and then crystallized at 20 °C for one month. The heat of fusion 
and the melting point were estimated to be 8.13 KJ/mol and 42 °C (as the peak 
temperature) by DSC measurements, respectively. It was confirmed by X-ray dif- 
fraction analysis that the crystalline molecular chain of this polymer is in a pla- 
nar zig-zag conformation [89, 90.] 

8.2 

Experimental 

CP/MAS and DD/MAS NMR Spectra. Figure 21 shows the CP/MAS and 
DD/MAS NMR spectra at 0 °C. This polymer has two different kinds of carbon 
atom in the monomeric unit, i.e. a- and (1-methylene carbons to the oxygen atom. 
In accordance with this molecular structure, two resonance Hnes at ca. 72 and 28 
ppm are observed in the CP/MAS spectrum (a). These downfield and upheld reso- 
nances are assignable to the a and (1-methylene carbon atoms, respectively. Here, 
both resonance lines seem to have upheld shoulders. In order to examine these 
upheld shoulders, we measured a DD/MAS spectrum by a single pulse sequence 
without CP. The spectrum (b) shows the spectrum obtained by a single pulse 
sequence with a repetition time of 5 s. Since 5 s is much shorter than the Tic of the 
crystalline component, we could see preferentially the noncrystalline contribution 
that was not clearly recognized in the CP/DD spectrum due to the lesser CP effi- 
ciency. In this spectrum, the resonance lines that were recognized in spectrum (a) 
as only upheld shoulders to the resonance lines of a- and (3-methylene carbons are 
distinctly recognized. Since their chemical shifts are very close to those of the 
respective carbons of this polymer in CDCI3 solution that are shown in (c) as a stick 
type spectrum, they are evidently assigned to the noncrystalline component. 
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Fig. 21. (a) CP/MAS *^C NMR spectrum of poly(tetramethylene oxide) at 0 °C. (fo) The partial- 
ly relaxed spectrum, (c) Stick-type spectrum in CDCI 3 at room temperature 
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Spin-Lattice and Spin-Spin Relaxation Times. In an attempt to examine the 
detail of these noncrystalline resonance lines. Tic and T 2 C were estimated using 
a similar method to that used in the previous sections. The results are summa- 
rized in Table 10. 

The resonance lines at 72.9 and 28.3 ppm are assigned to the crystalline com- 
ponents of a- and (3-methylene carbons because of their longer Tic values. These 
crystalline resonance lines are associated with two Tic values of ca. 209 and 9-10 
s. This shows that both methylene carbons possess two components with differ- 
ent Tic ’s, but this will not be discussed further, since the existence of plural Tic ’s 
is a normal finding for crystalline polymers as discussed in previous sections. 
On the other hand, the resonance lines at 70.9 and 27.0 ppm recognized for a- 
and P -methylene carbons are assignable to the noncrystalline component, 
because these chemical shifts are very close to those in the solution. These lines 
are associated with only one Tic of 0-15 or 0.14 s and two T 2 C values of 7.95 s and 
0.099 ms, or 8.22 s and 0.099 ms, respectively for the a- and p-methylene carbons. 
This suggests that the noncrystalline component involves two components,both 
associated with a same Tic and different T 2 C ’s. The noncrystalline component 
with a T 2 C of 7.95 or 8.22 ms is thought to form an amorphous phase and that 
with a T 2 C of 0.099 ms comprises a crystalline-amorphous interphase. In order 
to confirm this, we examined the elementary line shapes of each component and 
performed the line shape decomposition analysis of the equilibrium spectrum. 

8.3 

Phase Structure and Discussion 

Figure 22-(a) shows the DD/MAS spectrum in the resonance range of a-methyl- 
ene carbon at 0 °C. This spectrum represents the thermal equilibrium state of 
this sample, because it was obtained by a single pulse sequence with the repeti- 
tion time of 600 s longer than 5 times the longest Tic in the system. The spec- 
trum (b) is that of the crystalline component, which was obtained with use of 
Torchia’s pulse sequence [53]. In the equilibrium spectrum, the noncrystalline 
contribution (amorphous plus interfacial) can be seen upheld to the crystalline 
component. Figure 23 shows the elementary line shapes of the amorphous and 
crystalline-amorphous interphases that comprise the noncrystalline resonance. 
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Fig. 22. (a) Equilibrium spectrum of poly(tetramethylene oxide) at 0 °C in the range of the a- 
methylene carbon and (fo) its crystcdline component 





Fig. 23. (a) Partially relaxed spectrum of the a-methylene carbon at 0°C, obtained as an/2 sin- 
gle pulse sequence with = 0.8 s, and (b) transversally relaxed spectrum for 600 qs and (c) 
difference spectrum obtained by subtracting (b) from (a). Spectrum (a) represent all non- 
crystaUine contribution and (b) and (c) represent the contributions of the amorphous and crys- 
talline-amorphous interphases, respectively 



Spectrum (a) shows the DD/MAS NMR spectrum of the a-methylene car- 
bon that was obtained by a single pulse sequence with a repetition time of 0.8 s. 
This is 0.8 s is longer than 5 times the Tic of the noncrystalline component and 
much shorter than Tic of the crystalline component (cf. Table 10). Hence, this 
spectrum represents the contribution from the noncrystalline component that 
consists of amorphous and crystalline-amorphous interphases. Spectrum (b) is 
a partially relaxed spectrum transversely for 600 ps. Since 600 ps is much longer 
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than the of 99 |is for the interfacial component, and shorter than that for the 

amorphous component, this spectrum can be taken as the elementary line shape 
for the noncrystalline amorphous phase. The elementary line shape for the non- 
crystalline interfacial phase was constructed by subtracting spectrum (b) from 
spectrum (a), and is shown as spectrum (c) in the figure. The wide line shape 
located between the crystalline and amorphous resonances shows that the mol- 
ecular conformations in this phase are widely distributed over the permitted con- 
formations stationary in time. 

The line-decomposition analysis of the equilibrium spectrum of the a-meth- 
ylene carbon was carried out using the elementary line shapes thus obtained for 
the three phases. The result is shown in Fig. 24. The composite curve of the 
decomposed components reproduces well the experimentally observed spec- 
trum. The mass fraction of the crystalline component was estimated as 0.60 that 
is described in the figure. Based on the heat of fusion of 8.13 KJ/mol of this sam- 
ple and the value of 14.2 KJ/mol for the crystalline material of this polymer the 
crystalline fraction was estimated to be 0.57. Here the heat of fusion for the crys- 
talline material was obtained from the effect of diluent on the melting tempera- 
ture with use of the relationship developed by Flory [91 ] . The crystalline fraction 
estimated from the NMR analysis is in good accord with the value estimated from 
the heat of fusion, supporting the rationality of the NMR analysis. 

The mass fractions of the amorphous and crystalline-amorphous interphases 
are obtained to be 0.18 and 0.22 as described in Fig. 24. For this sample the long 
spacingwas estimated as 145 A°bysmall-angleX-ray scattering. With use of these 
mass fractions and the long spacing, the thicknesses of the crystalline, amor- 
phous, and crystalline-amorphous interphases were evaluated to be 87, 26, and 
16 A°, respectively, assuming a stacked lamellar crystalline structure. The thick- 
ness of the interphase estimated here as 16 A° is in good accord with the value 
deduced theoretically by Yoon and Flory [7] and the experimental result obtained 
by small-angle neutron-scattering analysis [38, 92]. Nevertheless, the interfacial 
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Fig. 24. Line shape analysis of the equilibrium spectrum of the a-methylene carbon shown in 
Fig. 22-(a). A, B, and C represent the crystaUine, crystalline-amorphous interfacial, and amor- 
phous components, respectively 
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thickness is significantly thinner than that of the bulk-crystallized linear poly- 
ethylenes with molecular weights of 20,000-1 10,000 whose thickness are 34 A° as 
described in Sect. 3. This implies that the transition region from the crystalline to 
the amorphous region of this polymer with the molecular chains including the 
C-O-C bond sequence could be somewhat thinner than that of polyethylene with 
C-C bonds. The molecular chains with the C-O-C bond sequence may be a little 
flexible compared to those with only a C-C bond sequence, if both molecular 
chains comprise planar zig-zag chain conformation. This effect may reflect the 
relatively shorter crystalline Tic value of this polymer in comparison with that of 
polyethylene. Compare the data in Table 10 with those in Table 1. 



9 

Isotactic Polypropylene Crystallize from the Melt 
9.1 

Introduction 

Isotactic polypropylene (iPP), which is highly crystalline polymer,has been stud- 
ied by several groups using solid-state high-resolution '^C NMR spectroscopy. 
Bunn et al. measured CP/MAS *^C NMR spectra of annealed and quenched a- 
form (monoclinic) and (3-form (hexagonal) samples at room temperature [93]. 
Doublets with 1:2 intensity were observed of CH 2 and CH 3 resonance lines of the 
annealed a-form sample. They assumed that these doublets were derived from 
the presence of two different packing sites in crystals on the basis of the arrange- 
ment of the molecular chain in a unit cell of a-form. They also suggested that 
disappearance of these doublets for quenched a- and |3-form samples was due 
to the decrease in crystalline perfection and the absence of the different sites, 
respectively. 

Lyerla et al. measured Tic over a wide temperature range from room temper- 
ature down to 105 K [94], and concluded that Tic ’s of not only CH 3 but also CH 
resonances depend on CH 3 rotational motion, and that the broadening of the CH 3 
resonance below -100 °C is also due to modulation of CH 3 rotational motion at 
the frequency of proton nutation in the presence of the decoupling field. Gomez 
et al. have also reported solid-state high-resolution NMR spectra of isotactic 

polypropylenes [95] . They used samples characterized by X-ray crystallography 
and reconfirmed the results obtained by Bunn et al. 

In this section we will discuss the molecular structure of this polymer based 
on our results mainly from the solid-state NMR, paying particular attention 
to the phase structure [24] . This polymer has somewhat different character when 
compared to the crystalline polymers such as polyethylene and poly(tetrameth- 
ylene) oxide discussed previously. Isotactic polypropylene has a helical molecu- 
lar chain conformation as the most stable conformation and its amorphous com- 
ponent is in a glassy state at room temperature, while the most stable molecular 
chain conformation of the polymers examined in the previous sections is planar 
zig-zag form and their amorphous phase is in the rubbery state at room tem- 
perature. This difference will reflect on their phase structure. 
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9.2 

Experimental 

The iPP sample was used after purification by Soxhlet extraction with toluene. 
The weight-averaged molecular weight of the purified sample was 2.30 x 10^ and 
the isotacticity was mm = 98%. The bulk-crystallized sample was obtained by 
isothermal crystallization at 140 °C for 6 days from 230 °C. The crystal structure 
was confirmed as a-form (monoclinic) by X-ray analysis. 

DD/MAS Spectra at Different Temperatures. Figure 25 shows 50 MHz DD/MAS 
NMR spectra over a wide temperature range together with a stick-type spec- 
trum of the isotactic sequence of this polymer in solution. First of all, for the 
spectrum at room temperature, three resonance lines assignable to CH 3 , CH 2 , and 
CH carbons are visible similarly to those in solution. However, CH 2 and CH res- 
onances of the bulk-crystallized sample exhibit upheld shifts from those in solu- 
tion, reflecting the difference in conformation between the solid and solution 
states. Moreover, CH 2 and CH 3 doublets are clearly observed, which are related 
to the contribution from the crystalline region as reported by Bunn et al. [93]. If 
these doublets stem from the presence of the different packing sites in a unit cell 
as they proposed, such a difference might disappear at higher temperatures and 
the doublets would change to singlet. In fact, with increasing temperature, the 
doublets of CH 2 and CH 3 resonances tend to singlets, as a result of downfield 
shift of the upheld component of each doublet. At 147 °C, that is higher than the 
crystallization temperature of this sample, and a few degrees lower than the melt- 
ing point, the resonances completely changed to singlets (above 87 °C sharp 




Fig. 25. DD/MAS NMR spectrum of isotactic polypropylene at different temperature with 
stick-type spectrum in solution 
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downfield resonance lines appear to the CH 2 and CH carbon resonances, this will 
be discussed later). However, no change in the spectrum at room temperature 
was observed even after heating at 147 °C. Moreover, no appreciable change was 
observed in the X-ray diffraction patterns in the heating and cooling courses in 
this temperature range. A possible explanation for the change of the NMR spec- 
trum with temperature and no change of the X-ray diffractions could be the pack- 
ing effect, which causes the CH 2 and CH 3 doublets as proposed by Bunn et al. 
[93], and by molecular motion in the crystalline regions, which becomes more 
active at higher temperatures and averages each crystalline carbon resonance 
over the two sites. This molecular motion is most likely a threefold jump rota- 
tion about the 3i-helical chain axis. Such a motion would produce little change 
in the X-ray diffraction patterns at higher temperatures as is observed. 

At higher temperatures in Fig. 25, sharp resonance lines become visible down- 
field to the CH 2 and CH resonances as pointed out above. Since their chemical 
shifts are very close to those of the CH 2 and CH resonances of the isotactic 
sequence of this polymer in solution, these resonance lines are assigned to the 
amorphous component. These assignments are also evidenced by examining the 
spin-lattice and spin-spin relaxation phenomena. 

Spin-Lattice and Spin-Spin Relaxation. The Tic’s of CH 2 and CH and CH 3 car- 
bons are summarized in Table 1 1. For the CH 2 and CH carbons, three different 
Tic’s are respectively observed, reflecting different molecular mobilities in the 
crystalline and noncrystalline regions. The longest, shortest, and medium Tic 
components can be assigned to the crystalline, amorphous, and crystalline- 
amorphous interphases, respectively, as examined below in detail. However, only 
two Tic values are recognized for the CH 3 carbon and the chemical shift is very 
close to that of the CH 3 carbon in the isotactic sequence in solution. The reason 
for this is that there is a threefold rotational motion of the CH 3 group through- 
out the crystalline and amorphous regions as in solution. The longer and short- 
er Tic’s (1.31 and 0.47 s) can undoubtedly be assigned to the crystalline and amor- 
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phous components, respectively. The Tic of 1-31 s is significantly shorter than 
Tic for the other carbons such as CH and CH 2 assigned to the crystalline com- 
ponent. This is because there is a threefold rotational motion of the CH 3 group 
even in the crystalline region as pointed out by Lyerla et al. [94]. In addition to 
the crystalline and amorphous phases, the presence of the crystalline-amor- 
phous interphase was ascertained as described later. What is the Tic for the 
methyl carbon in this phase? It should be 1.31 or 0.47 s; that is, the former if the 
threefold rotational motion occurs as in the crystalline phase and the latter if it 
occurs as in the amorphous phase. It maybe important to decide the degree of 
the threefold motion of methyl group in the crystalline-amorphous interphase. 
However, we could not reach any conclusions only using the Tic results. 

9.3 

Phase Structure and Discussion 

In an attempt to investigate the phase structure of this sample, the line shape 
analysis of the CH 2 resonance line in the DD/MAS spectrum at 87 °C that is 
shown in Fig. 25 was examined. The result is shown in Fig. 26-(a). The elemen- 
tary line shape of the crystalline phase was obtained as the line shape of the 
longest Tic component by Torchia’s pulse sequence [53] . It was a doublet and was 
represented approximately by two down- and upfield Lorentzians with an inten- 
sity ration of 2:1 (Spectrum A shown by dotted line in Fig. 26). Since all methyl- 
ene carbons in the a-crystalline form of this polymer are equivalent in the 
intramolecular helical conformation, the origin of the doublet could be attrib- 
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Fig. 26. Line shape analysis of the CH 2 resonance line of bulk iPP at 87 and 20 °C. A, B, and C 
indicate the crystalline, crystalline-amorphous interfacial, and amorphous components, 
respectively 
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uted to an intermolecular packing effect. In fact, the doublet with the ration of 
2:1 can be well understood as having originated from the existence of two pack- 
ing sites with different magnetic environments, examining the unit cell of the a- 
crystalline form. The elementary line shape of the amorphous phase was 
obtained using the spin-lattice and spin-spin relaxations as a component with 
Tic of 0-62 s (cf. Table 11). This elementary line shape, shown by a dotted line as 
Spectrum C in Fig. 26, could be well approximated by a Lorentzian centered at 
46.43 ppm. The relatively narrow line width and the chemical shift close to that 
of the CH 2 resonance line of this polymer in solution show the rubbery amor- 
phous state of this component. 

The line shape analysis of the resonance of the CH 2 carbon at 87 °C was car- 
ried out with use of the elementary line shapes thus obtained for the crystalline 
and amorphous components. However, well reproducible results as shown by the 
dotted lines in Fig. 26-(a) were obtained only when an additional broad Lorentz- 
ian (spectrum B in the figure) was introduced. The broad component B intro- 
duced here corresponds to the component with Tic of 7.1 s at a chemical shift of 
44.82 ppm shown in Table 1 1. Considering its broad line width and Tic, this com- 
ponent can be assigned to the crystalline-amorphous interphase. However, the 
chemical shift is very close to that of the crystalline phase and the Tic is as long 
as 7.1 s. For most crystalline polymers examined in the previous sections, the Tic 
of the crystalline-amorphous interphase is normally as short as a few tenths of 
a second, that is the same as the amorphous phase. On the contrary, the Tic of 
7.1 s assigned here to the interphase is appreciably longer than the 0.62 s of the 
amorphous phase. Such a character of the interphase is also recognized for the 
CH resonance as can be seen in Table 1 1 . For the CH carbon, the Tic of 4.9 s assign- 
able to the interphase is appreciably longer than the 0.32 s assignable to the 
amorphous phase. Taking into consideration the coincidence of the chemical shift 
with that of the crystalline phase, the reason for this will be that the helical mol- 
ecular chain conformation is still held to some extent even in the crystalline- 
amorphous interphase and the molecular motion in the Tic frequency range is 
somewhat restricted in comparison with the amorphous phase. 

Similar line shape analyses for the equilibrium spectra at different tempera- 
tures were performed. At room temperature, where the amorphous phase is in a 
glassy state, the determination of the elementary line shape of the amorphous 
component was a little difficult. However, excellent line-decomposition analysis 
was obtained by introducing a broader Lorentzian centered at the same chemi- 
cal shift as at higher temperatures. The result at room temperature is shown in 
Fig. 26-(b). Here the nature of the component line shapes A and B of the crys- 
talline and crystalline-amorphous interphases is similar to that in the spectrum 
at 87 °C. However, the component line shape for the amorphous component is 
quite different from that at 87 °C that is distributed over a very wide chemical 
shift range centered at the same chemical shift to that at higher temperatures. 
This reflects the glassy state of the amorphous phase. In the glassy state, the mol- 
ecular conformation in the amorphous phase will be distributed over all per- 
mitted conformations stationary in time and randomly in space. The wide com- 
ponent line shape of the amorphous component obtained here at room temper- 
ature well represents this molecular nature of the amorphous phase. 
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The mass fraction of each phase that was obtained by the line shape analysis 
of the CH 2 resonance line at different temperatures is summarized in Table 12 
with Tic- It can be seen that the mass fraction of the crystalline phase (degree of 
crystallinity) stays unchanged at 0.57 over the wide temperature range from 
room temperature to 1 10 °C, while the amorphous phase increases and the inter- 
phase decreases with increasing temperature. The Tic of the CH 2 carbon in each 
phase is mostly unchanged over the temperature range examined; 65~70 s for 
the crystalline phase, 0. 18~0.2 1 s for amorphous phase, and 7.0~7.5 s for the inter- 
phase. This shows that the molecular motion of each phase in the Tic time frame 
is almost the same either in the glassy or rubbery state. 

As pointed out above with relation to the data at 87 °C, the Tic of the crys- 
talline-amorphous interphase is appreciably longer than that of the amorphous 
phase, suggesting the retention of the helical molecular chain conformation in 
the interphase. We also note that a Tic of 65~70 s for the crystalline phase is sig- 
nificantly shorter than that for other crystalline polymers such as polyethylene 
and poly-(tetramethylene oxide), whose crystalline structure is comprised of 
planar zig-zagmolecular-chain sequences. In the crystalline region composed of 
helical molecular chains, there may be a minor molecular motion in the Tic 
frame, with no influence on the crystalline molecular alignment that is detected 
by X-ray diffraction analyses. Such a relatively short Tic of the crystalline phase 
may be a character of the crystalline structure that is formed by helical molecu- 
lar chain sequences. 



10 

Syndiotactic Polypropylene Gel 
10.1 

Introduction 

We studied the structure of bulk-crystallized isotactic polypropylene in the pre- 
vious section, in this section we will examine the gel structure of syndiotactic 
polypropylene (sPP) [25]. Many crystalline polymers form a gel from the con- 
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centrated solution. The network entities that hold the gel structure are reported 
as crystallites for crystalline polymers such as polyvinyl alcohol [96, 97] and 
polystyrene [98, 1 16] . The entities that form sPP gel must also be crystallites. The 
crystal structure and molecular conformation of this polymer in the bulk state 
have been widely studied by X-ray and electron diffraction analyses [99-109], 
solid-state high-resolution NMR [110-113] and vibrational spectroscopy 
[114-116]. With relation to the molecular conformation, three types of helical 
structure have been reported. The most stable form consists of a trans-trans- 
gauche-gauche molecular sequence (referred to as a ffgg form), but a metastable 
form with all trans molecular sequences (tttt form) is also produced by cold 
drawing. In addition to these two forms, the formation of a ttggttttttgg sequence 
has recently been reported when the cold drawn film is exposed to a solvent vapor 
such as toluene [108]. To elucidate the molecular conformations of the bulk-crys- 
talline sPP, solid-state high-resolution NMR spectroscopy has been found 
very effective, since the spectrum reflects sensitively the molecular conforma- 
tion. In this section, we will examine the molecular conformation of sPP gel, refer- 
ring to the above-cited molecular conformations in the bulk state. 

10.2 

Experimental and Discussion 

A sPP sample with a z-averaged molecular weight of 82.2 x 10® was used. Its 
racemic triad was estimated to be 95%. The gel sample was obtained by quench- 
ing a 13.6 wt% 0 -dichlorobenzene solution of the sPP in ice water at 150 °C. 

Molecular Conformation of sPP gel. Figure 27 shows the DD/MAS NMR 
spectrum of sPP gel. This spectrum was obtained by a single-pulse sequence 
(7t/2-FIDDD-T^€)n with the repetition time more than 5 times the longitudinal 
relaxation time Tic. Hence, this spectrum reflects the thermal equilibrium state 
of the gel. For comparison, the spectrum of the bulk ttgg crystal of this sample 
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Fig. 27. Equilibrium DD/MAS ’^C NMR spectrum of sPP/o-dichlorobenzene gel (13.6 wt%), 
obtained by a n/2 single pulse sequence with a repetition time of 300 s. In the lower part the 
equilibrium spectrum of the bulk sPP crystal in ttgg form is shown for reference. The arrows 
indicate the resonance of the amorphous component of each carbon 
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is also shown in the lower part of the figure. At first glance the spectrum for the 
gel is significantly different to that of the bulk crystal. Several additional reso- 
nance lines are recognized in the gel spectrum that are not recognized in the ttgg 
crystal. We have to assign these resonance lines to obtain information about the 
molecular conformation and phase structure of sPP in the gel form. All resonance 
lines in this spectrum are assigned to polymer carbons, because the resonance 
lines from the carbons of the solvent, o-dichlorobenzene, are thought to appear 
in the range of ca. 130 ppm. Referring to the spectrum of this polymer in solu- 
tion, the resonance lines at ca. 20 ppm are assigned to the methyl carbon, and 
those in the ranges of 25-30 ppm and 35-55 ppm to the methine and methylene 
carbons, respectively. 

To provide further detailed assignment of these lines, we examined the longi- 
tudinal and transverse relaxation times. The results are summarized in Table 13. 
The CH 3 , CH, and CH 2 resonance lines visible for the gel sample at 19.9,27.4, and 
46.4 ppm, respectively, are associated with Tic’s of 0.3, 0.2, and 0.2 s and T 2 c’s of 
14, 11, and 13 ms, respectively. The Tic and T 2 c’s recognized here are typical of 
those for the rubbery state of polymers and their chemical shifts are very close 
to those recognized in sPP solution. Hence, these resonance lines are assigned to 
the amorphous component of sPP in the gel form. We note that these lines are 
visible in the gel as sharp lines whereas they are only detectable in the spectrum 
of the bulk crystal as shoulders to the crystalline lines. The reason for this is that 
the gel contains a rubbery amorphous component whereas the amorphous com- 
ponent in the bulk crystal is in the glassy state. Since in the rubbery state the mol- 
ecular conformations are rapidly changed over all permitted conformations, the 
resonance line becomes narrow and visible. On the contrary, since the molecu- 
lar conformations in the glassy state are mostly fixed over permitted conforma- 
tions stationary in time and randomly in space, the resonance line widens and 
becomes invisible. 
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On the other hand, other lines are assignable to crystalline carbons, because 
these lines are generally associated with a relatively long Tic and very short T 2 C 
that are identical to those of the ttgg crystal. The CH 3 resonance line at 20.6 ppm 
assigned to the crystalline carbon is associated with an exceptionally short Tic 
of 0.3 s. This is caused by the threefold rotational motion of the methyl group 
that is permitted even in the crystlline region as discussed with relation to the 
Tic of methyl carbon of the bulk-crystallized iPP in the previous section. Here 
some significant difference is noted in the crystalline CH 2 resonance lines 
between the ttgg crystal and the gel. For the ttgg crystal two split lines are rec- 
ognized for the methylene carbon at 39.0 and 47.7 ppm. This split can be the result 
of the Y-gauche effect. It is semiempirically confirmed that the adjacent methine 
carbon separated by three bonds at the gauche position shifts the methylene res- 
onance by ca. 5 ppm upheld [109-111]. Figure 28 shows the helical molecular 
structure of the ttgg crystal of sPP. In the figure, the side view and the view from 
the c axis are shown. The closed and open circles in the main chain indicate the 
methylene and methine carbons, respectively. In the ttgg molecular sequence 



(a) (b) 




Irons gauche 



(c) molecular sequence in ttgg crystal form 




Fig. 28. Molecular structure of the ttgg crystEil form of sPP. a and b show the relationship 
between the methylene and methine carbons separated by three bonds, In trans (a) and gauche 
(b). c shows the side view and the view from the c axis of the molecular sequence in the ttgg 
crystal form. The closed and open circles in the main chain show the methylene and methine 
carbons, respectively. The methylene carbons marked as A and B are the methylene carbon 
with two gauche effects (the gt-gt carbon) and that without gauche effect (the tg-gt carbon), 
respectively. See the text 
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there are two kinds of methylene carbons as can be seen in the figure; one has 
tow gauche effects and the other no gauche effect. We refer to the former as the 
gt-tg methylene carbon and the latter as the tg-gt methylene carbon hereafter. In 
the figure, the methylene carbon indicated as A is the gt-tg methylene carbon and 
that indicated as B is the tg-gt methylene carbon. 

The peak at 47.7 ppm is assigned to the methylene carbon of the fg-gf sequence 
without such gauche methine carbons in the crystal, whereas the upheld 39.0 

ppm line is assigned the gt-tg methylene carbon with two gauche methine car- 
bons in the ttgg form. On the other hand, three distinguishable lines are recog- 
nized at 39.0,47.7 ppm, and 49.0 ppm to the crystalline methylene carbon in the 
gel. The intensity ration of the three lines is 3:2:1. The chemical shifts of the for- 
mer two, 39.0 and 47.7 ppm, are the same as those of the ttgg crystal. They are 
assignable to the gt-tg and tg-gt methylene carbons in the ttgg form as in the bulk 
crystal. The chemical shift of 49.0 ppm is identical to that of the methylene car- 
bon in the metastable tttt form that is observed in cold-drawn sPP. However, this 
line could not be assigned to the tttt crystal as discussed below. As pointed out 
at the beginning of this section, there is one more molecular chain sequence form, 
i.e., the ttggttttttgg form for this sPP polymer. The methylene chemical shift of 
this sequence could split into three, corresponding to the methylene carbons in 
the tt-tt, gt-tg, and gt-tt sequences. However, for the ttggttttttgg crystal, a reso- 
nance line for the methylene carbon should appear at 44 ppm from the gt-tt 
sequence with one •^-gauche effect. Since no resonance is recognized at 44 ppm, 
the triple splitting of the methylene resonance of the gel cannot be attributed to 
the existence of the ttggttttttgg crystal. 

Since two resonance lines at 39.0 and 47.7 ppm that correspond to those 
observed in the ttgg form and a resonance line at 49.0 ppm that corresponds to 
that in the tttt form are recognized in the gel spectrum, a coexistence of these 
two forms in the gel might be supposed. In an attempt to determine the possi- 
bility of the coexistence of the two forms in the gel, we measured the IR spec- 
trum that is sensitive to the molecular conformation. The number of normal 
vibrational modes depends sensitively on the molecular conformation based on 
the selection rule of the symmetry species. Kobayashi et al. confirmed the vibra- 
tional modes assignable to the ttgg conformation in the IR spectrum for the gel 
from a sPP/carbon disulfide system [117]. However, since we used o- 
dichlorobenzene as solvent, we examined whether the gel structure depends on 
the solvent. 

Figure 29 shows the IR spectra of the tttt and ttgg crystals and the gel. Here the 
spectrum of the gel was obtained by subtracting the contribution from o- 
dichlorobenzene. The differences in the number of detected IR active bands 
between the tttt and ffggcrystals reflect the number of monomeric units contained 
in the fiber period as well as the factor group symmetry. It can be seen that all 
vibrational modes of the gel are completely the same as the ttgg crystal. This indi- 
cates that the molecular conformation in the gel takes a ttgg sequence as report- 
ed for sPP/carbon disulfide gel [116]. The absorption bands characteristic of the 
tttt crystal, such as those at 1131 and 830 cm-1, cannot be recognized in the gel 
spectrum. This distinctly evidences that the crystal molecular conformation in 
the gel takes the ttgg form without coexistence of the tttt forms. This conclusion 
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Fig. 29. IR spectra of the ttttand ffgg^bulk crystals and sPP/o-dichlorobenzene gel (13.6 wt%). 
The spectrum of the gel represents only the contribution from the sPP polymer that was 
obtained by subtracting the contribution from the solvent 

is also supported by recent NMR experiments [118]. The triple splitting of the 

resonance line of the methylene carbon in the NMR spectrum was examined 
at high temperatures. It was confirmed that the resonance line at 49.0 ppm stays 
essentially unchanged upon elevating the temperature until melting. Since the tttt 
form is unstable at temperatures above 50 °C, this result also supports the con- 
clusion mentioned above. Therefore, all splitting lines of the methylene carbon in 
the NMR spectrum can be assigned to the sequence. 

Referring to the intensity ratio of 3:2:1 for resonance lines at 39.0,47.7, and 
49.0 ppm, it is concluded that the resonance of the methylene carbon in the crys- 
talline ttgg sequence in gel splits into two lines assignable to the methylene car- 
bons with two y-gauche effects and without y-gauche effect in the ratio of 1:1, 
the former at 39.0 ppm and the latter further splits into 47.7 and 49.0 ppm. The 
intensity ratio of 3:2:1 for lines at 39.0,47.7 and 49.0 ppm is thus plausibly under- 
stood. The splitting of the resonance line assignable to the methylene carbon 
without a y-gauche effect into 47.4 and 49.0 ppm cannot be explained by a dif- 
ference in the molecular sequence conformation. A possible explanation for 
splitting is given by taking a molecular packing effect into account. The reso- 
nance assigned to the tg-gt methylene carbon further splits into 47.7 and 49.0 
ppm due to the molecular packing effect while the resonance from the gt-tg car- 
bon stays as a singlet at 39.0 ppm. Sozzani et al. observed a peak at 49.0 (50.2 ppm 
estimated by them) in an as-polymerized sPP sample [111]. They did not assign 
it to the tt-tt conformation. Since the tg-gt methylene is located at the outer posi- 
tion of the helical structure of sPP, an enhanced interaction with neighboring 
helical molecular chains is expected. Therefore, there will be two distinct mag- 
netic environments possible for the tg-gt methylene to yield the splitting. The 
splitting width of 1.3 ppm observed here for the gel is comparable to that in the 
a-crystal of iPP where a splitting of 1.0 ppm has been reported as described in 
the previous section [24]. This splitting is understood as being due to the exis- 
tence of inequivalent sites for methylene carbons due to the molecular alignment 
in the crystal lattice. In the case of the sPP gel, the same effect is predicted. 
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Phase structure. It was confirmed in the previous section that the bulk iPP crys- 
tal consists of three phases; the crystalline, noncrystalline amorphous phase and 
crystalline-amorphous interphase. Hence, it is also assumed that the bulk sPP 
crystal forms a three-phase structure in a similar manner. The question here is 
whether the sPP crystal involves such a phase structure in forming a gel or not? 
In order to study this problem, we have analyzed NMR spectra of the sPP gel. 

The noncrystalline contributions to each resonance of CH2, CH and CH3 carbons 
in the DD/MAS NMR spectrum of the gel can be seen, as indicated by the 
arrows in Fig. 27, where their assignment as noncrystalline resonances was con- 
firmed by the spin-lattice and spin-spin relaxation times as described above with 
relation to the results in Table 13. We carried out the line-decomposition analy- 
sis of the resonance lines of the methine and methyl carbons, since these reso- 
nances are most pertinent for the present purpose because of the simplicity of 
the spectral shape. 

Figure 30 shows the component analysis of the resonances of the methine and 
methyl carbons in the equilibrium DD/MAS NMR spectrmn. Here a Lorentz- 
ian fimction is assumed for each component. The rationality for this assumption 
was confirmed by examining the elementary line shapes for each component using 
the differences in the Tic and T 2 C values in a similar way to that described in pre- 
ceding sections. The narrow Lorentzian components centered at 26.2 and 20.6 
ppm, and 27.4 and 19.9 ppm are assignable to the methine and methyl carbons in 
the crystalline and amorphous phases, respectively, as discussed previously (see 
Table 13). In addition to these components, broad Lorentzian components are rec- 
ognized centered at 26.6 and 21.1 ppm for the methine and methyl carbons. It was 





Fig. 30. Line shape analyses of the resonance lines of methine and methyl carbons in sPP/o- 
dichlorobenzene gel. A, B, and C indicate the crystalline, amorphous, and crystalline-amor- 
phous interfacial components, respectively. (This figure was obtained by revising Fig. 7 in Ref. 
25 whose horizontal chemical shift axis was somewhat shifted) 
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confirmed by minute examination of their longitudinal and transverse relaxations 
that these components are associated with T ic values as short as 0.2-0.3 s and very 
short T 2 C values comparable to those of the crystalline phase. In addition, these 
broad Lorentzians are widely distributed in a close range of the crystalline com- 
ponent of each carbon. Hence, these broad components can be assigned to the 
methine and methyl carbons in the transition phase from the crystalhne to amor- 
phous. In this phase the molecular conformations are widely distributed over the 
permitted conformations. The narrow line width of the amorphous components 
of both the CH and CH 3 carbons indicates that the molecular conformations are 
rapidly interchanged between the permitted conformations in the amorphous 
phase in the sPP crystal to form the gel. We conclude that the sPP crystal of the gel 
involves three phases, the crystalline, amorphous, and an interphase. 

The mass fractions of these three phases are shown in Table 14. the crystalline 
fraction is relatively small as 0.43 or 0.44. This low level of crystallinity may arise 
from relatively strong molecular entanglement due to the network structure. 

Phase Structure as Revealed from the Mobility of the Solvent. The phase struc- 
ture of the sPP crystal in the gel form, which was elucidated by the line-decom- 
position analysis of the DD/MAS NMR spectrum, will reflect on the mobili- 
ty of the solvent in the gel. The mobility of the solvent can be examined by the 
longitudinal relaxation of resonance lines assigned to the carbons of the solvent. 
Figure 31 shows the longitudinal relaxation for the line at 130 ppm of the 0 - 
dichlorobenzene. The open circles indicate the data of the pure solvent and the 
closed ones those of the solvent in the gel. As can be seen, the relaxation of the 
pure solvent evolves exponentially with a Tic of 3.0 s, whereas that of the solvent 
in the gel evolves nonexponentially. This indicates that there are some solvent 
molecules in the gel that differ in their mobility. We assume here that the longi- 
tudinal relaxation of each component of the solvent evolves exponentially. Then 
the longitudinal relaxation of the total solvent follows the relationship: 

f =Z^o, [l-exp(-T/Tii)] 

i 

Here r is the evolution time, / denotes the total longitudinal magnetization of the 
solvent at r and Io;the longitudinal magnetization of i-th component at r = °o. Ti, 
is the spin-lattice relaxation time that dictates the evolution of the magnetiza- 
tion of the i-th component. The Ej fo; is the eqilibrium magnetization of the total 
solvent at r = 00 . Since the equilibrium magnetization of each component is pro- 
portional to the mass fraction, /q; /^i Jo; provides the mass fraction of the i-th 
component. 
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Fig. 31. longitudinal relaxation of the 130 ppm resonance line of o-dichlorobenzene. • indi- 

cates data in gel and O data in the pure solvent 
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Using this relationship the longitudinal relaxation of the solvent shown in Fig- 
ure 31 was analyzed. The longitudinal relaxation could be well fitted as the sum 
of the relaxations of two components, the relaxation of each is dictated with the 
Tic’s of 2.2 and 9.2 s. It is known that in the gel of a polymer/solvent system the 
solvent can be generally classified into two components that differ in mobility, 
referred to as free and bound solvents. The two Tic’s of 2.2 and 9.2 s obtained 
here correspond to those for these free and bound solvents. The shorter Tic of 
2.2 s is almost equal to that of pure o-dichlororbenzene. Hence, the other com- 
ponent with Tic of 9.2 s can be assigned to the bound solvent. The mass fractions 
of these two components of the solvent that are obtained from /oi hi are list- 
ed in Table 15, together with the mass fractions of the amorphous and interphase 
of the sPP polymer in the gel. 

As can be seen from the Table 15, 70 and 30% of the solvent are respectively 
in the bound and free state. The mass fractions of the bound and free solvents 
are in rough accordance with those of the crystalline-amorphous interphase and 
the amorphous phase in the two noncrystalline phases of the polymer. This result 
suggests that the solvent exists in the two noncrystalline phases of the polymer, 
as the bound solvent in the crystalline-amorphous interphase and as the free sol- 
vent in the amorphous phase, leaving the crystalline phase pure. It is concluded 
that the sPP/o-dichlorobenzene gel involves three phases, ( 1 ) the pure crystalline 
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phase of the polymer, (2) the amorphous phase which consists of amorphous 
polymer molecules and free solvent, (3) the crystalline-amorphous interphase 
which consists of polymer molecules and bound solvent. 

In the amorphous phase, the molecular conformation and mobiltiy of poly- 
mer are similar to those in the rubbery state of the pure polymer and the mobil- 
ity of the solvent is almost the same as in the pure state. In the crystalline-amor- 
phous interphase, the molecular conformation and mobility of the polymer are 
severely restricted even if it contains the solvent molecules. Since the solvent mol- 
ecules are trapped among polymer molecular chains in the interphase, the flu- 
idity as a liquid is restricted. The molecular chain entanglement of the polymer 
in the interphase plays an important role in preventing the solvent molecules 
from flowing out of the network structure. A high percentage of the interphase 
region will be essential to hold the solvent in the network. On the other hand, the 
solvent molecules in the amorphous phase of the polymer are mostly free from 
surrounding polymer molecules. They could have a fluidity as pure solvent and 
could be segregated from the gel with ease. If the polymer chain does not exist 
to a significant extent in the interphase region, gelation will not occur since most 
of the solvent molecules are in a free state. Alternatively, the aggregation of poly- 
mer molecular chains will lead to a precipitation such as crystal growth from 
dilute solution. 



11 

Concluding Remarks 

We have reviewed briefly our studies of the solid-state structure of some crys- 
talline polymers. It has been found that bulk-crystallized linear polyethylene 
involves three phases that are distinctly different from each other in molecular 
conformation and mobility; the crystalline, the noncrystalline amorphous phase 
and the crystalline-amorphous interfacial phase. Such a three-phase structure is 
essentially universal for many semicrystalline polymers. In the case of polyeth- 
ylene, the amorphous phase can be absent for the solution-grown sample or ultra 
high modulus sample (high modulus polyethylene). However, even if crystallized 
from the melt at high pressure and high temperature or highly drawn at high 
temperatures under normal conditions, the samples generally involve a similar 
three-phase structure. In addition to polyethylene, we have briefly reviewed the 
study of poly(tetramethylene oxide) and polypropylene that, respectively, have a 
planner zig-zig and helical molecular chain conformation. It was found that these 
polymers involve essentially a similar three-phase structure, characteristic of the 
molecular structure, as well as the thermal and mechanical history of samples. 
Even the sPP crystal that forms a gel network exhibits a phase structure involv- 
ing a crystalline-amorphous interphase. 

The study of the phase structure of semicrystalline polymers, particularly that 
of the interfacial region, is of foremost importance, since it has an intimate rela- 
tionship with the macroscopic properties of polymers. Many techniques have 
been used to study this problem. To determine the mass fraction of each phase 
of the three-phase structure, such techniques as broad-line NMR or Raman 
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spectroscopy may be useful. However, these methods are not applicable general- 
ly to crystalline polymers other than polyethylene and they cannot be used to 
characterize the phase structure in a molecular level. For example, in the Raman 
spectroscopy the elementary line shape of each phase could not be obtained with 
ease. For polyethylene, the elementary line shapes of the crystalline and amor- 
phous components could be determined referring to the line shapes of the 1 00% 
crystalline and molten states [32]. However, it is not possible to obtain such ele- 
mentary line shapes for other polymers. As the result there are some conflicting 
conclusions, particularly about the crystalline-amorphous interphase [119]. 

On the other hand, in the solid-state high resolution NMR, elementary line 

shape of each phase could be plausibly determined using magnetic relaxation phe- 
nomenon generally for crystalline polymers. When the amorphous phase is in a 
glassy state, such as isotactic or syndiotactic polypropylene at room temperature, 
the determination of the elementary line shapes of the amorphous and crys- 
talline-amorphous interphases was not so easy because of the very broad line 
width of both the elementary line shapes. However, the line-decomposition analy- 
sis could plausibly be carried out referring to that at higher temperatures where 
the amorphous phase is in the rubbery state. Thus, the component analysis of the 
spectrum could be performed and the information about each phase structure 
such as the mass fraction, molecular conformation and mobility could be 
obtained for various polymers, whose character differs widely. 
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This review summerizes the recent advances in characterization of some speciai intractabie 
macromoiecuies in soiution by iaser iight scattering. Since both static and dynamic iaser iight 
scattering (LLS) are theoreticaUy weii estabiished, we focus the discussion on experimentai 
detaiis, such as the design of a high-temperature LLS spectrometer, the sampie ciarification, a 
novei differential refractometer, and some newly developed methods for data analysis which 
include a combination of static and dynamic LLS leading to molar mass distribution deter- 
mination and LLS calibration of gel permeation chromatography. 
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List of Symbols and Abbreviations 



Sect. 1 

LLS 

GPC 

FFF 

<M> 

P 

M 

Mn 

Mw 

Mz 

fn(M) 



laser light scattering 
gel permeation chromatography 
field flow fractionation 
average molecular weight 
an integer number 

molecular weight for monodisperse species 
number-average molecular weight 
weight-average molecular weight 
intensity-average (z-average) molecular weight 
differential number distribution 



Sect. 2 



I(t) 

n(t) 

Rw(6) 
dissolution 
dISsolvent 
dIS standard 

<I> 

Rw,standard(6 ) 

n 

^standard 

a 

or dRg> 

q 

K 

Na 

dn/dC 

Xo 

c 

A2 

QELS 



scattering intensity at time t 

scattering photon counts at time t 

excess Rayleigh ratio at angle 0 

time-averaged scattering intensity of solution 

time-averaged scattering intensity of solvent 

time-averaged scattering intensity of standard 

time-averaged scattering intensity 

Rayleigh ratio of standard at angle 0 

refractive index of solution 

refractive index of standard 

detection geometry parameter 

root-mean square z-average radius of gyration 

scattering vector 

optical constant 

Avogadro number 

specific refractive index increment 

wavelength of the laser light in vacuum 

concentration of polymer solution 

second virial coefficient 

quasi-elastic light scattering 
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G(2)(t,q) 

t 

t 

PCS 

A 

P 

I(t,q) 

CONTIN 

G(D 

r 

E(t,q) 

<r> 

D 

kd 

f 

G(D) 

<D> 

G 



intensity-intensity time correlation function at q and 

delay time t 

time variable 

delay time variable 

photon correlation spectroscopy 

normalized first-order electric field time correlation 

function at q and delay time t 

baseline 

coherence parameter 

scattered intensity at time t and q 

Laplace inversion program 

line-width distribution 

line-width for monodisperse species 

scattered electric field at time t and q 

average line-width 

translational diffusion coefficient 

diffusion second virial coefficient 

dimensionless parameter 

translational diffusion coefficient distribution 

average translational diffusion distribution 

proportionality factor 



Sect. 3 

PPTA or Kevlar poly(l,4-phenyleneterephthalamide) 

Teflon and Tefzel fluorocarbon polymers 

PEL-PCL polyfethylene terephthalate-co-caprolactone) 

An refractive index difference between solvent and solution 



Sect. 4 

ko 

OCD 

fw(M) 

t^DLS 

'^'^^w.calcd 

hlw, measd 

ERROR(aD) 

ERROR(kD) 

[R] 



F WjCum' 

SEC 

V 

A 

B 

C(V) 



(M) 



scaling constant in D=kDM “® 
scaling constant in D=kDM"“'^ 
differential weight distribution 

weight-average molecular weight calculated on the basis of 
G(D) from DLS 

weight-average molecular weight measured from static LLS 

error function of an 

error function of kp 

intrinsic viscosity 

scaling constant in [q]=k,^M “t 

scaling constant in [q]=kr,M “t 

cumulative weight distribution 

size-exclusion chromatography 

elution volume 

calibration constant in V=A+B log (M) 
calibration constant in V=A+B log (M) 
elution volume distribution 
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weight -average molecular weight calculated on the basis of 


■^'^^WjCalcd 


C(V) from GPC 


Sect. 5 


hlw.app 


apparent weight-average molecular weight 


V 


refractive index increment for the whole sample 


v(M) 


refractive index increment for copolymer with molecular 
weight M and weight distribution fw(M) 


Wa(M) 


weight fraction of A for a given polymer chain with 
fixed M and F^fM) 


Wb(M) 


weight fraction of B for a given polymer chain with 
fixed M and F^fM) 


Va 


refractive index increment of polymer A 


vb 


refractive index increment of polymer B 


Wa 


weight fraction of A 


Wb 


weight fraction of B 


Vsl 


V in solvent 1 


Vs2 


v in solvent 2 




Va in solvent 1 


Vb.sI 


Vb in solvent 1 


Va,s2 


Va in solvent 2 


Vb,s2 


Vb in solvent 2 


fw,app(M),sl 


apparent weight distribution in solvent 1 


f 

%,app(M),s2 


apparent weight distribution in solvent 2 


Mw,L 


low weight-average molecular weight polymer 


Mw,H 


high weight-average molecular weight polymer 


Xl 


weight fraction of low molecular weight polymer 


Xh 


weight fraction of high molecular weight polymer 


GlIO 


line-width distribution for low molecular weight polymer 


GH(r) 


line-width distribution for high molecular weight polymer 


Al 


area covered by Gl(F) 


Ah 


area covered by Gh(F) 


Ar 


area ratio Al/Ah 


R 


radius of a colloidal particle 


9 


density of particle 


(Mw)cal 


calculated weight-average molecular weight 


b 


thickness of the hydrodynamic layer 



1 

Introduction 

Differing from small molecules, typical polymers have an average molar mass of 
~10^ g/mol or higher and a wide distribution of chain length. For a given type of 
polymer, its properties, even its appearance, are greatly influenced not only by 
its average molecular weight, but also by its molecular weight distribution. 
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Therefore, the development and application of a polymer often require a precise 
characterization of these quantities. 

A number of methods including laser light scattering (LLS) are routinely used 
to determine the average molecular weights and molecular weight distribution 
of a polymer. They include end-group chemical analysis, vapour pressure 
osmometry, membrane osmometry, sedimentation equilibrium, static (classic) 
LLS and very recently developed matrix-assisted time-fly mass spectroscopy as 
absolute methods in the sense that these do not require a calibration with a set 
of narrowly distributed polymer samples with known molecular weights. The rel- 
ative methods include viscometry, gel permeation chromatography (GPC), field 
flow fractionation (FFF) and dynamic LLS. 

The average molecular weight <M> of a polydisperse polymer sample is gen- 
erally defined by 

<M> = 1 “ fn(M)MP dM / f„{M)MP"^ dM (1.1) 



if p is an integer, where fn(M) is the number distribution of molecular weight M. 
Thus, P=1 for the number-average molecular weight (Mn), P=2 for the weight- 
average molecular weight (M^) and p=3 for z-average molecular weight (M^). 
For example, Mn is measured by end-group analysis and osmometry; Mw.by GPC, 
sedimentation equilibrium and static LLS; and Mz, by sedimentation equilibri- 
um and dynamic LLS. In practice, the ratio of Mw/Mn is called the polydispersi- 
ty index and is conveniently used to characterize the distribution width of a giv- 
en polymer sample. 

Static light scattering as a classic and absolute method has been long and 
widely used to characterize both synthetic and natural macromolecules [1-4]. 
In the last two decades, thanks to the advances in stable laser, ultra -fast electronics 
and personal computers, laser light scattering (LLS), especially dynamic LLS 
(denoted here by DLS), has gradually changed from a very special tool for phys- 
ical chemists to a routine analytical tool in polymer laboratories or even to a daily 
quality-control device in production lines [5-8]. The LLS instruments commer- 
cially available nowadays are normally capable of making both static and 
dynamic measurements simultaneously. 



2 

Basic Principles of Laser Light Scattering 

When a monochromatic, coherent light is incident into a dilute macromolecule 
solution, if solvent molecules and macromolecules have different refractive 
index, the incident light is scattered by each illuminated macromolecule to all 
directions [9, 10]. The scattered light waves from different macromolecules 
mutually interfere, or combine, at a distant, fast photomultiplier tube detector 
and produce a net scattering intensity I(t) or photon counts n(t) which is not uni- 
form on the detection plane. If all macromolecules are stationary, the scattered 
light intensity at each direction would be a constant, i.e. independent of time. 
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However, in reality, all macromolecules in solution are undergoing constant 
Brownian motion, and this fact leads to fluctuation in I(t). The fluctuation rate 
can be related to the translational diffusion of the macromolecules. The faster 
the diffusion, the faster the fluctuation will be. 

2.1 

Static Laser Light Scattering 

In static LLS, the angular dependence of the excess absolute time-averaged scat- 
tered intensity, known as the excess Rayleigh ratio, Rw(6)> is normally measured, 
where Rw(6)~(<I^solution~<l^solvent)/<I^standard " Rw,standard(6) ’ (n/Hstandard)* 
with <I> and n denoting the time-averaged scattering intensity and the refrac- 
tive index, respectively, and 1 < a < 2 depending on the detection geometry. Thus, 
if the scattering volume is selected by a slit, a=l, and if the scattering volume is 
selected by a pinhole much smaller than the beam diameter, a=2. If the solution 
is very dilute, Rw(6) at a relatively small scattering angle 6 can be related to the 
weight-average molecular weight Mw, the second virial coefficient A 2 and the 
root-mean square z-average radius of gyration (or written as <Rg> for 

simplicity) by the expression [9, 10] 

~ (1 -I- -<Rg>q^) -H 2 A 2 C -H... (2.1) 

Mw ^ 

where C is the mass concentration of the polymer component, K=4Tr^n2 (dn/dC)^ 
/(NaAo ) and q=(4Tm/Xo) sin(0/2), with Na, dn/dC, n and Xq being the Avogadro 
constant, the specific refractive index increment, the solvent refractive index, and 
the wavelength of the laser light in vacuum, respectively. With Rw(6) measured 
at a series of C and 0, we can determine Mw, Rg, and A 2 by use of a Zimm plot, 
which allows both 0 and C extrapolations to be made on a single grid [11]. Fig- 
ure 1 shows a Zimm plot for an alternating copolymer of ethylene and tetraflu- 





sin^(e/2) 1000 C 

Fig. 1. Typical Zimm plot for an alternating copolymer of ethylene and tetrafluoroethylene (Mw 
= 5.4 X 10^ g/mol, Rg= 45.4 nm andA 2 = 1.97 x 10 mol mL/g^) in diisobutyl adipate at 240 °C 
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oroethylene (Mw=5.4xl0^ g/mol, Rg=45.4 nm and A 2 =l. 97x10“^ mol mL/g^) in 
diisobutyl adipate at 240 °C [ 12] . It should be noted that Eq. (2.1 ) is valid under 
the restriction that the polymer solution exhibits no adsorption, no fluorescence, 
and no depolarized scattering. As for rigid and nearly rigid rods causing depo- 
larized scattering, readers can refer to the excellent review article of Russo and 
the references therein [ 13] . As for the correction of adsorption and fluorescence, 
readers are advised to refer to the characterization of Kevel in concentrated sul- 
furic acid by Qing et al. [14-16]. 

2.2 

Dynamic Laser Light Scattering 

DLS measures the intensity fluctuation instead of the average light intensity (this 
is where the word 'dynamic' comes from), and its essence may be explained as 
follows. When the incident light is scattered by a moving macromolecule or par- 
ticle, the detected frequency of the scattered light will be slightly higher or low- 
er than that of the original incident light owing to the Doppler effect, depending 
on whether the scatterer moves towards or away from the detector. Thus, the fre- 
quency distribution of the scattered light is slightly broader than that of the inci- 
dent light. This explains why dynamic light scattering is sometimes called qua- 
si-elastic light scattering (QELS). The frequency broadening (~1 0^-10^ Hz) is so 
small in comparison with the incident light frequency (~10^^ Hz) that, if not 
impossible, it is very difficult to detect it. However, it can be effectively recorded 
in the time domain through an intensity-intensity time correlation function 
G^^)(t,q) in the self-beating mode. Thus, dynamic light scattering is also known 
as photon correlation spectroscopy (PCS). 

G^^^(t,q) can be related to the normalized first-order electric field time corre- 
lation function jgA* (t,q)| by [9, 10] 



g(^^ (t,q) = < I(t,q)I(0,q) > = A[1 +|3 



, 0 ) 



(t,q) 



( 2 . 2 ) 



where A is the measured base line, (3 is a parameter depending on the coherence 
of the detection optics, and t is the delay time. 

Eor a monodisperse sample, gA^Eq) is theoretically represented by 



g^*^(t,q) 



G 



-rt 

e 



(2.3) 



where G and T are the proportionality factor and the line-width respectively. Eor 
dilute solutions, E measured at a finite scattering angle is related to C and q by 
[17] 



r = q2D(l + kdC)(l + f<R^>,q2) 



( 2 . 4 ) 
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Here D is the translational diffusion coefficient of the solute molecule at C ^ 0 
with C the mass concentration of the solute, kd the diffusion second virial coef- 
ficient, f a dimensionless parameter depending on polymer chain structure and 
solvent, and <Rg^> the mean square radius of gyration of the polymer chain. 
Hence, for C and q small enough, Eq. (2.3) may be approximated by 

g^*^(t,q) = (2.5) 



The proportionality factor G may depend on many characteristics of the solute 
polymer, but since D is the only polymer-dependent variable in Eq. (2.5 ), it is rea- 
sonable here to treat G as a function of D only for a homologues series of poly- 
mers. Thus, for a polydisperse polymer sample with a continuous distribution of 
molecular weight M, Eq. (2.5) maybe generalized as 

g^^Huq) = rGlDle'l'^MD (2.6) 

•'0 



where G(D) is called the translational diffusion coefficient distribution. This 
equation is the basic of the entire discussion in the present paper. Note that since 
g^^Hhq) unity as t goes to zero, we have 

rG(D)dD = l (2-7) 

•'0 



that is, G(D) is normalized. Equation (2.6) indicates that once g(^^(t,q) is deter- 
mined from G^^^(t,q) through Eq. (2.2), G(D) can be computed by Laplace inver- 
sion [ 1 8-24] . Among many others, the program called CONTIN [24] is one of the 
most widely used and accepted for this computation. However, due to the band- 
width limitation of photon correlation instruments, some unavoidable noises, 
and a limited number of data points, the data obtained for g^^^(t,q) do not always 
provide information necessary and sufficient to determine G(T) uniquely. This 
difficulty is well-known as the ill-posed problem. Thus, in practice, reducing the 
noises in the measured G^^^Lq) becomes more important than choosing a pro- 
gram for data analysis. 

Figure 2 illustrates the g^^Hhq) data for chitosan (Mw=l. 06x10^ g/mol and 
<D>=5.92x 10 * cm^/s) in aqueous 0.2 M CH 3 COOH/O.I M CH 3 COONa at 25 °C, 
0=45° and C=4.96xl0“^g/ml. Here,<D> is the average diffusion coefficient giv- 
en by 

<•00 

<D>=[ G(D)DdD (2.8) 

•^0 



Figure 3 shows the G(D) function for the same system as in Fig. 2 at 0 ^ 0 and 
C^O. 
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t /(ms) 

Fig. 2. Typical normalized intensity- intensity time correlation function for chitosan (M„ = 1 .06 
X 10^ g/mol and <T> = 2.19 ms) in 0.2 M CH 3 COOH / 0.1 M CH 3 COONa aqueous solution at 
T = 25 °C, e = 45° and C = 4.96 X lO '* g/mol 




Fig. 3. Typical translational diffusion coefficient G(D) for chitosan (M«= 1.06 x 10^ g/mol and 
<D> = 5.92 X 10 ® cm^/s) in 0.2 M CHsCOOFl / 0.1 M CHsCOONa aqueous solution at T = 25 
°C, 0 — ^ 0 and C — ^ 0 



3 

Experimental Section 
3.1 

Solvent Selection 

If a macromolecule can be dissolved in more than one solvent, the choice of the 
solvent for laser light scattering measurement should be made generally accord- 
ing to the following three criteria: 1) it is colorless so that the adsorption cor- 
rection can be avoided, 2) it has a higher contrast, i.e. a higher value of the spe- 
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cific refractive index increment dn/dC, so that the signal-to-noise ratio is 
increased, and 3) it is less polar and less viscous so that the solution may be clar- 
ified more easily. 

Sometimes in practice, we may have no choice of solvent for a given polymer. 
For example, poly(l,4-phenyleneterephthalamide) (PPTA or Kevlar) is only sol- 
uble in very strong acids which are viscous. In such cases, ultracentrifugation 
instead of filtration has to be used to remove dust particles from the solution 
[ 14-16, 25] . As for copolymers, the selection of proper solvents is even more dif- 
ficult, because at least two solvents which satisfy the above mentioned three cri- 
teria are needed. For this reason, reported characterization of copolymers is 
quite limited [26-29]. 

One of the challenging problems in the application of LLS methodology is to 
study static and dynamic properties of fluorocarbon polymers, such as Teflon 
and Tefzel (registered trademarks of Du Pont), which defy the characterization 
owing to their insolubility in ordinary solvents. It is this unique solubility that 
makes these fluorocarbon polymers very useful in many applications, but great- 
ly annoys those who wish to explore their solution properties. However, it was 
found that some solvents were capable of dissolving them at high temperature, 
but the finding was not enough to solve all characterization problems because 
the technique for clarification and measurement of the solutions at high tem- 
perature remained to be established. 



3.2 

High-Temperature Spectrometer 

In order to characterize a polymer soluble only at high temperature, many diffi- 
culties had to be overcome before light-scattering measurement was actually 
made. Thus, a special high-temperature LLS spectrometer was first designed and 
developed at State University of New York in Stony Brook, and then its technique 
was transferred to both Du Pont and BASF [30-35]. 

Figure 4 shows a schematic diagram of the high-temperature LLS spec- 
trometer at Du Pont. In it, a thermally-controllable plate ( 1 3 ) used as a heat sink 
isolates the rotary table (12) from the outer thermostat (3) by means of two sets 
of stainless steel standoffs (14). The outer brass thermostat (3) is isolated from 
the room with 0.5-in. -thick porous silicone rubber. This arrangement creates 
an oven that allows the temperature to keep easily in the 320-340 °C range with 
temperature fluctuations of less than 0.2 °C. A glass (vacuum) jacket of 2.25- 
in. o.d. isolates the inner thermostat from the oven. The vacuum jacket reduces 
the temperature gradient in the light scattering cell. The inner thermostat has 
a separate temperature controller and a miniature platinum resistance ther- 
mometer that can be connected to a digital voltmeter through the vacuum jack- 
et by means of ceramic feed throughs. With this design, short-term (20 min) 
control of ±0.05 °C, intermediate-term (60 min) control of ±0.1 °C, and long- 
term control of ±0.5 °C can be achieved at 340 °C even in the absence of a vac- 
uum. Long-term temperature stability depends partially on room-temperature 
fluctuations even in the presence of the outer thermostat and the isolation 
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Fig. 4. Schematic top and side views of the high-temperature thermostat and detection system 
of the light scattering spectrometer: ( 1 ) silicon rubber insulation; (2) heating wires for the brass 
thermostats; (3) outer brass thermostat with fluid circulation facilities; (4) vacuum glass jack- 
et for thermal isolation made of precision polished glass of 2.25-in. o.d. with Kovar seals at 
both ends of a stainless-steel temperature shield with precision polished glass of 2.25-in. o.d.; 
(5) inner brass thermostat, which has a separate temperature controller and a thermometer 
and can accommodate a light scattering cell up to 27-mm o.d; (6) Glan-thompson polarizers; 
(7) fluid circulation paths; (8) lens; (9) field aperture; (10) optical fiber bundle; (11) rotating 
plate for multiple detectors; (12) RT-200 Klinger rotary table with 0.01° step size; (13) cooling 
plate to isolate the outer thermostat from the rotary table; and (14) stainless steel standoffs for 
thermal isolation 



between the two thermostats. Figure 5 shows typical temperature fluctuations 
of the inner thermostat at 340 °C. To mention more, a high-temperature LLS 
detector coupled with GPC has recently been developed and the determination 
of the molar mass distribution of poly([phenylene sulfide) in 1 -chloronaphtha- 
lene at 220 °C has been made possible with it [36]. The advantage of on-line 
coupling LLS with GPC is obvious, since GPC is a fractionation method and 
LLS allows an absolute molar mass measurement and hence makes the cali- 
bration of GPC columns. 
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Fig. 5. Typical temperature fluctuations in the inner thermostat at 340 °C. Intermediate-term 
(1 h) temperature fluctuations were ±0.1 °C. Long-term (10 h) temperature fluctuations were 
± 0.5 °C 




3.3 

Solution Preparation at High Temperature 

In order to prepare and clarify a polymer solution at a temperature higher than 
200 °C, two different apparatuses were developed [30, 33]. Importantly, they are 
able (1) to dissolve a polymer under an inert atmosphere without losing solvent 
and without building up inner pressure due to solvent evaporation, and (2) to 
transfer the solution into a filtration device by a remote control. 

Figure 6 shows a specially designed dissolution/filtration apparatus which can 
be placed inside a small oven. Known weights of a polymer sample and a filtered 
solvent, as well as a small glass-enclosed magnetic stirrer, are placed in A at room 
temperature. The solution vessel (A) is then connected to the precleaned filter 
(B). After degassing followed by introduction of nitrogen, both stopcocks are 
closed and the oven is heated to the desired temperature to dissolve the polymer 
while the solution is stirred. When the polymer is considered to have been com- 
pletely dissolved, the solution vessel is turned 180° by means of the seal glass joint 
(J). This allows the polymer solution to be transferred from A into B without 
exposure to air. A gentle pressure using nitrogen is applied to force the polymer 
solution to pass through the fine-grade sintered glass filter (f) and to move 
directly into the precleaned dust-free cylindrical light scattering cell (C). In this 
way, dust-free polymer solutions can be successfully prepared, keeping the tem- 
perature high. 

Figure 7 shows another type of dissolution/filtration apparatus. Sleeve A (with 
no bottom) is joined to the shaded stopper, which is connected to a reflux con- 
denser by means of a greaseless glass joint. Cup (B), with a magnetic stirrer (E) 
sitting on top of the fine-grade sintered glass filter (F 2 ), allows solution clarifi- 
cation. First, a dust-free solvent and a polymer are placed in B of the argon-filled 
apparatus which has the preattached dust-free light-scattering cell (D). The 
reflux condenser flushed with argon is then inserted taking care that the appa- 
ratus with the polymer and solvent is under an inert atmosphere at room tern- 
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Fig. 6. High-temperature dissolution/filtration apparatus. The entire apparatus is placed in a 
high-temperature oven controlled at 250 ± 2 °C. (A) Solution vessel, where known weights of 
polymer and solvent as well as a small glass-enclosed magnetic stirrer are introduced. (B) a 
fine-grade sintered glass filter, connected to A and C by means of clean seal glass joints (J) 
(14/20, Wheaton Scientific), (f) Fine-grade sintered glass filter. (C) Cylindrical light scattering 
cell of 27-mm o.d. with a clean seal glass joint 



perature. The entire apparatus is set in a small oven and the temperature is raised 
to the desired value. When the polymer is completely dissolved, an argon pres- 
sure is applied through the sintered glass filter (Fi) to let the polymer solution 
move from B to C. The additional pressure difference plus gravity will filter the 
polymer solution directly into the dust-free light-scattering cell (D). After the fil- 
tration process is completed, the additional argon pressure is released. The stop- 
cock above the upper filter is closed during a light-scattering experiment. Briefly, 
extreme care must be taken when the test solution is prepared and subjected to 
light scattering measurements at temperatures near the boiling point of the sol- 
vent because the pressure build-up could cause an explosion. This points to the 
absolute necessity of having a pressure releasing mechanism. 
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Fig. 7. Separate components of a high-temperature dissolution/filtration apparatus. The 
assembled apparatus, as shown in Fig. 1 of Ref. 30, is placed in a high-temperature oven. (A) 
Cylindrical insert (with on bottom) with a diameter ~2 mm less than that of the solution ves- 
sel (B); (C) Filtration section with a fine-grade sintered glass filter (Fi) and ground glass joints 
to a light scattering cell (D) and a ground glass joint adapter for the water cooled condenser 
which is located outside the temperature-controlled oven. (E) Magnetic stirring bar. The shad- 
ed area denotes volume reduction so that the volume accessible by vapour phase is no more 
than a few times the fluid phase. The miniature water-cooled condenser has a coarse-grade 
sintered-glass filter (Fa) so that the entire system is always isolated from external dust. The 
greaseless stopcock above Fa is for operating the apparatus as a closed system, for introduc- 
tion of low vacuum in order to degas the solvent before dissolution, for filling the apparatus 
with inert gas, such as argon, in order to cdleviate chemical decomposition, and for releasing 
a possible pressure build-up at high temperature if chemical decomposition takes place. The 
entire apparatus is portable and can be inserted into the high-temperature light-scattering 
spectrometer with the light-scattering cell (D) and part of the filtration component (C) con- 
trolled at a given high temperature 
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3.4 

Differential Refractometer 

One of the most important parameters in static LLS is the specific refractive 
index increment (dn/dC), defined as /fmc^o(9n/dC)T,p;t. Since this parameter is 
not an intrinsic property of the polymer, the conditions of fixing temperature T, 
pressure P and wavelength of light in vacuum A are needed in its definition. Note 
that, according to Eq. (2.1 ), an error of E% in dn/dC will lead to an error of 2E% 
in the derived Mw 

The refractive index increment Dn is normally measured by using either a dif- 
ferential refractometer or an interferometer. In the former, light is refracted at 
the boundary between the sample and a reference liquid. Commonly, the beam 
displacement is measured directly and then converted to An by multiplying by a 
calibrated constant which can be obtained by using a solution with an accurate- 
ly known refractive index difference An [37]. This is not an absolute method since 
the constant has to be calibrated under the same conditions as in the light scat- 
tering measurements. In the latter, two light beams with identical geometrical 
paths traverse two different optical paths. One passes through the sample and 
the other through the reference liquid. This method relies on the interference of 
the two beams. Its details can be found elsewhere [38, 39] . In a high-temperature 
LLS measurement, a deformed cylindrical light scattering cell is preferred to the 
conventional divided differential refractometer cuvette in which the exit laser 
beam is refracted by the solution/air interface [30]. 

Figure 8a shows the design of a recently developed and commercially available 
refractometer (ALV GmbH, Langen, Germany). A small pinhole (P) with a diam- 
eter of 400 pm is illuminated with a laser light. The illuminated pinhole is imaged 
to a position-sensitive detector (PD) (Hamamatsu S 3932) by a lens (L) located at 
an equal distance from the pinhole and the detector, where the distance is four 
times the focal-length (f= 1 00 mm) of the lens. Thus, this novel refractometer uses 
a (2f-2f) design instead of a conventional (If) design which uses parallel incident 
light beams and makes the distance between the detector and the lens equal to 
only one focal length. A temperature-controlled refractometer cuvette (C) (Hell- 
ma 590.049-QS) is placed just in front of the lens. It is a flow cell and has a volume 
of ~20 mL, which is divided by a glass plate at -45° into two compartments. The 
pinhole, the cuvette, the lens and the detector are rigidly mounted on a small opti- 
cal rail. The refractometer has dimensions of only 40 cm in length, 1 5 cm in width 
and 10 cm in height, and the length can be easily reduced to 20 cm with another 
lens if necessary. The output voltage (-10 to 10 volts) from the position-sensitive 
detector is proportional to the displacement of the light spot from the center of 
the detector, and can be measured by a digital voltmeter or an analog-to-digital 
data acquisition system and a personal computer. 

Figure 8b shows the basic principle and the light path of the refractometer, 
where 0', 0", 0"', and the cuvette are drawn enlarged to make the details clear. If 
both the compartments are filled with a solvent (i.e. n=no), the illuminated pin- 
hole will be imaged at point O. However, if the solvent in one of the compart- 
ments is replaced by a dilute polymer solution with a slightly different refractive 
index (i.e. n=no-i-An), the light will be bent first by the glass plate, then by the 
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Fig. 8 a. Schematic view of a novel differential refractometer (produced by ALV GmbH, Lan- 
gen, Germany), which consists of a pinhole (P), a differential refractometer cuvette (C), a lens 
(L, f = 10 cm) and a position sensitive detector (PD). All components are rigidly mounted on 
a 40 cm long optical rail, b Light path in which one compartment of the cuvette contains a sol- 
vent with refractive index n and the other contains a solution with slightly different refractive 
index n = no -i- An. The cuvette and angles 6', 6" and 6"' (actually very small, ~ 0,0 1 radian) are 
enlarged to make the light path distinct 



cuvette wall and finally by the lens. The image is shifted away from the point o 
by a distance of Y. Figure 7b shows that 

Y=Yi-t-Y 2 +Y 3 =c tan(0')-l-(2f-X-c)tan(e")+2f tan(0"') (3.1) 

and 

f tan(0")=f tan(0'")-i-c tan(0')-i-{2f-X-c) tan(0") (3.2) 

where c, X and 0 are constants. Snell's law gives 

noSin(9O-0)=(no+An)sin(9O-0-0') (3.3) 

and 

(no+An)sin(0')=sin(0") (3.4) 

where 0', 0" and 0"' are actually so small because An is in the order of 10“^ RI 
units that we may set sin(0')=0', sin(0")=0", tan(0')=0', tan(0")=0" and 
tan(0"')=0"'. Combination of Eqs (3.1)-(3.4) leads to 

Y=KAn (3.5) 

where K=[X-i-c"(l-l/no)]tan(9O°-0).For agiven optical set-up and solvent, X, c, 
0, no, and hence K are constants. Equation (3.5) shows that the signal is propor- 
tional to An, and the larger the value of X, the higher the sensitivity (Y/An) is. 
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This means that the cuvette should be placed as close as possible to the lens in 
the experimental set-up. 

In the (2f-2f ) design, the detector and the pinhole (acting as a light source) are 
placed at the exact imagining positions along the optical axis of the lens. This con- 
figuration is optically equivalent to placing the detector directly behind the pin- 
hole, so that the laser beam drift is eliminated. In comparison with the conven- 
tional differential refractometer, this novel design has made the measurement of 
An much easier and provides reliable and accurate values for dn/dC since it is sta- 
ble and all the results can be recorded and averaged on a computer. Figure 9 shows 
the concentration dependence of An for a 13% PET-PCL copolymer in three dif- 
ferent solvents. The lines represent the least-square fits to the data points. 

The refractometer with its present dimensions can be easily installed into any 
existing laser light-scattering spectrometer together with the laser source, the 
thermostat and the computer, as exemplified in Fig. 10. The optical glass plate 




Fig. 9. Possible arrangement of the novel differential refractometer with an existing laser light 
scattering spectrometer 



Incorporatiun of a differential refractometer 
with a laser light-scattering spectrometer 




Fig. 10. Concentration dependence of the refractive index difference (An) between the poly- 
mer solution and solvent for a 13% PET-PCL copolymer. The lines represent the least-square 
fits to the measured data 
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placed in the laser light path at 45° reflects laser light by about 4%, and the reflect- 
ed light is used as the light source. With this design, the light scattering and the 
refractive index increment can be simultaneously measured under the identical 
experimental conditions of wavelength and temperature. The details of this nov- 
el spectrometer have been described elsewhere [40]. 



4 

Data Analysis 

The methods of analyzing data for the concentration and angular dependence of 
the time-average scattering light intensity and the intensity-intensity time cor- 
relation function can be found in many LLS books and related literature. In this 
section, we will mainly concern ourselves with how to combine static and 
dynamic LLS results to characterize special polymers in regard not only to the 
average molar mass, but also to the molar mass and composition distributions. 

4.1 

Conversion Between Translational Diffusion Coefficient 
Distribution and Molar Mass Distributions 

Though not involving fractionation, DLS allows estimation of the molecular 
weight distribution of a polymer. The principle is as follows. For a polydisperse 
polymer consisting of n homologous species, it is well known that G(Di) for 
species i at vanishingly small C and q is given by 

n 

G(Di) = MjWi/^MjWj (4.1) 

i=l 

where Wj denotes the weight fraction of species i with molecular weight Mi. For 
a continuous distribution of molecular weisht this eives 

G(D)dD = Mf„(M)dM/M„ (4.2) 



where fw(M) denotes the weight distribution of molecular weight M. Thus, we 
get 

f„(M) = (M„/M)G(D)(dD/dM) 



Empirically we have for a series of homologous polymers [41] 

D = kijM““° 

where k^ and are constants. Experimentally, for a flexible polymer, 
0.5<aD<0.6 in a good solvent and aD=0.5 in a Elory © solvent; for a rigid rod- 
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like chain, aD=l; and for a semi-rigid worm-like chain, 0.6<an<l. With Eq. (4.4), 
Eq. (4.3) is transformed to 



f„(M) = koM„M 



(4.5) 



which indicates that fw(M) can be determined if G(D) is obtained by Laplace 
inversion of Eq. (2.6) and the values of the parameters k^ and are available 
from separate sources. This is the basic idea of the method which allows infor- 
mation about fw(M) to be derived by DLS. 

Since, as noted above, the success in determining G(D) is actually not in the 
choice of a computer program for Laplace inversion but reducing the noise lev- 
el in measured g^^^(t,q). Thus, it is crucial that the solution is cleaned (i.e. "dust- 
freed") very thoroughly before it is subjected to laser light scattering measure- 
ments. Eor example, in studies conducted by the author, efforts were made to 
ensure that the relative difference between the measured and calculated base- 
lines did not exceed 0.1%. The error analysis related to the above problem can 
be found elsewhere [42, 43]. 



4.2 

Scaling of Translational Diffusion Coefficient D with Molar Mass M 
4.2.1 

Using a Set of Narrowly Distributed Standards 

The most straightforward method for calibrating the relationship between D 
and M is to measure both D and M for a set of monodisperse samples with dif- 
ferent molecular weights. In reality, the monodisperse samples have to be 
replaced by narrowly distributed standards made available either by relevant 
living polymerization or by fractionation of a broadly distributed sample. 
However, only a few kinds of polymers, e.g. polystyrene and poly(methyl 
methylacrylate), can actually be prepared so as to have a sufficiently narrow 
molecular weight distribution (M„/M,,~l.l),and the fractionation is very time 
consuming. Thus, the straightforward calibration of the D vs M relation is not 
always practical. 

Eigure 11 shows the plot of log(D) versus log (M) for a set of narrowly dis- 
tributed polystyrene standards in toluene at 20 °C [44] . The line represents a least- 
square fitting to the data points, giving D(cm/s)=3.64xlO“^M Using this 
relation, we were able to estimate the molar mass distribution of polystyrene by 
using only DLS [45] . However, as noted above, in reality, only a few kinds of poly- 
mers can be prepared to have narrow enough distributions of molecular weight. 
Hence, we often have to satisfy ourselves with more broadly distributed samples 
having different average molecular weights. This means that special analytical 
methods have to be developed to calibrate or scale the translational diffusion 
coefficient D and molar mass M by using from broadly distributed samples. Ideas 
on this theme are described in the following sections. 
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Fig. 1 1 . Double logarithmic plot D vs M for polystyrene in toluene at 20 °C. The line represents 
a least-square fits to the data points, giving D(cm/s) = 3.64 x (Ref. 42) 



4.2.2 

Using Two or More Broadly Distributed Samples 

As can be easily shown, it follows from Eqs.(4.4) and (4.5) that 

(4.6) 

(4.7) 



=(ko)''“'’ / J~G(D)D'^“ME 

The quantity on the right-hand side is denoted by . i.e. 

Mw!4ic =(kD)'^“° / j”G(D)D^^“^dD 



For any given set of and ap, it can be calculated when G(D) is determined 
from DLS measurements. We denote it for a polydisperse sample i by ); . 

Then we get 




poo * / poo 4 / 

Gi(D)D^'“°dD / Gj(D)D^'“°dD 
JO •’ Jo * 



(4.8) 



where G(D)j is for sample i. The right-hand side can be calculated from measured 
G(D)i and G(D)j for any ao, and if the chosen value of an happens to be equal to 
that for the system under study, the resulting value of (Mw,°aicd)i /(Mw,°aicd)j 
should agree with the value of (Mw)i / (Mw)j which can be determined by SLS. In 
reality, the desired will be reached by using a computer program which seeks 

a minimum of ERRORfan) defined by 
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ku / W 



Fig. 12. Plotof ERROR(kD) calculated with datafor five Chitosan samples with different weight 
average molecular weights. Here the minimum of the dashed curve corresponds to = 0.665 
± 0.015 and ko = (3.14 ± 0.20) X lO "* 



ERROR(ap) 



y y r (Mw)j 

t t 




(4.9) 



where N is the total number of the polydisperse samples examined. 

Next, with the ap value so determined, we compute for each of the 

N samples from Eq. (4.7) by varying ko and seek a ko value which minimizes 
ERROR(ko) defined by 



ERROR(kD) = |;[-^^;5|^i 1]^ (4.10) 

Eigure 12 depicts a plot of ERROR(ko) at three values of ap calculated from 
SLS and DLS data for five chitosan samples with different Mw It is seen that for 
each chosen an ERROR(kD) shows a sharp minimum,but the location and height 
of the minimum varies significantly with ap and the minimum becomes the 
smallest at aD=0.665 and kD=3. 14x 10 k With these values, the molecular weight 
distributions of chitosan samples have been successfully characterized [18]. 

4.3 

Combination of LLS with Other Off-line Methods 
4.3.1 

Intrinsic Viscosity 

If only one broadly distributed sample is available, we have to resort to another 
method to determine the relation between D and M. One of them is to estimate 
an from the Mark-Houwink equation for intrinsic viscosity. It is known that the 
intrinsic viscosity [q] can be scaled with M by the Mark-Houwink equation, i.e. 
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M, / (g/mol) 

Fig. 13. Comparison of cumulative weight distributions F„^(;um(M) [= Jm ] 

obtained by DLS and SEC (size exclusion chromatography) for a linear poly- 
ethylene in 1,2,4-trichlorobenzene at 135 °C 

[r\]=k^M“n, and according to Flory and also to de Gennes [41, 46], ap may be 
equal to=(a,^+l)/3. With estimated from by this relation, Mw from static 
LLS and G(D) from dynamic LLS can be used to determine as described 
above. Chu et al. [47,48] successfully applied this method to linear polyethylene 
in 1,2,4-trichlorobenzene at 135 °C for which ac was estimated from the report- 
ed value of ar^=0.72. [49] 

Figure 13 shows the resulting cumulative weight distribution Fw,cum(M) 
[=1^ fw(M) dM] andcomparesitwiththeresult obtained by high-temperature 
size exclusion chromatography (SEC). The agreement is reasonably good, but it 
should be noted that the weight distribution from LLS is usually narrower and 
more skewed toward the high molar mass than that from GPC because the scat- 
tered light intensity is proportional to the square of molar mass so that higher 
molar mass species weigh more in LLS. 

4.3.2 

Gel Permeation Chromatography 

The static laser light scattering apparatus used as an on-line GPC detector has been 
popular for a while. Here, we illustrate another but less known method of com- 
bining the results from (gel permeation chromatography) and DLS. The basic prin- 
ciple is as follows: There is a similarity between these two tools in that the trans- 
lational diffusion coefficient D obtained by DLS and the elution volume V in GPC 
are related to the hydrodynamic size of a given macromolecule. In a first approx- 
imation, if the hydrodynamic size is proportional to the molar mass, we have 

V=A-bBlog(M) (4.11) 

where A and B are constants similar to ko and ap. It should be noted that this 
approximation simplifies but does not affect the following discussion. Combin- 
ing of Eqs. (4.4) and (4.11) leads to 
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V=A+Blog(D) (4.12) 

where A =A+B log(kD)/an and B =-B/an. Furthermore, we get from Eq. (4.12) 

V2 = A2+2AB log(D) + BMog2(D) (4.13) 

Averaging both sides of Eqs. (4.12) and (4.13) over the concentration profile C(V) 
in GPC, we obtain 

<V> = A+B <log(D)> (4.14) 

and 

<V^> = A^+2AB <log(D)>+B^ <log^(D)> (4.15) 



where 

l»oo 

< V > = VC(V) dv 

Jo 

and 

<V^>=r V^C(V) dV (4.16) 



while 

r log(D)C(V) dV 
< log(D) > = 

r C(V) dV 

Jo 

and 

r log^(D)C(V) dV (4.17) 

< log2 (D) > = ^ 

r C(V) dV 

•'O 

Though not yet theoretically proved, it is usually assumed that the differential 
area C(V)dV under a GPC curve is proportional to the differential mass of the 
polymers dW that are contained in the differential elution volume dV. Since dW 
oc fw(M)dM, we have 

C(V)dV oc f„(M)dM (4.18) 

If C(V) is normalized, this gives 

C(V)dV=fw(M)dM (4.19) 

Combining Eqs.(4.3), (4.4) and (4.19), Eq. (4.17) can be rewritten 

r log(D)G(D)D^^“‘’dV 

< log(D) > = — 

r G(D)D^^“‘>dV 
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and 



r log^(D)G(D)D^^“°dV 
< log2 (D) > = 

r G(D)D^'“‘’dV 

Jo 



(4.20) 



After Eq. (4. 1 9 ) is multiplied by M, both sides are integrated over the entire range 
of M to yield 

l»oo 

M„=£MC(V)dV (4.21) 



Elimination of M from the right-hand side using the relation D=kDM and 
Eq. (4.12) transforms Eq. (4.21) to 

dV (4.22) 



which is combined with Eq. (4.6) for Mw to give 



r io(A-V)/(aDB)c(v) dV rG(D)D^^“MD 

Jo 



= 1 



(4.23) 



This equation contains only one unknown parameter ap. For a chosen an, we 
first calculate <log(D)> and <log^(D)> using Eq (4.20), then calculate A and B 
by solving Eqs.(4.14) and (4.15) with <V> and <V^> computed from the GPC 
diagram, and finally calculate the left side of Eq. (4.23). By iteration, we can find 
a value of an which may minimize the difference between the left and right sides 
of Eq. (4.23). For the ap so obtained, we can calculate ko from either Eq. (4.6) 
or (4.22) by using Mw determined directly from static ELS and C(V) from SEC 
or G(D) from dynamic ELS. With A, B, kp and ap, we are ready to calculate A 
and B. In this way, we can calibrate not only the M vs V relation, but also the M 
vs D by a single process on only one broadly distributed sample. This method 
has been tested and applied in the characterization of gelatin in water [50, 51]. 



5 

Applications 

When the relation between D and M is established, we can easily convert G(D) 
obtained by dynamic ELS into a differential molecular weight distribution, such 
as fw(M). We have successfully applied the above methods to various kinds of 
polymeric and colloidal systems, such as for Kevlar [15, 23], fluoropolymers 
(Tefzel & Teflon) [12,30-35, 52], epoxy [53-55], polyethylene [56, 57], water-sol- 
uble polymers [18, 50-51, 58, 59], copolymers [60-62], thermoplastics [63-65] 
and colloids [66-72]. Three of those applications are illustrated below. 
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5.1 

Segmented Copolymers 

We consider a copolymer sample consisting of monomers A and B. The sample is 
generally polydisperse in both molecular weight and chain composition. We sup- 
pose that the copolymer species i is characterized by the molecular weight Mj and 
the composition wai which is the weight fraction of A in the chain of that species. 
It is assumed that no composition heterogeneity exists in the chains of the same 
length. For homopolymers the refractive index increment (at infinite dilution) 
does not depend, in a good approximation, on the molecular weight of the chain, 
but is equal to that of the entire sample. For copolymers this is not the case, and, 
according to the theory of light scattering, Eq. (4.1 ) for Gi may be replaced by 

Gi=(Vi/v)^MiWj/;^MjWj (5.1) 

j 

where v; is the refractive index increment due to the copolymer species i, v that 
of the entire sample, and w; the weight fraction of the copolymer species i. When 
the molecular weight distribution maybe treated as continuous. Equation (5.1) 
is generalized to 



G(D)dD = (v(M)/v)^Mf„(M)MdM (5.2) 



where, as before, fw(M) denotes the weight distribution of M, and v(M) is the 
refractive index increment due to the chains of molecular weight M. Note that 
v(M) depends on wa(M), which is the continuous version of wai. We assume that 
this dependence is represented by 

v(M) =VaWa(M)-i-Vb[1 -wa(M)] (5.3) 

where Va and Vb are the refractive index increments of the homopolymers con- 
sisting respectively of A and B. 

Corresponding to Eq. (4.3) for homopolymers, we introduce fw,app(M), an 
apparent weight distribution function of M,by 

fw,app(M) = (M„ /M)G(D)(dD/dM) (5.4) 



With Eq. (4.4) (assumed to hold for copolymers too), this gives 

fw,app(M)/M^ =(l/kD)''““ G(D)d1+(2/“°^ (5.5) 



Therefore, fw,app(M)/Mw can be calculated from DLS determination of G(D) 
along with Eq. (4.4) when ku and are known separately. On the other hand, 
substituting Eq. (5.2) together with Eq. (5.3) into Eq. (5.4), we get 
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fw,app(M)/M„ =V ^[VAWA(M) + VB(l-WA(M))f f„(M) (5.6) 

Now, we choose two solvents 1 and 2 for a given copolymer. Since fw,app(M), v, 
Va, and Vb vary with solvent, their values in solvent i are denoted by fw,app(M)^'\ 
v^, Va^'^ and Vb^'^. Then it follows from Eq. (5.6) that 

^w,app (^) I Wa(M)Va^ ^ + [1 ~ Wa(M)]Vb^ ^ I 
fw,app®(M) [v® Wa(M)Va® + [1-Wa(M)]Vb^'^| 

The ratio on the left-hand side can be obtained as a function of M since 
fw.app(M)®/Mw can be determined, as described above, from experimental infor- 
mation. Thus, Eq. (5.7) allows determination of wa(M), the chain composition 
distribution when all other parameters on its right-and side are measured by dif- 
ferential refractometry. Once wa(M) is known, we are ready to compute v(M) 
from Eq. (5.3), fw(M) from Eq. (5.6), and finally Mw [59]. 

Figure 14 shows PLS determined - G(D)s by Eq. (5.3) for low-mass (O) and 
high-mass (□) segmented copolymer poly(ethylene terephthalate-co-caprolac- 
tone)s (PET-PCL) containing 13% PET in tetrahydrofuran (THE) at 25 °C. 
Repeating the measurements on these samples in another solvent chloroform 
should lead to a new set of G(D)s, which allows Eq. (5.7) to be used to calculate 
wpet(M). 

Figure 15 shows the results from such calculations for low mass (O) and high 
mass (□)13% PET-PCL samples. We see that the PET content increases with 
increasing M for M< ~4x 1 0"^ and approaches a constant value of ~ 1 4% in the high 
molar mass range. For the 58% PET-PCL sample, the composition distribution 
is nearly constant. The composition of the high-molar mass 13% PET-PCL sam- 
ple overlaps with that of the low mass 13% PET-PCL sample in the same molec- 




Fig. 14. Apparent weight distributions calculated from the translational diffusion coefficient 
distributions corresponding to low-mass (O) and high-mass (□) copolymer segmented 
poly(ethylene terephthalate-co-caprolactone) (PET-PCL) containing 13% PET in tetrahydro- 
furan (THE) at 25 °C 
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Fig. 15. Estimate of the chain composition distributions for low mass (O) and high mass 
(□)13% PET-PCL samples 



ular weight range. This indirectly indicates that the estimation of the composi- 
tion distribution is reasonable. The lower content of PET in the low molar mass 
range can be attributed to the two-step synthesis [59-62]. 

5.2 

A Polymer Mixture Containing Individual Chains and Clusters 

If a mixture is made of individual polymer chains and polymer clusters, the mea- 
surement of static LLS will lead to an apparent weight-average molar mass Mw,app 
which is expressed by 



Mw.app ~ (5.8) 

where the subscripts "L" and "H" denote low molar mass linear polymer chains 
and high molar mass polymer clusters, respectively, and xp and xh are their 
weight fractions with xl+xh= 1. If the linear chains and clusters are significant- 
ly different in the hydrodynamic size, dynamic LLS will detect two distinct peaks 
in the measured G(D), with one peak corresponding to the linear chains and the 
other to clusters. 

Figure 1 6 shows G(D) of a simulated polymer mixture at two scattering angles 
("O", 14° and 17°). The mixture consists of two polystyrene standards hav- 
ing distinctly different weight average molar masses (3.0x 1 0^ and 5.9x 1 0® g/mol) 
and a high mass polystyrene which is used to simulate the polymer cluster [66]. 
The area ratio Ar of the two peaks is expressed by 

A,=^. = ,5.,) 

r Gh(D) dD ^w,H^H 
Jy 

with Y the cutoff translational diffusion coefficient between Gp(D) and Gh(D). 
In practice, the values of Ar at finite scattering angles must be extrapolated to 
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q^O. Figure 17 shows this extrapolation of the for the two peaks in Fig. 16. 
With Eqs. (5.8) and (5.9), Mw.app from static LLS and Ar from dynamic LLS allow 
and Mw,hXh to be computed. In principle, by knowing any one of the four 
parameters (Mw,l> Mw.h> xl and xh), we should be able to determine the remain- 
ing three parameters. This method has been thoroughly tested with mixtures of 
polystyrene standards [64] . As for a particular polymer mixture, we should find 
a way allowing determination of one of the four parameters. For example, in the 
study of polymer association, we can determine the Mw,l of starting individual 
polymer chains, and in the study of the gelation process, we can use a filtration 
method to remove large microgels, so that the weight fractions of xp and xh can 
be subsequently determined. This method has been used to characterize novel 
thermoplastic polymers with phenolphthalein in their backbone chains [61-63]. 




Fig. 16. Translational diffusion coefficient distributions G(D) of a simulated polymer mixture 
at two scattering angles (“□”, 17° and“0”, 14°). The mixture contains two polystyrene stan- 
dards of distinctly different weight average molar masses (3.0 x 10^ and 5.9 x 10® g/mol) and 
a high mass polystyrene 




Fig. 17. q^-dependenceof the area ratio (Ar = Al/Ah) for the two peaks of the translational dif- 
fusion coefficient distribution shown in Fig. 16 
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5.3 

Polymer Colloids 

Combination of static and dynamic laser light scattering is also useful to deter- 
mine not only the size distribution but also the particle structure of polymer col- 
loids such as the adsorbed surfactant layer thickness [73] and the formation of 
nanoparticles [74, 75]. A recently developed method of determining the density 
of polymer particles is outlined below to illustrate the usefulness of laser light 
scattering as a powerful analytical tool for investigating more sophisticated col- 
loidal problems [76-78]. 

For a colloidal particle of uniform density its molar mass M is proportional 
to the cube of its radius R, i.e. 



M=(4/3)tiR39Na (5.10) 

where 9 is the particle density and Na the Avogadro constant. The diffusion coef- 
ficient D of the particle (at infinite dilution) is related to the Stokes radius Rh by 

D = (kBT/6Titx)(l/Rh) (5.11) 

where kfi is the Boltzmann constant, T the absolute temperature, and ri the sol- 
vent viscosity. We assume that Rh is larger than R by the thickness b of the sol- 
vated layer, i.e. 



Rh = R + b (5.12) 

Substituting this into Eq. (5.1 1) and using Eq. (5.10) for R, we obtain 

D = F(kBT/6Ttq) (4 ti9Na)‘^^ (5.13) 



with 



F = l/[H-b(4Tr9NA/M)i/3] (5 j4) 

Comparing Eq. (5.13) with the relation D=kDM and considering b<<R, we 
find approximately 

a =1/3 (5.15) 

kD = F(kBT/6Trq)(47X9NA)^^^ (5.16) 

Thus, with M in Eq. (5.14) replaced by M^, it follows from Eq. (4.6) that 



Mw = 



[l + b(47tpNA/M^)'^^f 



4ttpNA )(j^)3^| ( 5 . 17 ) 

3 67tq 




132 



ChiWu 



This equation contains two unknown parameters (b and 9), and if we know one 
of them, the other can be calculated from Mw and G(D). With this idea, it was 
found that the average density of the polystyrene microspheres made of a few 
uncrosslinked chains is slightly lower than that of bulk polystyrene or conven- 
tional polystyrene latex. 



6 

Conclusions 

This review has shown that static and dynamic laser light scattering (LLS) com- 
bined provide a very powerful method for polymer characterization. LLS has 
advantages over other polymer characterization methods, which include ultra- 
centrifugation and chromatography, in such features as speed, non-perturbation 
and extreme dissolution conditions (high temperature or strong acid). The most 
important advantage is that the calibration is independent of the particular LLS 
instrument used. However, the LLS method for the determination of mass dis- 
tributions described in this paper is disadvantageous in that its resolution is not 
as high as the fractionation methods, especially for samples whose mass distri- 
butions have closely packed peaks. The LLS method should play a definite role 
in circumstances where polymers intractable by conventional characterization 
methods have to be treated. In principle, dynamic LLS can be used together with 
other polymer characterization methods which take advantage of the depen- 
dence of the hydrodynamic volume on molecular weight. 
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The Smoluchowski coagulation equation derived back in 1916 is usually linked with the dif- 
fusivity and size of aggregating particles. It can also be used as a versatile tool for mean-field 
kinetic analyses. In this paper it is shown to be an efficient tool for studying the relationships 
between the reactivity of functional groups in monomers and the size distribution of polymer 
species in step and chain growth polymerizations. The Smoluchowski coagulation equation 
and its modifications are applied as models of kinetically controlled growth reactions (with 
irreversible elementary reaction steps). Advantages and limitations of this method of model- 
ing polymerization processes are discussed. 

Keywords: Size distribution, cross-linking polymerization, aggregation, classical modeling. 
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List of Symbols and Abbreviations 



A,B 

Ck 

d 

AE 

f 

H, G, E, A, B 
Ki,j 

ki,j 

Mk 

N 

Na 

nk 

P 

Pi 



r 

RA2,R’B2, 

ARB,RA^ 



Si 

t 

-Wk 

x,y,z 

a, (3, A, B, C 

K,kA,kB 

X 

W 

a 

T 

co,a, 0,T, z 



types of functional groups 
concentration of fc-mers 
space dimensionality 
activation energy 
functionality of a unit 
generating (counting) functions 

coagulation kernel, the rate constant for the reaction between i- 
andj-mer 

rate constant for the reaction between i- and j-mer 

fc-th moment of size distribution (k=0,l,2, ...) 

total number of units in the system 

the Avogadro constant 

number (mol) fraction of fc-mers 

conversion degree 

aggregate of size i, polymer molecule of polymerization degree 
i, f-mer 

number average degree of polymerization 
weight average degree of polymerization 

types of monomers (A and B denote functional groups) 

component of the product kernel 

time 

reaction time at the gel point 
number of ways of assembling a k-mer 
dummy variables in generating functions 
numerical constants 

ratios of rate constants, relative rate constants 
cyclization parameter, degree of homogeneity, 
contribution to activation energy 
symmetry factor 

units proportional to time, rescaled time 

parameters in generating functions beyond the gel point 

exponents 



1 

Introduction 

The modeling of a polymerization process is usually understood as formulation 
of a set of mathematical equations or computer code which are able to produce 
information on the composition of a reacting mixture. The input parameters are 
reaction paths and reactivities of functional groups (or sites) at monomeric sub- 
strates. The information to be modeled maybe the averages of molecular weight, 
mean square radius of gyration, particle scattering factor, moduli of elasticity, 
etc. Certain features of polymerizations can also be predicted by the models. 
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Among them is the gel point conversion, if multifunctional units are present, as 
well as accompanying divergence of viscosity, onset of equilibrium elasticity 
modulus, etc. By comparing the results of modeling with experiment, one can 
verify to what extent the chemistry is affected by physical interactions which are 
practically always active in polymerizations. 

Note that all the output parameters mentioned are directly related to system 
connectivity. The meaning of the word chemistry is in fact reduced to connec- 
tivity build-up since each chemical reaction between two functional groups 
leads to the formation of a link. 

This paper deals with one of the mean-field methods of modeling the con- 
nectivity build-up that can be applied to polymerization processes. As in the 
other mean-field methods of modeling, certain physical features such as con- 
centration fluctuations or fluctuation coupled diffusion control of reaction steps, 
etc., are neglected. 

The paper concentrates on the use of the Smoluchowski coagulation equation 
which is presented in Sect. 2. This equation is usually linked to the effects of dif- 
fusion and size of colloid particles in aggregation processes. Here, however, we 
refer to it as to probably the first known kinetic model applicable for the analy- 
sis of parallel-consecutive reactions. In Sects. 3-5 the Smoluchowski coagulation 
equation is shown to be a useful tool for kinetic modeling the step growth poly- 
merization of one or copolymerization of two bifunctional monomers. Section 6 
describes the use of the Smoluchowski equation to model the polymerization of 
multifunctional monomers. The equivalence of the model and those based on 
assemblage of structures from fragments using statistical arguments is dis- 
cussed in Sect. 7. An example of modeling a kinetically-controlled polymeriza- 
tion process which yields results differing from those obtained by using statisti- 
cal approaches is presented in Sect. 8. Post-gelation analysis which is available 
for the simplified systems only is described in Sect. 9. In Sect. 10 attempts at tak- 
ing into consideration intramolecular reactions are reviewed. Finally, in Sect. 1 1 , 
some general features of the Smoluchowski coagulation equation, mostly those 
loosely related to polymer science, are presented. 



2 

The Smoluchowski Coagulation Equation 

Consider an irreversible second order aggregation process where an aggregate 
of size (or mass) i and an aggregate of size (mass) j combine to form an i+j aggre- 
gate: 



(i,j = l>2,...) (1) 

The size of aggregates is an integer quantity. It is therefore convenient to con- 
sider it as a multiplicity of the size of the unit aggregate Pi, or, if applied to poly- 
mers, as the polymerization degree. Adopting polymer terminology, we shall 
refer to an aggregate of size i as an z-mer. 
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To characterize the process defined by Eq. (1), Smoluchowski applied the 
empirical [1] Guldberg-Waage [2] law of mass action. In his papers [3,4] pub- 
lished in 1916 and 1917 he derived and examined the set of differential kinetic 
equations equivalent to 



dcj 

dt 



-Cl 






i=l 



dc 



C2^K2jCi 



(2) 



dc 



QCo ^ 

-^ = KuCiC2-C3 2,K3,iCi 



i=l 



Nowadays, the set at Eq. (2) is known better in its compact form 

i=l 1=1 

which is called the Smoluchowski coagulation equation. 

Solutions of the Smoluchowski equation are of interest in all branches of 
physics where aggregation processes take place. There is a vast literature on its 
application including several reviews [5-7]. Here, we concentrate on those 
applications which can be used to model the polymerization processes. 

The Kij in Eqs. (l)-(3) is the rate constant of aggregation. The set of rate con- 
stants for all i and;, often called the coagulation kernel, is an infinite symmetric 
matrix with nonnegative elements. In order to keep the form of as simple as 
possible we shall assume that 

- the system is closed, i.e., the total number of units, IV, is constant, 

- the system is uniform with no concentration fluctuations affecting the reac- 
tion, 

- Ki^j depends only on i and;. 

The concentration of a k-mer, Ck, can be expressed in arbitrary units. It is conve- 
nient, however, to express its concentration as 



number of k-mers 

Ck = (. 4; 

N 

The set of cfs constitutes the distribution of aggregate sizes. The Smolu- 
chowski equation (Eq. 3) defines its change with time. 

In any kinetic analysis, the time, t, comes naturally as an independent vari- 
able, but in Eq. (3) it is, in fact, a variable merely proportional to the real time. It 
can be regarded as a rescaled time with the scaling factors depending on the actu- 
al units of both coagulation kernel and concentration as well as on the type of 
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polymerization reaction considered. Suppose we changed the definition at Eq. 
(4) of concentration to become simply the number of moles of k-mers. Then, the 
new time units would be related to the old ones with the scaling factor N/Na (Na 
is the Avogadro constant). 

The form of the Smoluchowski coagulation equation at Eq. (3) is significant. The 
first term on the right hand side describes the rate at which k-mers are formed from 
smaller components of matching sizes. The factor 1/2 takes care of the symmetry. 
It prevents the terms fCj_,CjC, and K^CiCj both being counted. For i=j, the 1/2 reflects 
the fact that formation rate of 2i-mer is only a half of the rate at which /-mers dis- 
appear from the system [8] . The second negative term of Eq. (3) is the rate at which 
k-mers disappear from the system in reactions of growth. Since the aggregation is 
irreversible, there is no term for the decomposition of k-mers into smaller compo- 
nents. Unless otherwise stated, we will consider formation of acyclic (tree-hke) 
aggregates only. In polymer parlance, intramolecular reactions will be disregarded. 

The definition at Eq. (4) implies that the zeroth moment of the k-mer size dis- 
tribution, Mq, is the reciprocal of the number average degree of polymerization 
in the system 

~ 1 

= ( 5 ) 

i=l n 

and the first moment provides the normalization of the distribution 

Mi = XiCi=l (6) 

i=l 

In the following sections it will be shown how the Smoluchowski equation, or 
rather its modifications, can be used to predict these two and higher moments 
of size distributions for various polymer systems. 



3 

Random Step Growth Homopolymerization of a Bifunctional Monomer 

This classical polymerization model [9] is sometimes referred to as an RA 2 or 
ARB system, with A and B denoting functional groups. The groups may react 
either with themselves {RA 2 system) or A groups may react with B groups only 
{ARB system). The ‘linear’ polymerization described in this section is in fact a 
special case of models 1 or 3 of Table 4 in Sect. 6, but it is dealt with separately 
because of its special importance to polymer science. 

The functional groups in reacting species (i-mers) are equally reactive irre- 
spectively of i. The reaction between an z-mer and j-mer can be written as 

aki- 

Pi+P^ ^Pi+j (i,j = l,2,..) 



( 7 ) 
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Since kij does not depend on i orj,it can be absorbed, together with the sym- 
metry factor o (equal to 4 for RA 2 or 2 for ARB), into time units and the Smolu- 
chowski coagulation equation would assume its simplest form with the kernel 
equal to one: 

^ = (8) 

i=l i=l 



The solution of Eq. (8) was presented by Smoluchowski himself [4] . The time 
dependence of the zeroth moment of the distribution is particularly easy to cal- 
culate from Eq. (8). Simply sum up the terms in Eq. (8) for all k’s to obtain 



dMp 

dt 






since 



/ \2 
00 

Vi=l 



k 

(A 

k=l i 



The solution of Eq. (9) is 



(9) 



( 10 ) 



Mo(t) = ^ (11) 

t-l-2 

A convenient method of extracting higher moments of the distribution direct- 
ly from Eq. (8) (to be fully exploited in following sections) is to define the poly- 
nomial generating function [10] containing concentrations of species, c,=c,(t), 
and the dummy variable x with no physical significance 

00 

H(t,x) = ^Cix‘ (12) 

i=l 

Now, multiplication of both sides of Eq. (8) by followed by summing up for 
all k’s yields 

^ = (13) 



The solution of Eq. (13) is 



H = 



4x 

(t+2)[t(l-x) + 2] 



( 14 ) 
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Note, however, that since the successive moments of distribution are simply 
given by 



(9H 

(9x 



= Mi; 



X=1 



dx^ 



= M 2 -Mi; 



lx=l 



dx^ 



-M3-3M2+2M1 



lx=l 



etc., we can convert Eq. (13) without even solving it into a set of ordinary differ- 
ential equations for the moments. Thus, by differentiating Eq. (13) with respect 
to X followed by substituting x=l we obtain 



d tdx 



„ dn dM, „ , 

H Ma ^ = 0=>Mi =1 

dx “o'x dt ' 



(16) 



since H(t,l)=Mo, 



d^Yi 

didx^ 






-hH 



dx^ 



■M„ 



^2H^ d(M2-Mi) 

dx^ dt 



= Mf => M2(t) = t-I-1 



(17) 



<9^H <9^H d(M3 - 3 M 2 +2Mi) 

=3— --hH Mo - = 

dtdx^ dx^ dx^ dx^ dt 

3Mi (M 2 - M, ) ^ M 3 (t) = i(3t^ + 6 t -H 2 ) 



etc. 

From the moments we can calculate the number, weight, and z-average 
degrees of polymerization which are defined thus: 





(19) 



respectively. The first averages of polymerization degree calculated from Eq. (19) 
are listed in Table 1. They are shown as functions of (rescaled) reaction time as 
well as functions of conversion degree, p. The latter is, as usual, defined as the 
fraction of all functional groups in the system that has reacted to bonds. The rela- 
tion between conversion and time is derived by noting that two functional 
groups have to react to form a bond and that each A:-mer has exactly k-1 bonds. 
Hence 
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Table 1. Number, weight, and z-average degrees of polymerization for the random homopoly- 
merization of bifunctional monomer calculated using the Smoluchowski equation (Eq. 8) and 
expressed in terms of time (t) and conversion degree (p) 



P = 




= Mi-Mo 



( 20 ) 



or 

p = — (21) 

t-b2 

By using Eq. (21) we find out that the moments of the distribution obtained 
by using the Smoluchowski equation are exactly the same as those derived for 
the same random homopolymerization using statistical arguments [9]. 

The explicit distribution of molecular sizes can now be obtained by solving 
recursively the Smoluchowski equation at Eq. (8). For C\ we get 



which instantly gives 



Similarly, by solving 



dCj 

"d? 



= -qMo 



Cl 



4 

(t-b2)^ 



dc2 

dt 



-^Cj C2M0 



(22) 



(23) 



(24) 



we obtain 



4t 

(t-b2)^ 



(25) 



etc. It is not difficult to show that carrying out with solving consecutive equa- 
tions we would obtain the entire distribution [11] 
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Cu =- 



4t 



k-i 



(t+2) 



k+i 



(26) 



(Obviously, the same distribution is obtainedby expanding Eq. (14) into a pow- 
er series around x=l.) This is the Flory most probable distribution. One can eas- 
ily verify it by replacing c/t with number (mol) fraction of k-mers, nu, by using 
the relation 



Ck-nkMfl ^27) 

Then, by replacing time with conversion using Eq. (21) we arrive at the famil- 
iar [9] 

Hk = a-p)p''“' (28) 

Thus, for the random homopolymerization of a 2 -functional monomer, the 
kinetic analysis leads to the same molecular size distribution as that derived by 
using statistical arguments. Conditions for such an equivalence will be discussed 
in the following sections 



4 

Homopolymerization of a Bifunctional Monomer Reacting 
with Substitution Effect 

Many monomers react with the so called substitution effect. To deal with it quan- 
titatively we shall adopt the approach of Gordon and Scantlebury [ 12] . Suppose the 
monomer molecule contains /similar fimctional groups and every group has, or 
has not, the same reactivity. The zeroth shell substitution effect is simply no effect 
at all. Each functional group reacts independently of whether or not other groups 
have reacted. If, however, the reactivity of remaining/- 1 groups changes as soon 
as the first has reacted, we have the first shell substitution effect (ESSE). Clearly, with 
ESSE the reactivity of remaining/- 2,/- 3,... unreacted groups may still change as 
the successive groups react. In other words, with ESSE the unit’s substitution 
degree determines the reactivity of its unreacted groups. An example of the ESSE 
that is particularly appeahng to chemists seems to be the way the primary amino 
group reacts. Initially, both amino protons have an equal chance of entering a reac- 
tion, but after one has reacted the reactivity of the other proton usually changes. 

The second shell substitution effect (SSSE) is operative when the reactivity of 
a group depends not only on the substitution degree of the unit it belongs to, but 
on the substitution degree of all units that are first neighbors to that unit. Unless 
stated otherwise, we shall assume for simplicity that all the unreacted groups of 
a unit have an equal chance of reacting. 

When an RA 2 monomer reacts with the substitution effect then the applica- 
tion of kinetics analysis will lead to a different molecular distribution to that 
resulting from a simple statistical analysis. This was convincingly demonstrated 
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by Kuchanov [6]. For the kinetics analysis one has to distinguish three types of 
reaction: 



P,+Pi 

Pi+Pj 

A\r 

P, +Pj 22 ^Ph 

i,j = 2,3,... 



(29) 



The analysis becomes simpler if we assume, following Gordon and Scantle- 
bury [13], that the contributions from reacting units to the activation energy of 
reaction between their reactive groups are additive. Let AE be the activation 
energy of the reaction between groups on unsubstituted units (monomers). 
When a monomer molecule reacts with a unit having one reactive group already 
reacted, let the activation energy change by the value of ip. Clearly, when each of 
the reacting units has one of its groups already reacted then the activation ener- 
gy will be AE + 2xp. With this assiunption we may absorb the constant 4^22 into 
the time units and deal only with the relative rate constant defined thus: 



K = 12/ = expl 

/k22 






'RT 



so that 



= KT 



(30) 



(31) 



'-22 



The Smoluchowski equation becomes now slightly more complicated. It can 
be written in the form 



dcj 



2 2 

= -K Cf -KC 



i=2 



1 

dt 2- 



K^Ci -/CC,C2 -C2 






(32) 



i=2 



dt 2 



k-2 



-Ck) 



i=2 



i=2 



When multiplied by x‘ and summed up over all k it becomes 



b>H 

dt 



2 



0 



( 33 ) 
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with 



where 



= JCCjX -l-^CjX* 

i=2 



Ho(t) = H(t,l) 



(34) 



(35) 



Note that the generating function in Eq. (33) without the tilde is the same as 
that given by Eq. (12) when x=l. The complete equation at Eq, (33) also 
becomes Eq. (13) when x=l,but it is difficult to solve analytically for any x. The 
moments of the distribution, however, can be calculated numerically as func- 
tions of time since 



H = H-l-(/c-l)ciX 



(36) 



By using the first line of Eq. (32) written in the form 



^ = -,cc,Ho (37) 

dt * ® 

one may extract the moments by successive differentiation of Eq. (33) with 
respect to x followed by putting x=l. 

Comparison of the weight average degrees of polymerization (equal to sec- 
ond moments) calculated by using Eq. (33) and a simple statistical model (cas- 
cade theory [14, 15]) is presented in Fig. l.The nonequivalence of the two meth- 
ods can clearly be seen. The number fractions of the monomer, once-, and dou- 
bly-reacted units were calculated from Eq. (33). The fractions were treated as 
probabilities for the units of being in their respective reaction states. The „sta- 
tistical“ weight average degree of polymerization was then calculated by using 
the probabilities following the cascade theory standard procedures [14, 15]. 



5 

Multicomponent System of Bifunctional Monomers 

An example of multicomponent system that can be dealt with by using kinetics 
equations of the Smoluchowski type is provided by a step growth alternating 
copolymerization of two bifunctional monomers [16]. This system requires lit- 
tle more laborious, but quite straightforward algebra. 

Let the reaction involve monomers RA 2 and R’B 2 - Functional groups A react 
with functional groups B. One can distinguish as many as nine elementary reac- 
tion steps. This is because each molecule contains either an even or odd number 
of units. Those with odd numbers of units may be terminated with either two 
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0 0,2 0.4 0.6 02 I 

Comnenim 

Fig. 1. The weight average degree of polymerization, P,y, in a linear polymerization with 
monomer reacting with a substitution effect. Solid line - calculated according to Eq. (33), bro- 
ken line - calculated by using cascade theory for the specified relative rate constant k 

groups A or two groups B and hence they may react differently. The monomers 
are degenerate cases of odd molecules. 

Table 2 summarizes all the reactions that may be distinguished in the sys- 
tem. On the right to each reaction there are shown codes of molecules that react 

Table 2. The list of reactions taking place during step growth alternating copolymerization of 
two bifunctional monomers 
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and that are formed. The codes are simply the pairs of numbers of units RA 2 
and R’B 2 in each molecule, respectively. This is not the only way of coding mol- 
ecules for the purpose of kinetics analysis. In fact the method of coding 
depends on the kind of information sought for in the analysis. Cheng and Chiu 
[17] were able to obtain detailed information on reacting molecules in a mul- 
tistage reaction system by using codes in the form of matrices or multicom- 
ponent vectors. 

Kinetics equations are then written for the rates of change of concentrations 
of k,k-mers {k,k-l mers and k-l,A:-mer). Equation (38) gives just the equation for 
molecules of even degree of polymerization (A:,A:-mers). The numbers of the cor- 
responding reactions in Table 2 are shown over each term. 



. (IV) (V) 

dci dcuu 

^ + 4kjQCQjC|; |;_j 4-2kj 



(II) (III) 

-2kQjCj_QCi^_l^ -2kjgCQj[Cij_i; -2kjjC|^ j 



(VI) 


(IX) 


k-1 


k-1 




k-j ^^Ti-l.jl-k-i+l.k- 


i=i 


j=2 


(VII) 


(VIII) 1 



CO CO CO 

Ss.i+X"io-i+Ss-io 

j=l j=2 i=2 



-II 



Hn 



(38) 



The specific rate of dimer formation (k=l), in reaction I in Table 2, has to be 
considered separately. It reads 



dc 



14 



dt 



- 4kooCi_oCoj 



form. 



(39) 



Before we present some analytical solutions for simplified cases let us con- 
sider the general case where monomers react with substitution effect. (An exam- 
ple of such a monomer might be, say, 2, 6-toluilenediisocyanate). In reactions I-IX 
there are four rate constants. We can define contributions to rate constants from 
the units, thus 

and k, = -^ (40) 

kii k„ 



The physical meaning of these constants is the following: kx (or ks) describes 
the reactivity of a functional group of the monomer relative to that of a group 
on a unit at the end of a chain. 

Furthermore, we can assume that the contributions to the activation energy 
of reaction from the substitution degree of reacting units are additive. Thus, the 
activation energy for the reaction between i substituted unit A and; substituted 
unit B (i,j=0,l) is 
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AEi,j=AEo,o+iV^A+jV^B 



(41) 



Accordingly, the ratio of rate constants can be written as the contribution 
product: 

^-k k 

k.. ~ ^ ^ (42) 



It is now only a matter of algebra to arrive at the following (master) set of kinet- 
ics equations [16] 



|^ = AE-A(Eo+2Bo) 

|^ = BE-B(Eo+2Ao) 

— = ^e2+2AB-E(Eo+Ao + Bo) 
The functions in Eqs. (43)-(45) are defined thus: 



(43) 

(44) 

(45) 



A(t,x,y) = kACi_oX + Xs.H^^^^ ^ 
1=2 

B(t, x,y) = keCo^y + 

i=2 



(46) 



(47) 



E(t,x,y) = ^Cj_jxV^ (48) 

i=i 

and the lack of a tilde over symbols A and B signifies that Ra and ks have been 
set to 1. (They become functions identical with those in Eqs. (53) and (54), which 
follow). The subscript ‘0’ denotes the value of a function at x =y = 1. It is not dif- 
ficult to solve the set of Eqs. (43)-(45) numerically for any Ra and ks- 

A somewhat simpler set of equations can be obtained for the reaction of 
monomers where reactivities of functional groups remain unchanged. When we 
set 

1^00 = ^10 =^01 = k,j =4 ( 49 ) 



the system of kinetics equations of the type of Eq. (38) reduces to the follow- 
ing three equations describing the rates at which molecules appear or disap- 
pear: 
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and 



dt 



^^k-l,k 

dt 



k-l 


oo oo 


^j,j^k-j,k-j-l ~ *-k,k-l 


Xs.i+2Xs-i.j 


j=l 


j=i i=i 


k-l 


- 


^^^j,j^k-j-l,k-j ~^k-l,k 




i=i 


J=1 j=l 



(50) 



(51) 



^<^k,k 

dt 



k-i k 

j=i j=i 



-c 



k,k 



Xsi+SSH 



i=i 



i=i 




(52) 



including ci_o and Cq,i since monomers can now be regarded as molecules A.. .A, 
and Equations (50)-(52) can be converted into a more manageable form 
by using the ‘untilded’ functions A and B, mentioned previously. 



A(t,x,y) = 5^Cj_j_ix'y' * 

j=i 



(53) 



B(t,x.y) = 5^Cj_i_jX' y (54) 

1=1 

to obtain the set of differential equations equivalent to Eqs. (43)-(45), but with 
all tildes removed. 

The moments of size distribution are given by 



k=l k=l k=l 



(55) 



The terms required to calculate these quantities are readily available from 
the sets of ordinary differential equations obtained by successive differentia- 
tion of Eqs. (43)-(45) with respect to x and/, followed by setting x=y=l . In fact, 
enough information is provided by differentiation of A, B, and E with respect 
to X, alone. Following the definitions of moments and averages, one can derive 
appropriate identities that eliminate the need to calculate other derivatives. 
Thus, in order to calculate number and weight average degrees of polymeriza- 
tion, it is sufficient to solve (usually numerically) the following sets of differ- 
ential equations 
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dt 



2AqBq; 1 ° “ 2AqBq5 ® — 2AqBq Eq +Aq +Bq| (56) 



and 
dA 



- EjjAq 2A^Bq5 -E^Bq 2BxAq; 



dE 



dt 



^ -2Aj(Bq + 2B^Ao E^^Ao + Bo) 



dB. 



“ ^xxAq +2Ej;A3; 2A,(3;Bq; ^ ^ ~ E B q l" 2 E^ B ^ 2 B^Aq; 

dE, 



dt 



- - Ex +2 AxxBq + 4AxBx +2BxxAq - Exx(A(, + Bg) 



(57) 



(58) 



(The subscript indicate the differentiation variables and‘0’ no differentiation 
before x andjy are set to 1.) It is not difficult to solve the sets analytically for the 
stoichiometric mixture (Ao= Bo = 1/2). The solutions together with explicit dis- 
tribution functions are collected in Table 3. 

It has been shown [16] that the functions in Table 3 lead, as they should, to the 
same simple expressions for the number and weight average degrees of poly- 
merization as those in Table 1. 



Tables. Explicit solutions of the sets of equations (Eqs. 56-58) and the distribution functions 
of polymer species for the step growth alternating copolymerization of two bi-functional 
monomers [16] with the stoichiometric monomer mixture, (- (Q)=c (0) = 4- 



A(t,x,y) = 



2x 



(2-ht)^-t^xy 



B(t,x,y) = 



2y 



(2 + t)^ -t^xy 



E(t,x,y) = - 



4txy 



(2+t)|(2 + t)^-t^xy 



Ao(t) = 



2(1 + t) 



Bo(t) = 



2(1 + t) 



Eo(t) = 



t 



(l+t)(2 + t) 



Ax(t) = 



(2 + t)^ 
8(1 -tt)^ 



Bx(t) = - 



8(H-t)2 



Ex(t) = 



t(2+t) 
4(1 -tt)^ 



Axx(t) = 



t^ (2-tt)^ 
16(1 -tt)^ 



Bxx(t) = - 



t’ 



16(14- 1)^ 



Exx (t) = 



t (2-tt) 
8(1 -tt)^ 



2t 



2(k-l) 
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6 

Random Homopolymerization of a Monomer with More Than Two 
Functional Groups 

In this section we deal with the model step growth polymerization processes 
involving monomers having more than two functional groups. The Flory-Stock- 
mayer model of homopolymerization of an /-functional monomer [18, 19] is 
included (see Table 4). This model is fundamental for understanding gelation 
processes [20]. Other processes relevant to polymer science which can be ana- 
lyzed by using the Smoluchowski equation (Eq. 3) are also briefly described. All 
are random polymerization processes, i.e., functional groups react at random 
without any substitution effect and obey the Flory assumption of the indepen- 
dence of functional group reactivity on the size of species to which they are con- 
nected. 

The Smoluchowski coagulation equation has effectively been applied to mod- 
el random irreversible homopolymerization of the monomer types presented in 
Table 4. As can be seen, in all cases the coagulation kernel has the bilinear form: 

Kj_j=A-l-B(i-l-j)-l-Cij (59) 

Spouge [21] has shown that for any kernel of this form the Smoluchowski equa- 
tion can be solved analytically, i.e., the distribution of molecular size of polymer 
molecules at any extent of reaction prior to the gel point (if any) or at least 
moments of the distribution can be expressed explicitly. 

Of particular interest are the special cases with the product kernel 



where 



Ki,i=SiSj 


(60) 


jj =ai+p 


(61) 



Table 4. The types of monomers undergoing homopolymerization and corresponding bilin- 
ear kernels in the Smoluchowski coagulation equation [21], The functional groups A react 
with each other in models 1 and 2 and with groups B only in models 3 through 5 
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and a and are arbitrary real or complex numbers. Application of the Smolu- 
chowski equation to models 1 and 2 belonging to this class have been reviewed 
several times [7, 22, 23]. Furthermore, it has been shown [24] that models with 
iCy =i +j (model 4 in Table 4) can be converted into model 2. The general solution 
for Kij=A+B{i +j) was provided by Drake and Write [25] and Treat [26] . Model 3 
was first studied by Flory [9] who used statistical arguments. It describes the sys- 
tem leading to highly branched structures which can be consider as an one-step 
method of manufacturing dendrimers. Model 5 seems not to be relevant to poly- 
mer science and the reader is referred to the original paper by Sponge [27]. 

The Smoluchowski coagulation equation with the kernel given by Eq. (60) 

dc 

= (62) 
i=l j=l 

can be transformed into the following partial differential equation 



fb>HV 




l^(9x J 








(63) 



with the function 



and 



H(t,x) = ^CjX* 



i=l 



(64) 



Hx(t) = 



I 



(65) 



Alternatively, by using the function [23] 








oo 

G(t,x) = ^Cje*‘'' 
i=l 




(66) 


one arrives at 










<?G e"^*pGV ^G^ 

dt 2 i^(9x J dx ^ 


(67) 


with 


1=1 


x=0 


(68) 
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As in previous sections, Eq. (63) is derived from Eq. (62) by multiplying it by 
x*k followed by summing over all k’s. Use is made of the identities 



Sk=Si+Sk-i-/5 



and 



k=l 




(69) 



(70) 



Eor differentiation of Eq. (61 ) with respect to x followed by putting x=l 
yields 



dH, 

dt 2 



(71) 



Hence, 



Hx=- 



(3t+ 



(a+P) 



(72) 



For p=0, the result is Hx = a. 

The zeroth moment of the distribution is now obtained from Eq. (63) by set- 
ting x=l to get 



dHp _ 1 ^2 

dt 2 



(73) 



and, upon solving 



Hn 



a + p 



-at 



a + p 



+ pt 



(74) 



or, for p=0. 



Ho=l-a ^t/2 



(75) 



Equation (63) can be solved by using the method of characteristics [6]. For the 
Flory-Stockmayer model where s; = (/- 2)i -i- 2 is the number of unreacted func- 
tionalities in an i-mer the solution is [28] 



H(t,x) 



e f_£ e 2(f-l) ft 
^ 2^ 1 + ft 



(76) 
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with the parameter ^ related to x and f through 

' ' (77) 

(The solution is obtained by solving the set of three ordinary differential equa- 
tions: dx/ df= -xHx - X\ dX/dt = - XHx; and dH/ df = - where X=3ff/ with 

the initial conditions x(0) = and H{0) = Hx= f/{l+ft), as follows from Eq. 
(72) after substituting a =/- 2 and p = 2.) The explicit form of the distribution 
of i-mers can be obtained from Eqs. (76) and (77) by the Lagrange expansion. 
The resulting distribution function for the Flory-Stockmayer model (model 1) 
is presented in Table 5. 

The aggregation processes modeled by the Smoluchowski equation with the 
kernel Kij =ij (model 2) or, more generally with K^={ij)‘^ are relevant to many 
branches of physics [29-33]. The exponent n> is a geometric factor reflecting the 
mechanism of aggregation. If the rate of aggregation is taken to be proportional 
to the surface area of clusters, then m reflects the dependence of the surface area 
on dimensionality, d. For compact clusters m=l-l/d. For aggregation or poly- 
merization processes where steric hindrances become effective or intramolecu- 
lar bonding may not be neglected, the factor co should be taken as smaller than 
1. Unfortunately, explicit solutions of the Smoluchowski equation are unknown 
for the cases with 0 < &)< 1. For the case with co=l (model 2 in Table 4) the clus- 
ters aggregate with the rate proportion al to their volume (or mass). This case is 
also similar to the random evolution of graphs, the so called Erd6s-Ren)ri process 
[34] . Both produce the same size distribution of components (aggregates or con- 
nected subgraphs). The differences between the two models have recently been 
discussed [35]. 

To solve Eq. (63) with a= 1 and p=0 {K^ =ij, model 2) one may again apply the 
method of characteristics to obtain 



H(t,x) = ^-i| h 


(78) 


x = ^exp[(l-(^)t] 


(79) 



Table 5. Moments and explicit distributions generated by the Smoluchowski coagulation equa- 
tion (prior to gel point) for models 1 and 2 in Table 4. [For both models Mi=l] 
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(Now, the set of differential equations is dx/df = -x^X +x; dx/dt = xX^ - X; and 
dHIdt = - x^X^/2 with the initial conditions x(0) = X(0) = 1; and H(0) = Hx 

= 1.) The result of the Lagrange expansion of Eqs. (78) and (79) is presented in 
Table 5. As can be seen, the second moments (as well as higher ones) diverge at 
certain t. This divergency corresponds to gelation. The post-gelation relation- 
ships (including the moment expressions more general than those in Table 5 ) are 
discussed in Sect. 8, below. 



7 

The Equivalence Between Kinetics and Simple Statistical Models 



The Smoluchowski coagulation equation deals with irreversible aggregation 
processes. There is no term in the equation describing backward reactions. Nev- 
ertheless, in 1943, Stockmayer [36] showed that for homopolymerization of an 
/-functional monomer the distribution of polymer species generated by the 
Smoluchowski equation is identical with that calculated by using statistical 
arguments. The statistical approach involves generation of the (most probable) 
distribution from the probabilities of functional groups being reacted. Accord- 
ing to Whittle [37, 38], the statistically generated distribution is that in which 
there exists both local and global equilibrium between unreacted functional 
groups and those reacted to form bonds. As Sponge [21] has shown, for a dis- 
tribution to be the equilibrium one the following combinatorial identity must 
hold: 



(k- 



_l)Zk 

k! 



k-l 






W; W 



k-i 



i! (k-i)! 



(80) 



In Eq. (80) Wk is the number of ways of assembling a k-mer from distinguish- 
able units. The identity simply states that the rate of assembling a k-mer from 

components i and k - i (possible in I | = — — ways) equals the rate of break- 

{ij i!(k-i)! 

ing one of A: - 1 bonds in the A-mer. For this identity to be valid, the kernel Xy 
must be equal to the number of ways the components can be linked together. All 
kernels in Table 4 meet this requirement. In model 1 (Flory-Stockmayer), it is the 
number of ways the functional groups can react; in model 2, it is the number of 
units in the reacting components between which a link may be formed, etc. 
Sponge has also generalized this requirement by proving [39] that the identity at 
Eq. (80) is valid if and only ifKij has the bilinear form given by Eq. (59). The case 
with a constant value of iCy, dealt with in Sect. 3, also belongs to the same class 
of kernels. 

The existence of Wk related to the distribution of aggregate sizes and ful- 
filling the identity at Eq. (80) implies that the aggregation process leads to an 
equilibrium distribution. For the Flory-Stockmayer model (model 1 in Table 
4) the number of assembling an acyclic A-mer from distinguishable units is 
[40] 
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f'^(fk-k)! 

'^''“[(f-2)k+2]! (81) 

while for Kij=ij 



Wk 



= k '^“2 



(82) 



It is now easy to see that if these Wk fulfill Eq. (80), the explicit distributions c/t 
in Table 5 necessarily fulfill the relation 

k-l 

(k-l)Ck=i^Ki^_jCiCk_i (83) 



The same recurrence equation is obtained by solving successive terms of the 
Smoluchowski equation with 2 s,c; = Hx equal to //(I +ft) and 1 for models 1 
and 2, respectively. ’ ' 

In the past, the equivalence between the size distribution generated by the 
Smoluchowski equation and simple statistical methods [9, 12, 40-42] was a 
source of some confusion. The Spouge proof and the numerical results obtained 
for the kinetics models with more complex aggregation physics, e.g., with a pres- 
ence of substitution effects [43, 44], revealed the non-equivalence of kinetics and 
statistical models of polymerization processes. More elaborated statistical mod- 
els, however, with the complete analysis made repeatedly at small time intervals 
have been shown to produce polymer size distributions equivalent to those gen- 
erated kinetically [45]. Recently, Faliagas [46] has demonstrated that the kinet- 
ics and statistical models which are both the mean-field models can be consid- 
ered as special cases of a general stochastic Markov process. 

8 

Gelling Systems with Substitution Effect 

Modeling of a polymerization system becomes more and more complicated as 
the number of different kinds of monomers and reactive groups increases. To the 
complexity add substitution effects that are quite common for monomers of 
practical importance. Furthermore, the substitution effects may be considered 
as a method of relaxing the Flory postulate of the independence of reactivity of 
the size of molecules. They model conversion dependence of the reactivity in the 
system. 

For simplicity, we start with a single tri-functional monomer RA^, the func- 
tional groups of which react with each other with the first shell substitution 
effect. It is now convenient to write down six types of reactions [44] , one for each 
pair of reacting units of substitution degree 0, 1, and 2. They are presented in 
Table 6 together with the appropriate rate constants. Similarly as in previous sec- 
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Table 6. The types of reactions in model polymerization of RA3 monomer with functional 
groups reacting with the first shell substitution effect. The subscripts at the symbols of 
reagents correspond to the codes of molecules (y-mers, see the text) and those at the rate con- 
stants indicate the degrees of substitution of reacting units. The left and right column corre- 
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tions, the constant kafj (a, j5=0, 1,2) describes reactivity of functional groups in 
a pair of reacting units of substitution degree a and p. All unreacted groups on 
a unit have the same reactivity which depends solely on its substitution degree. 

Table 6 presents the reactants and products of reactions. The pairs of num- 
bers in subscripts are codes for molecules. The first code number for an (y)-mer 
is the number of units in polymer molecule with exactly one reacted functional 
group (substitution degree 1) while the second is the number of units with two 
reacted groups (substitution degree 2). The symbol ?o,o stands for the monomer. 
Note that there is no need to specify the number, say k, of units with all three 
groups reacted; they do not enter reactions anymore and, as long as the mole- 
cule has no cycles, the two parameters provide all information on its size n. The 
following two relations hold: the obvious one n=i +j + k and 2{n - l)=i + 2j + 
3k. The second relationship sums up all reacted groups that form the n-1 bonds 
in an acyclic molecule (two groups per each bond). The parameter k may now be 
eliminated to yield 



n = 2i + j-2 (84) 

Kinetic equations describing the rates of formation and further reactions of 
(y )-mers are more complicated. They contain up to ten different terms. Consid- 
erable simplification and relatively easy access to average parameters of reacting 
systems (such as average degrees of polymerization) can be achieved by apply- 
ing the assumption on the additivity of activation energies similar to that adopt- 
ed in Sect. 3. The rate constants can then be split into the product form: 
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(a,/?=0, 1>2) where /c^stands for the substitution independent part of the constant. 
The number of rate constants can now be reduced to just two [44] : 




With 9k*kg absorbed into time units, the relative rate constant with which a 
monomer reactive group reacts is simply unity, while the relative rate constant 
for mono- and disubstituted units is jci and K 2 , respectively. 

Note that for the random case, equivalent to the Flory-Stockmayer model 
(model 1 in Table 4), jci=2/3 and K2=1/3. 

For the function defined as 



00 00 

H(t,x,y) = ci(t)+ ^^Ci^ (t)[K-ixf [jC2y]^ (87) 

i=2 i=o 

(ci(t) = CQfl{t)) the set of kinetics equations for the reactions in Table 6) can be 
transformed into the single partial differential Smoluchowski-like equation in 
the form 



<9H ,1 


(9H <9h’ 


2 

( TT TT \ 


(9H (9h" 


QjI 

II 

to |l- 


K^XC^ + K2Y — + — 
dx dy 


“(Ci+Hx+Hyj 


ci+^ix— +?f2y^ 
dx dy 



(88) 



where 



Hx=Hx(t) = ||| 



\=VK^,y=llK2 i=2 j=0 



(89) 



(0 = 1^1 



x=llK2,y=llK2 i=2 j=0 



(90) 



Equations (89) and (90) define functions of time obtained by differentiation 
followed by substitution of x andy by 1 /ki and 1 IK 2 , respectively, so that the dum- 
my variable terms cancel in the derivative. i j 

As in the previous sections, Eq. (88) is obtained by multiplying by j 

all possible terms that can be written down for formation, and disappearance, of 
an i,j-mer followed by summing up the result over all i and j. 

In order to verify Eq. (88) note that the term corresponding to the rate of for- 
mation of (i,j)-mer by reaction of two monosubs tituted units (reaction 4 in Table 
6) is 



i-i j-i 
r=2 s=0 



( 91 ) 
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(in parantheses are the numbers of functional groups available for reaction). 

This term yields, after multiplying by and summing up over i 

and j, 



1 

2 






<9x 



(92) 



The disappearance term, say, the one corresponding to reaction 5 in Table 6 



reads 


00 00 






r=2 s=l 


(93) 


and turns into 


of Eq. (88). 


“X'lX^— Hy 
dx ' 


(94) 



It is not difficult to extend the reasoning to the general case of polymer- 
ization with the first shell substitution effect approximation. By defining the 
function 



H(t,Xi,..,Xf_i) = Ci(t) + X S •• X ]' ['^2^2 1 (95) 

i=2 j=0 z=0 

where the Kr stand for rate constants relative to that of monomer 






(96) 



one arrives at the general Smoluchowski-like equation in the form 

n2 



(9H 1 
di “2 



<9H <9H ,9H 

JClXiCi -h JC2X2 — + + K-f_iXf_i -h 

<yX^_2 (y X^_2 



-(ci -h Hi -h .. -h Hf_2 + Hf_i) 



<9H (9H (9H 

Cl + X’lXi -I- .. -I- X-f_2Xf _2 — -H K-f _jXf _1 

<yXj O' X^_2 C'X£_j^ 



(97) 

with Hr(t) defined by relations analogous to Eqs. (89) and (90). An equation 
equivalent to Eq. (97) was first derived by Kuchanov [43,47]. If there is no sub- 
stitution effect we have 



jc^ =(f-r)/f ;(r = l,2,..,f-l) 



( 98 ) 



and the system becomes the general Flory-Stockmayer one. 
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An interesting new model is provided by Eq. (97) with k^= 1 for all r. In the lim- 
it of/^ oo it becomes the RAoo model, fully equivalent to that described by Eq. 
(3) with Kij=ij. At bounded /it describes the evolution of /-trees, i.e., a process 
similar to that dealt with by Erdos and Renyi [20, 34], but with substitution 
degree restriction imposed on the vertices (units) and cycle formation disal- 
lowed. 

Explicit solutions of Eqs. (88) or (97) are not known, even for the Flory-Stock- 
mayer model. The concentrations of individual species Cy(f) for that model are 
available, but only by recursive solution of appropriate kinetics equations (before 
they are summed up to yield Eqs. (88) or (97)). The function ci(f) -i- Hx(t) + Hy(t) 
which is needed to solve the equations is 

c,-tH„-tH„=— (99) 

y t-i-3 

Physically, it is the fraction of functional groups available to reaction (and 
hence equivalent to Hx in Eqs. (65) and (72)). Equation (99) reveals that for this 
version of the FS model the time units are „nine times longer" than those from 
direct application of the original Smoluchowski equation. 

Equation (88) was used to demonstrate [44] the differences in distributions of 
molecular sizes in pre-gel stages of an RA^ polymerization with substitution 
effects as calculated according to a statistical and a kinetics model. The moments 
of distribution and the gel points were calculated by the numerical solution of a 
few ordinary differential equations that were derived from Eq. (88) and compared 
with analogous quantities calculated from a statistical model. 

To illustrate the application of Eq. (88) we will derive equations needed to 
calculate the second moment of molecular size distribution. The r-th moment of 
the size distribution which, taking into account Eq. (84), is defined thus: 



M, (t) = Cl (t) + ^ ^ (2i + j - 2)'' C; j (t) (100) 

i=2 i=o 

can be expressed as linear combinations of successive derivatives of H with 
respect to x and/or y of at most order r, calculated at x=\!k\ andy=l/K 2 . 

By raising the term right after the sum signs in Eq. (100) to the second power 
0 and regrouping the result we get 



^2 =Ci-H^^[4i(i-l)-l-4ij-l-j(j-l)-4i-3j-h4]cij = 
i i 




XX 



H Hxy H yHyv Hx Hy -h 4Hq -3Cj 

k \ 'f'l K'2 



(101) 



The subscripts of H indicate variables with respect to which differentiation 
has to be made in Eq. (87) (0 for no differentiation) before the dummy variables 
are set equal to 1 /ki and I/K 2 as in Eqs. (89) and (90). The required derivatives 
are obtained from the Smoluchowski-like equation (88) by 
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i. setting x=l/ K\ and j=1/ K 2 to get (upper dot denotes time derivative) 



Ho=-l(q+H,+H^) 


(102) 


ii. setting x=0 andy=0 to get 




Cl =-Ci(ci+H,+Hy) 


(103) 


iii. differentiating appropriately Eq. (88) with respect to dummy variables x and 
y and then setting x=Hk\ andy=l/jC 2 to get 


Hx=-ki(H,-Ci)(ci+H,+H^) 


(104a) 


Hy=-k2(H^-H,)(ci+H,+H^) 


(104b) 


Hxx = (kiCi -H ) -2kiH,o,(ci - hH,; - iHy) 


(104c) 


Hxy =(kiCi -i-Hxx +Hxy)(k2Hx +Hxy -i-Hyy)-[(ki -i-k2)Hxy -k2Hxx] 


|(ci+H,+Hy) 

(104d) 


Hyy = (k2H, + H,^ + Hyy )' - 2k2(Hyy - H,^ )(ci + H, + Hy ) 


(104e) 


Equations (102)-(104) are the Riccati differential equations that have no solu- 
tions in quadratures for arbitrary k\ and K2. They can, however, be easily solved 
numerically to obtain the desired second moment of the distribution (and hence 
the weight average degree of polymerization) as a function of time, according to 



Eq. (101). 

Only for the FS model (ki= 2/3, K2=1/3) have explicit solutions of Eqs. 
(102)-(104) been found. The functions Hxx, Hxy, Hyy have, as expected, a singu- 
larity at f=3 of the form (3 - t) h The singularity corresponds to the gel point. At 
fc=3, where tc is the time where the gel point is, these functions simultaneously 
diverge. No modifications of Eq. (88) that would be valid beyond this point are 
known. 

For arbitrary ki and K 2 no explicit expression for the gel point exists. In sta- 
tistical models, the gel point can be found as the time (or conversion) when the 
fraction of fully reacted units reaches one third of the fraction of monosubsti- 
tuted units. Simply, when there is not enough terminating units to saturate all 
branches extending from those with all three groups reacted, the infinite (gel) 
molecule necessarily forms. Note that ci, HxIki, and Hy/K 2 are just the fractions 
of units with 0, 1 and 2 reacted functional groups, respectively. The fraction of 
those with all three groups reacted is 1 - ci - H^/ki - Hy/k 2 . It is interesting that. 
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except for the FS case, the simple condition for gelation deduced from statistical 
models does not hold for the kinetics one (Eq. 88). The functions Hxx> Hxy, Hyy 
are found to diverge before that condition was met. 

Modeling with the Smoluchowski-like equation generalized to take into 
account FSSE is not limited to the simple RAf polymerization. A kinetics 
approach similar to that described in this section have been used to study 
crosslinking reactions of epoxy resins with components introduced into the sys- 
tem at different times [17]. Kinetic equations analogous to Eq. (101) have been 
derived [48] for an RA 2 + R’B^ system as well as for systems containing 3-func- 
tional monomers having functional groups of intrinsically different reactivities 
[49]. 

The effect of changes in reactivity of monomers due to substitution on the 
process of polymerization, mostly on the position of gel point was studied by 
using Monte Carlo calculations applying essentially the same approach as that 
presented in this section [50-52]. 



9 

Post-Gelation Relationships 

In principle, the Smoluchowski coagulation equation and its modifications pro- 
vide the correct distributions of aggregate sizes at positive values of the para- 
meter t. For some kernels the solutions are limited to the range of 0 < f < tc. At tc 
the second and higher moments of distribution diverge. In this section the 
changes in distribution for t> tc are discussed. 

It can be seen from Table 5 that for models 1 and 2 the gelation time is at 

C= — — andC=l (105) 

" f(f-2) 

respectively. (Simply the denominator in Mj is zero.) If we look at the explicit 
distributions in Table 5, we find that the concentrations for every k> 1 initially 
increase with time and then peaks and decays. The value of tc is the limiting time 
where the molecule of size reaches its maximum concentration 

Ic ~ limtmax(^) (106) 

k — >00 



( fmax is obtained by setting dct/dt=0. ) By that time, the concentrations of all finite 
(sol) molecules have already passed their peak values and decrease because of 
the growth of the largest species. The largest molecule can then be identified with 
the gel. Since k^ the expression for its concentration cannot be evaluated. 

At the critical time when gel appears, tcy the first moment of the distribution, 
which is constant for t<tc and reflects mass conservation in the system, starts 
to decay, corresponding to the loss of units to the gel molecule. The initial rate 
of the decay is easy to deduce. Since, for model 1, direct differentiation of the 
explicit form of concentration (Table 5) yields 
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dck k-l-[(f-2)k + 2]ft] 
dt t(l-hft) 

the rate of decay of the first moment is 

dMi f^(f-2) 

dt , — dt f-1 

tf k=l 



(107) 



(108) 



(More precisely the upper limit of summation in Eq. (108) should be set to 
some value L and allow to grow to infinity only after tc=llf{f- 2) is substituted.) 
Similar reasoning for model 2 leads to the initial decay rate equal [7] to -1. 

The Smoluchowski equation was thought of as failing beyond the gel point 
[22]. Dusek [53] and Ziff [54] have demonstrated this not to be true for the Flo- 
ry-Stockmayer model. At the same time Leyvraz and Tschudi [32] presented an 
elegant and more general derivation of the solutions of the Smoluchowski equa- 
tion both prior to and beyond the gel point. 

To keep the mathematics as simple as possible, here we sketch the arguments 
of Ziff [54] and Ziff and Stell [28], who carefully analyzed the post-gelation 
behavior of model 1. Ziff and Stell distinguish three versions of its post-gel 
behavior. 

The Flory model is the version where the equivalence between kinetics and 
statistical descriptions is extended to the post-gel stage of polymerization. Con- 
sequently, the functional groups are assumed to continue to react at random with 
no distinction on whether they belong to sol molecules or to gel. To analyze this 
version one can use the explicit form of function H. As usual, the moments are 
available through successive derivatives of H (Eq. 76) with respect to x calculat- 
ed at x=l. We may rewrite Eq. (77) in the form 

i=c+ft(c-c'“‘) 

where ^ stands for ^ at x=l. 

Equation (109) has two roots in the half-open range (0,1]. The trivial root, ^ = 
1, and the root given by the equation 

( 110 ) 

obtained from Eq. (109) by removing the factor 1 - f . The second root, however, 
becomes accessible only past the gel point. Bearing that in mind we can now 
express the moments of the distribution by using Eq. (76) with the lowest posi- 
tive root of Eq. (110). The equations will then apply to the whole system prior the 
gel point and to the sol fraction only beyond gelation. Thus, 



Mo=H(t,l) = ^ '--(f-2)C' 



( 111 ) 
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(Substitute ^ from Eq. (109) into the second term of Eq. (76)), 



and 






-3(f-2)Mi-2Mo 


3 




X=1 



2(l+ft)C^-C^~^ 

l+ft-ft(f-l)f ^-2 



( 112 ) 



(113) 



where the subscript ^ (^^ ) indicates that the functions in Eqs. (76) and (77) are 
differentiated once (twice) with respect to 

The concentration of functional groups in the sol is 

H (114) 

1 + ft 

etc. 

The reader familiar with the cascade theory will notice that the root 0 < ^ < 1 
is related to the extinction probability, v, i.e., the probability of a unit chosen at 
random to belong to a sol molecule [14, 55]. This probability is [53] 



y = C 



f-i 



(115) 



The second version of the post-gel analysis is equivalent to that of Stockmay- 
er [19]. He assumed that the distribution of sizes within sol past the gel point 
remains unchanged and the ratios of concentrations can still be deduced from 
those given in Table 5. Ziff and Stell [28] concluded that this assumption is equiv- 
alent to a disregard of reactions between gel and sol. Consequently, the concen- 
tration of each k-mer as well as that of the sol fraction decays according to the 
equation 

,, , ,(f-i)b 

Ck(t) = Ck(tc)^— ^ (116) 

Ziff and Stell [28] themselves proposed a third version of the post gel analy- 
sis which differed from Elory’s one by disregarding the reactions taking place 
within gel molecule (no intramolecular reactions both in sol and gel). 

The analysis for model 2 of Table 4 (Kij=ij) using the Elory model leads to the 
following moment expressions: 

c 

1-Ct 



Mo = C-Ct/2; Mi=C; M2 = 



(117) 
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where ^ is the real root of the equation 

(1-Ot-tlOgf =0 

Beyond the gel point (t > 1), the nontrivial root 1 1) should be taken into 
account. 

Of interest also are the results of applications of the Smoluchowski equation 
for systems with more complex aggregation physics than that provided by bilin- 
ear kernels. Leyvraz and Tschudi [32] conjectured that for the kernel iCy=(r)')® 
gelation occurs only when cu > 1/2. The post-gelation behavior of a general sys- 
tem with a multiplicative kernel given by Eq. (60) has been analyzed by van Don- 
gen et al. [56]. By assuming that the distribution past the gel point could be 
expressed through that at t=tc, thus 



with 



Ck(t) 



Ck(tc) 

H-b(t-t^,) 




(119) 

( 120 ) 



they were able to obtain the recurrence relation very similar to Eq. (83) 



where 



i+j=k 



(121) 



(Sk-l)n,,=i ^SiSpinj ^22) 

i+j=k 

Post-gel distribution functions or their Laplace transforms were derived for 
several special cases including Sk=k^- Unfortunately, no analytical method of 
predicting tc has been found to exist. A numerically derived relationship between 
tc and CO have recently been published [57]. With time units scaled so that tc=l 
for model 1 it reads 



1 

8.902 10"^ -1.3034(y-k2.2144ffl^ 



(123) 



10 

Attempts at Taking Into Account Cyclization Reactions 

The Smoluchowski coagulation equation describes the rate of formation of 
acyclic aggregates. Only then it describes the evolution of a Markovian distribu- 
tion [34]. Strictly speaking the Smoluchowski equation simply disregards any 
cycle formation. For polymers this is true for models with high functionality 
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or BRAo<,) where the ‘wastage’ of functionalities in cycle formation can be 
considered negligibly small. There have been attempts to include cyclization 
reactions into polymerization models, but, except perhaps the linear polymer- 
ization models, they are usually very crude or just not realistic. 

The set of kinetics equations for linear polycondensation can be written as [58] 



P -tP 

J.J I ij 

Pi- 



ck. . 






■+) 




(124) 



where the first equation is the same as that in Eq. (7) and the second one 
describes the cyclization reactions. Q stands for a cycle of size i. The rate con- 
stant in the second reaction is usually made dependent on the size of molecule. 
The standard assumption is [59] that it has the form 



k(^>~r 



■ 3/2 



(125) 



This assumption is equivalent to considering the polymer molecule as a 
Gaussian chain. For a Gaussian chain the probability of the two ends colliding in 
three-dimensional space is proportional to its length to the power -3/2. For the 
Kuhn (or freely-jointed chain) model the same assumption maybe taken for suf- 
ficiently long chains [60]. For linear polymers in good solvents, no similar sim- 
ple assumption can be adopted. To study cyclization one has to resort to more 
sophisticated mathematical treatments (see, e.g. [61]). 

Since in irreversible linear polymerization the cyclic molecules once formed 
do not react further, the kinetics equations can be written thus (here Cj is the 
concentration of cycles): 



dCk 

dt 

dCk 

dt 



k-l 




-Ck£ci-kfc^ 



i=l 



(126) 



A kinetics approach to model polymerizations of RAo^ and RAf monomers 
including cyclization was used by Lu and Bak [62] . They formulated and solved 
four systems of equations, in which cyclization terms were added to the Smolu- 
chowski equation: 

dc °° 

(Jt = ^27) 

with the parameter A either a constant or linear function of the size of k-mers 
(taken as Xsk with Sk=ock + P). The disadvantage of the model is that it means that 
the species that undergo cyclization are automatically removed from considera- 
tions. An explicit gel time was derived for the kernel Kij=ij and size independent 
A. It reads (for the monodisperse initial conditions): 
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tc - — -log(l-A)-l+— + — + .. + ^ — --h.. (128) 

A 2 3 1 “h 1 

Equation (128) can be derived by converting Eq. (127) into the partial differ- 
ential equation of the same type as Eq. (63) and solving the ordinary differential 
equation for the second moment 

(129) 

dt 

Note, however, that for constant A, it is the monomer which at small t contributes 
to cyclization the most, which is physically unrealistic. 

It is not very difficult to extend formally the treatment presented in Sect. 8, 
namely the Smoluchowski-like equation (Eq.88),to model, besides the substitu- 
tion effect, the ability of functional groups to react intramolecularly. For the sim- 
plest case of RA-i homopolymerization, a crude method [63] is to code the mol- 
ecules with four indices two of which count the units with two or three reacted 
functional groups that are engaged in cycles. The Smoluchowski-like equation 
reads then 



(9H 1 
di ~ 2 



<9H 

K1UH0+K2X— + Y—+— 

du dv d X 



i2 



(Ho + H,+H^ + H,) 
[(k-2v)^-(k:iU)^ Ai 



<9H 






Hq -t-STjU^:^ h/C, V^:^ l-K’-, X 

d^H 






d u^ 



-i-2[yK:2 v-/CiUK'2x]A2 



d^U 

dud X 






(130) 

d^K 

dx^ 



and describe the time evolution of the generating function 



H(t,u,v,x,y) = Ho(t)-t IIII [k-2vJ [K-2x]‘[y]“ 

i j 1 m 



where 



(131) 



H„(t) 



dn 

da 



U=l/Vi,V=l/K’2,X=l/K'2.y=l 



(a=u,v,x,y) 



(132) 



and u, V, x andy are the auxiliary variables labelling the types of units (see Fig. 
2). As many as three parameters A, (z=l,2,3) model the extent of intramolecular 
reaction, each for a pair of differently substituted units (they may be considered 
functions of time,size of molecules, dilution of the system, etc.). In one approach 
they were assumed constant and estimated by using graph theoretical methods 
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U - ^ V X 



Fig. 2. A single molecule formed in RA3 polymerization. The molecule contains various types 
of units. The labels of units are the same as the dummy variables in Eq. ( 130 ): u - monosub- 
stituted unit with two unreacted groups (solid lines), v - disubstituted unit, % - disubstituted 
unit engaged in cycle, y - trisubstituted unit engaged in cycle 

[64] to be in the ratio A3M2 = 1>5 and A2M1 = 2. Unfortunately, Eq. (130) has not 
been solved or examined numerically, as yet. 
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Other Remarks on the Smoluchowski Coagulation Equation 

The continuous form of the Smoluchowski coagulation equation 

— ^^ = AjK(x-u,u)c(x-u,t)c(u,t)du-c(x,t)jK(x,u)c(u,t)du (133) 

0 0 



was first given by Muller [65] . Here, x and u are real numbers standing for aggre- 
gate sizes and K{x,u) is a collision frequency function. Another method of mod- 
eling aggregation processes is the time-dependent Markovian process known as 
the Lushnikov process [66]. The analysis by Faliagas [46] is in fact an elegant 
application of LushnikoVs model to /-functional polymerization. 

Considerable attention has been devoted to aggregation processes, the physics 
of which is well described by the Smoluchowski equation with so-called homo- 
geneous kernels (i.e., with Kai,aj=oc^Kij). Here, parameter A is the degree of 
homogeneity. To describe the size distribution of aggregates in a finite system, 
functions have been proposed [67] in the form 



Ck(t) 



K 



fi— 

s(t) 



7 



(134) 



where /( ) is a scaling function, or, more specifically, as 



Ck(t) 



[s(.)f 




f 

exp 

V 




(135) 
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where s{t) is a scaling parameter usually identified with the z-average size of 
aggregates, which behaves according to 

t^ A <1 

s(t)~- e‘; A = 1 (136) 

(t„-tf';A>l 

The scaling function (Eq. 135) holds for both time t and aggregate size k suf- 
ficiently large [68] . For particular kernels, the exponents, z and z’, have been cal- 
culated. They are summarized in Table 7. 

As mentioned previously, the most relevant to polymer science seems to be 
the product homogeneous kernel ffy=(z)')®, (the homogeneity degree X=2co). 
According to van Dongen and Ernst [69], exponents in the scaling relation giv- 
en by Eq. (135) are 

T = 0) + - (137) 



They also pointed out that the scaling function at the gel point reduces to (cf. 
[31]). 

Ck(tc)~k ^ ( 135 ) 

The weight average degree of polymerization changes in the vicinity of gel 
point according to [70] 

Mw~(tcr-tn (140) 



3 — 2CO 
2CO-1 



For negative co, the coagulation leads after long time to a log normal distrib- 
ution [71]. 



Table ;.TIie ex|ioiwiils in llie liiiiclioii fEq. I.f 51dacribijif llie j5yiii|ilolk: disiribiilioii ol 
ill j |iioccm iiiodclf^ i>.'illi llir5iiioliiclim>'TkicojfiiLjlioii A]iijlioii 
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Critical times at gelation and exponents 0 in Eq. (135) for several values of co 
in the range 0.55-1 have recently been calculated numerically by Krivitsky [72]. 

On the other hand, in the Erdos-Renyi process, the size of a „typical“ largest 
component in a system comprising N units has been proved to be proportional 
to log N up to the gel point. Beyond the gel point it (i.e., gel molecule) has a size 
proportional to N itself [20, 34]. 

An interesting modification of the Smoluchowski equation leading to the so 
called generalized Smoluchowski equation [73] can be obtained by adding the 
terms describing the rates at which clusters or polymer molecules undergo frag- 
mentation: 

k-1 oo k— 1 oo 

^ K; k_iCk + ^ Ki kCj+k q^2) 
i=l i=l i=l i=l 

is the rate constant of the splitting of an (z+jl-mer into i- and j-mers. The third 
sum is the rate at which k-mer dissociate into smaller species, while the fourth is 
the rate at which it is formed from bigger ones. In general, the fragmentation 
terms do not complicate seriously the mathematical treatment of combined 
aggregation-fragmentation models [74-76]. 
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